function [sol,eval]=f553(sol,options)
m(1)=sol(1);
m(2)=sol(2);
m(3)=sol(3);
%失效概率矩陣
q=[0.01 0.05 0.10 0.18;
0.08 0.02 0.15 0.12;
0.04 0.05 0.20 0.10];
%約束條件
g1=51-(m(1)+3).^2+m(2).^2+m(3).^2;
g2=20*sum(m+exp(-m))-120;
g3=20*sum(m.*exp(-m/4))-65;
%計算加懲罰項的適值
if ((g1>=0)&(g2>=0)&(g3>=0))
  multi=1;
  for i=1:3
     summ=0;
     for j=2:4
        summ=summ+q(i,j).^(m(i)+1);
     end      
     multi=multi*(1-(1-(1-q(i,1)).^(m(i)+1))-summ);
  end
  eval=multi;
else
%取M=500
  eval=-500;
end