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