[K,L]=meshgrid(0:10.0:300);
P=1.01.*L.^0.75.*K.^0.25;
surf(L,K,P)
colormap(spring)
shading interp
box on
view(60,20)
figure()
[c,h]=contour(L,K,P);
clabel(c,h);

[x,y]=meshgrid(-1.0*pi:0.1:pi);
z=(x.^2+3*y.^2).*exp(-1.0.*x.^2-y.^2);
surf(x,y,z)
colormap(winter)
shading interp
box on
view(60,20)
figure()
[c,h]=contour(x,y,z);
clabel(c,h);

[x,y]=meshgrid(-2.5*pi:0.5:2.5*pi);
z=(sin(x).*sin(y))./(x.*y);;
surf(x,y,z)
colormap(hsv)
shading interp
box on
view(60,20)
figure()
[c,h]=contour(x,y,z);
clabel(c,h);

[x,y]=meshgrid(-3.0:0.1:3.0);
z=1.0-(x.^2.0+y.^2.0);
surfc(x,y,z)
colormap(copper)
shading interp
box on
view(60,20)
figure()
[c,h]=contour(x,y,z);
clabel(c,h);