function S = penal(x)
global r
S = x(1)*x(2)+r*(x(1).^2+x(2).^2+x(3).^2-25).^2;
endfunction
