x=linspace(0,100,1000);
y=linspace(0,100,1000);
[plx,ply]=meshgrid(x,y);
z=plx+plx*ply+4*ply;
mesh(plx,ply,z)