[X,Y]=meshgrid(0:0.2:2*pi,0:0.2:2*pi);
for k=1:15
surf(X,Y,0.8*k.*sin(X).*cos(Y)) ;
axis([0 2*pi 0 2*pi -20 20]);
M(k)=getframe;
end;

movie(M,2,3) %video jede 2x, 3 FPS frames per second, 3 snímky za sekundu

filename = 'muj3d.gif';
for k=1:15
im = frame2im(M(k)); %RGB formát
[imind,cm] = rgb2ind(im,256); % indexed image - 256 barev 0-255
  if k == 1 
          imwrite(imind,cm,filename,'gif', 'Loopcount',inf); % opakuje se
      else 
          imwrite(imind,cm,filename,'gif','WriteMode','append'); % přidá se další
      end 
end;

v = VideoWriter('muj3d.avi');
open(v);
writeVideo(v,M);
close(v);