[X,Y] = meshgrid(1:0.5:10,1:20);
Z = X*X+Y*Y;
surf(X,Y,Z)