clear;
 clf;
 
 deff('[y] = pression(x,t)', 'y = 1+sin(t)*sin(x)')
 deff('[y] = vitesse(x,t)', 'y = -sin(t+%pi/2)*cos(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));
 p=2*ones(nx,4); // pour fixer les extremes de l'echelle de couleur,
 p(:,4)=zeros(nx,1); // les niveaux 0 et 2 sont affichés en permanence
 vx(:,1) = 0:%pi/9:%pi; // position des vecteurs vitesse
 nx=max(size(vx));
 vy1 = 0.75*ones(2,nx);
 vy2 = 0.25*ones(2,nx);
 vy3 = -vy2;
 vy4 = -vy1; // 4 lignes de flèches
 xset('colormap', graycolormap(256)); // niveaux de gris
 
 for i=0:9
    for j=0:9
       t=(i*10+j)*0.02*%pi; // temps
       clf;
       p(:,2)=feval(x,t,pression);
       p(:,3)=p(:,2);
       vx(:,2)=vx(:,1)+0.35*feval(vx(:,1),t,vitesse);
       scf(0);
       subplot(2,1,1)
       grayplot(x,y,p,rect=[0,-1,%pi,1],axesflag=0);
       subplot(2,1,2)
       plot2d(vx(:,1),vy1(1,:),...
          rect=[0,-1,%pi,1],style=0,axesflag=0) // pour l'echelle
       xpoly([0,%pi],[0.95,0.95]);
       xpoly([0,%pi],[-0.95,-0.95]);
       xarrows(vx',vy1,arstyle=1); // vecteurs vitesse
       xarrows(vx',vy2,arstyle=1);
       xarrows(vx',vy3,arstyle=1);
       xarrows(vx',vy4,arstyle=1);
       nom='test\test'+code2str(i)+code2str(j)+'.gif'; // nom du fichier
       xs2gif(0,nom,1); // sauvegarde
    end;
 end;