clear

[x,y] = meshgrid(-50:0.5:50);
yy1=1+(x+1)/4;
yy2=1+(y+1)/4;
l=size(x,1);
for i=1:l
    for j=1:l
        if x(i,j)>10
            u1(i,j)=100*(x(i,j)-10)^4;
        elseif x(i,j)<-10
            u1(i,j)=100*(-x(i,j)-10)^4;
        else
            u1(i,j)=0;
        end
        if y(i,j)>10
            u2(i,j)=100*(y(i,j)-10)^4;
        elseif y(i,j)<-10
            u2(i,j)=100*(-y(i,j)-10)^4;
        else
            u2(i,j)=0;
        end
    end
end
z=0.1*((sin(3*pi*x)).^2+(x-1).^2*(1+(sin(3*pi*y)).^2)+(y-1).^2)+u1+u2;
mesh(x,y,z);
title('Generalized Penalized Function','FontSize',16);