u = linspace(-5,5,20);
v = linspace(-5,5,20);
[u,v]=meshgrid(u,v);
X = u.^2;
Y = v.^2;
Z = u+2.*v; 
surf(X,Y,Z)
axis square