frames=60;

for f=1:frames
    navier_stokes_visualization(f-1);
    hold off;
    M(f)=getframe;
end