clear all

x1 = 0:pi/50:2*pi;
x2 = 0:pi/50:3*pi;

for i=1:length(x1)
    for j=1:length(x2)
        
        y(i,j) = sin(x1(i))*cos(x2(j));
        
    end
end

t = 0:pi/20:4*pi;
z = sin(t);

colormap(jet)

%contour(x2,x1,y,20)

for k=1:length(t)

   mesh(x2,x1,y*z(k))
   % surf(x2,x1,y)
   % shading interp
   % hidden off
   xlabel('x_2')
   ylabel('x_1')
   zlabel('Y(x1,x2)')
   axis([0 3*pi 0 2*pi -1 1])
   view(30,45)
   pause(0.1)

end