%
% Problem Epage 202
% 
% Written by:
% -- 
% John L. Weatherwax                2008-02-20
% 
% email: wax@alum.mit.edu
% 
% Please send comments and especially bug reports to the
% above email address.
% 
%-----

clear all; 
close all; 
clc; 

ezsurf( 'sin(sqrt(x^2+y^2))/sqrt(x^2+y^2)',[-6*pi,+6*pi]); 
% add some lighting
view(0,75); 
shading interp; 
lightangle(-45,30);
set(findobj('type','surface'), ...
    'FaceLighting', 'phong', 'AmbientStrength', 0.3, 'DiffuseStrength', 0.8, ...
    'SpecularStrength', 0.9, 'SpecularExponent', 25, ...
    'BackFaceLighting', 'unlit' ); 
axis off; 
saveas( gcf, 'shading_options_example', 'epsc' );