x=linspace(-%pi,%pi,5) [x,y]=field(x,x); z= 3*sin(x).*cos(y) [xb,yb,zb]=beziersurface(x,y,z); xsetech([0,0,1.0,0.5]); plot3d3(x,y,z); xtitle('A first surface '); xsetech([0,0.5,1.0,0.5]) plot3d2(xb,yb,zb,-1,35,45," ",[4,2,3]); xtitle('The bezier interpolated surface (n=10)'); xsetech([0,0,1,1]);