function [sol,eval]=f552(sol,options)
x1=sol(1);
x2=sol(2);
r1=0.1;
r2=0.8;
%約束條件
g1=x1-2*x2+1;
g2=x1.^2/4-x2.^2+1;
%加懲罰項的適值
if (g1>=0)&(g2>=0)
  eval=(x1-2).^2+(x2-1).^2;
else
  eval=(x1-2).^2+(x2-1).^2+r1*g1+r2*g2;
  eval=-eval;
end