function f=wolunfun(x)
f=x(1)*(x(2)+x(3)*20)/2
function[c,ceq]=woluncon(x)
c(1)=15-atand(x(3)/x(2));
c(2)=atand(x(3)/x(2))-30;
c(3)=1.949*10^6*exp(-1.7762*2*x(2)/(x(2)+20*x(3)))*(1/(x(1)^3*(x(2)+20*x(3))^3))^0.5-244;
c(4)=114800/(x(1)^3*x(3)*(0.5+(x(2)+1)^0.5))-45.08;
c(5)=((118552/(x(1)*x(2)))^2+(34519/(x(1)*x(3)))^2)^0.5*(1.179*10^(-8)*x(1)^3*x(3)^3)/(x(1)*x(2)-2*x(1)-0.4)^4-0.001*x(1)*x(2);
ceq(1)=x(1)-6.3;
ceq(3)=x(3)-2;
x0=[8,10,2];
lb=[2,7,1];
ub=[10,25,3];
options=optimset('Display','iter','LargeScale','off');;
[x,fval,exitflag,output,lambda]=fmincon(@wolunfun,x0,[],[],[],[],lb,ub,@woluncon,options)