[x,y] = meshgrid(1:200);z = 200*exp(-((x-100).^2+(y-100).^2)/2);mesh(x,y,z);
[x,y] = meshgrid(1:200);h = 200*exp(-((x-100).^2+(y-100).^2)/2);surf(x,y,z);
