f = @(x,y) x .* sin(x) .* sin(y) ;
n = 40 ;
vx = linspace(0, 2*pi, n) ;
[x, y] = meshgrid(vx, vx) ;
z = f(x,y) ;
fs = surf(x, y, z) ;
xlabel('x') ;
ylabel('y') ;
zlabel('z') ;
ax = gca ;
exportgraphics(ax,'surf.pdf','ContentType','vector') ;
return
fs.EdgeColor = 'none' ;
fs.FaceColor = [0 .8 0] ;
fs.FaceLighting = 'gouraud' ;
light ;