t=%pi*(-10:10)/10; deff('[z]=surf(x,y)','z=sin(x)*cos(y)'); rect=[-%pi,%pi,-%pi,%pi,-5,1]; z=feval(t,t,surf); contour(t,t,z,10,35,45,'X@Y@Z',[1,1,0],rect,-5); plot3d(t,t,z,35,45,'X@Y@Z',[2,1,3],rect); xtitle('plot3d and contour ',' ',' ');