function y=fm(x)
    %y=8*x(1)^2+4*x(1)*x(2)+5*x(2)^2;
    %y=4*x(1)^2+3*x(2)^2-4*x(1)*x(2)+x(1);    
    %y=(x(1)-4)^2+(x(2)-5)^2;
    %y=(x(1)-100)^2+(x(2)-150)^2+(x(3)-620)^2;
    %y=100*(x(2)-x(1)^2)^2+(1-x(1)^2)^2;
    
    y=100*(x(2)-x(1)^2)^2+(1-x(1))^2+90*(x(4)-x(3)^2)^2+(1-x(3))^2+10.1*((x(2)-1)^2+(x(4)-1)^2)+19.8*((x(2)-1)*(x(4)-1));