n=12;
m=moviein(n);
t=0:2*pi/n:4*pi;
x=0:pi/12:4*pi;
nj=length(x);
for i=1:n
   for j=1:nj
      y(j)=sin(x(j)-t(i));
   end
   plot(x,y);
   axis([0,4*pi,-1.5,1.5])
   m(:,i)=getframe;
end
movie(m,20)