clear;
 clf;
 
 deff('[y] = note1(x,t)', 'y = sin(t)*sin(x)')
 deff('[y] = note2(x,t)', 'y = 2/3*sin(t)*sin(x) + 1/3*sin(2*t)*sin(2*x)')
 deff('[y] = note3(x,t)', 'y = 0.57*sin(t)*sin(x) + 0.29*sin(2*t)*sin(2*x) + 0.14*sin(3*t)*sin(3*x)')
 
 x=0:0.1:%pi;
 
 for i=0:9
    for j=0:9
       t=(i*10+j)*0.02*%pi;
       clf;
       y1=feval(x,t,note1);
       y2=feval(x,t,note2);
       y3=feval(x,t,note3);
       scf(0)
       subplot(3,1,1)
       plot2d(x,y1,rect=[0,-1,%pi,1],axesflag=0)
       subplot(3,1,2)
       plot2d(x,y2,rect=[0,-1,%pi,1],axesflag=0)
       subplot(3,1,3)
       plot2d(x,y3,rect=[0,-1,%pi,1],axesflag=0)  
       nom='test'+code2str(i)+code2str(j)+'.gif';
       xs2gif(0,nom,1)
    end;
 end;