function plos_graf
x=linspace(-pi/2,pi/2,20);
y=x;
[X,Y]=meshgrid(x,y);
Z=cos(X.^2+Y.^2);
subplot(2,2,1)
mesh(X,Y,Z)
subplot(2,2,2)
surf(X,Y,Z)
subplot(2,2,3)
contour3(X,Y,Z)
subplot(2,2,4)
contour(X,Y,Z)
colorbar
end

