[shengsheng]

改写如下：
function Untitled
x=[4,14,0.3];
ib=[2.5,5,0.25];
ub=[10,inf,0.33];
options=optimset('largescale','off');
[x,fval]=fmincon(@gear,x,[],[],[],[],ib,ub,@mycon)
function f=gear(x)
f=0.6672*x(1)^3*x(2)^3*x(3)*[1+(1-x(3))^2+(1-x(3))]-0.4858*x(2)^(-2)-0.188;
%约束条件的M文件：
function [g,q]=mycon(x)
g(1)=160057/[x(1)^3*x(2)^3*(1-0.5*x(3))^2*x(3)]^(-0.5)-664;
g(2)=398121*x(1)^(-3)*x(2)^(-2)*(1-0.5*x(3))^(-2)*x(3)^(-1)-480;
g(3)=0.12+4.858*x(2)^(-1);
g(4)=2.5-x(1);
g(4)=0.25-x(3);
g(6)=x(3)-0.33;
q=[];