%This script runs myfun
%fancy plots
%MA 302

[X,Y]=meshgrid(linspace(-5,5,50));

c=2.5;
figure
Z = myfun(X,Y,c);
colormap hsv
surf(X,Y,Z,'FaceColor','interp',...
   'EdgeColor','none',...
   'FaceLighting','gouraud')
daspect([5 5 2.5])
axis tight
view(-50,30)
camlight left


c=5;
figure
Z = myfun(X,Y,c);
colormap hsv
surf(X,Y,Z,'FaceColor','interp',...
   'EdgeColor','none',...
   'FaceLighting','gouraud')
daspect([5 5 2.5])
axis tight
view(-50,30)
camlight left
