>> t=0; y=1:10; x=y;[x,y]=meshgrid(x,y); for n=1:100; t1=3.*sin(x).^2+cos(2.*y)+t;t=t1; end, grid, m