x=linspace(0,100,1000);
y=linspace(0,100,1000);
[plegmax,plegmay]=meshgrid(x,y);
z=plegmax+plegmax*plegmay+4*plegmay;
mesh(plegmax,plegmay,z)
view(0,90)