x=linspace(-1,1,300);
y=linspace(-1,1,300);
[X,Y]=meshgrid(x,y);
Z=X.^2+Y.^2;
Z1=2*(X.^2+Y.^2);
Z3=linspace(0,4,300);
Z2=Z3';
for i=1:299
    Z2=[Z2,Z3'];
end
size(Z)
Y1=X.^2;
mesh(X,Y,Z);
hold on
mesh(X,Y,Z1);
hold on
mesh(X,X,Z2)
hold on
mesh(X,Y1,Z2)
%xis on