f ='120*x(2)+(x(2)^2*(25-x(2))+x(1)^2*(25+x(2)))/(2*x(1))';
[x,fval,exitflag]=fmincon(f,[0;0],[],[],[],[],[12;9],[15;12])
我用的是matlab的fmincon函数求解,这里x(1)=x,x(2)=y
x =
12
9
fval =
1338
exitflag =
1
根据结果 x=12 y=9,此时的函数值为1338,exitflag = 1 表明结果收敛于X解
fun=@(x)(120*x(2)+(x(2)^2*(25-x(2))+x(1)^2*(25+x(2)))/(2*x(1)));
[x,f]=fmincon(fun,[13 10],[],[],[],[],[12 9],[15 12])
x =
12 9 %x=12,y=9
f =
1338