clear;
clf;

deff('[y] = pression1(x,t)', 'y = 1+sin(t)*sin(x)')
deff('[y] = pression2(x,t)', 'y = 1+sin(2*t)*sin(2*x)')
deff('[y] = pression3(x,t)', 'y = 1+sin(3*t)*sin(3*x)')

x=0:0.05:%pi; // une demi-onde par tuyau
y=[-1,-0.9,0.9,1]; // pour la carte de gris
nx=max(size(x));
p1=2*ones(nx,4); // pour fixer les extremes de l'echelle de couleur,
p1(:,4)=zeros(nx,1); // les niveaux 0 et 2 sont affichés en permanence
p2 = p1 ; p3 = p1 ;

xset('colormap', graycolormap(256)); // niveaux de gris

for i=0:9
   for j=0:9
      t=(i*10+j)*0.02*%pi; // temps
      clf;
      p1(:,2)=feval(x,t,pression1);
      p1(:,3)=p1(:,2);
      p2(:,2)=feval(x,t,pression2);
      p2(:,3)=p2(:,2);
      p3(:,2)=feval(x,t,pression3);
      p3(:,3)=p3(:,2);
      scf(0);
      subplot(3,1,1)
      grayplot(x,y,p1,rect=[0,-1,%pi,1],axesflag=0);
      subplot(3,1,2)
      grayplot(x,y,p2,rect=[0,-1,%pi,1],axesflag=0);
      subplot(3,1,3)
      grayplot(x,y,p3,rect=[0,-1,%pi,1],axesflag=0);
      nom='test\test'+code2str(i)+code2str(j)+'.gif'; // nom du fichier
      xs2gif(0,nom,1); // sauvegarde
   end;
end;