%編寫目標函數
function y=target67(x)
y=-x(3)/(x(3)+x(4));

%編寫非線性約束函數
function [c,ceq]=fun67(x)
%無非線性不等式約束
c=[];
   
r0=3;
r1=10*x(1)/(0.3333+x(1));
r2=4.1667*x(2)/(0.6667*(1+x(3)/(1+0.5*x(1)))+x(2));
r3=6.5*x(2)/(0.6429*(1+x(4)/(7*(1+0.6667*x(3)))))+...
    x(2)*(1+x(4)/(3.5+2.3333*x(3)));
r4=6.25*x(3)/(0.1875+x(3));
r5=3.516*x(4)/(0.3333+x(4));
%非線性等式約束
ceq(1)=r0-r1*x(5);
ceq(2)=r1*x(5)-r2*x(6)-r3*x(7);
ceq(3)=r2*x(6)-r4*x(9);
ceq(4)=r3*x(7)-r5*x(9);

%優化變量邊界約束
lb=[0.1285 0.1603 0.2699 0.4668 0.2 0.2 0.2 0.2 0.2]';
ub=[0.1571 0.1959 0.3199 0.5706 5.0 5.0 5.0 5.0 5.0]';
%無限型約束
A=[];
b=[];
Aeq=[];
beq=[];
%設定初始值
x0=[0.1428 0.1781 0.2999 0.5187 0.5 0.5 0.5 0.5 0.5]';
%調用優化函數
[x,fval]=fmincon('target67',x0,A,b,Aeq,beq,lb,ub,'fun67')