clear

%[x,y] = meshgrid(-10:0.05:10);
[x,y] = meshgrid(-100:1:100);
z=((sin(sqrt(x.^2+y.^2))).^2-0.5)./((1+0.001*(x.^2+y.^2)).^2)-0.5;
mesh(x,y,z);
title('J. D. Schaffer''s Function','FontSize',16);
%axis([-3.14,3.14,-3.14,3.14,-1,1]);