Bezier surface test

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]);