clear

[x,y] = meshgrid(-65.56:0.6556:65.56);
a=[-32 -16 0 16 32 -32 -16 0 16 32 -32 -16 0 16 32 -32 -16 0 16 32 -32 -16 0 16 32 ;-32 -32 -32 -32 -32 -16 -16 -16 -16 -16 0 0 0 0 0 16 16 16 16 16 32 32 32 32 32];
l=size(x,1);
for i=1:l
    for j=1:l
        yy=0;
        for ii=1:25
            yy=yy+1/(ii+(x(i,j)-a(1,ii))^6+(y(i,j)-a(2,ii))^6);
        end
        z(i,j)=1/(1/500+yy);
    end
end
mesh(x,y,z);
title('Shekel''s Foxholes Function','FontSize',16);
axis([-70,70,-70,70,0,600])