大家帮帮忙，很急的，谢谢拉，我怎么也算不出来，如果谁知道请把程序和结果发到我邮箱或者这里。邮箱taospring@163.com 谢谢大家了

f='(0.1576-(0.1501*x(1)+0.8499*x(2)))^2+(0.1184-(0.1576*x(1)+0.8424*x(2)))^2+(0.1443-(0.1184*x(1)+0.8816*x(2)))^2+(0.8424-(0.1501*x(3)+0.8499*x(4)))^2+(0.8816-(0.1576*x(3)+0.8424*x(4)))^2+(0.8557-(0.1184*x(3)+0.8816*x(4)))^2';
其中x（1）＋x（3）＝1；x（2）＋x（4）＝1
x（i）大于等于零；i＝1，2，3，4。

谢谢大家了，教教我啊，我晚上在线等。谢谢

function f=mmfun(x)
f=(0.1576-(0.1501*x(1)+0.8499*x(2)))^2+(0.1184-(0.1576*x(1)+0.8424*x(2)))^2+(0.1443-(0.1184*x(1)+0.8816*x(2)))^2+(0.8424-(0.1501*x(3)+0.8499*x(4)))^2+(0.8816-(0.1576*x(3)+0.8424*x(4)))^2+(0.8557-(0.1184*x(3)+0.8816*x(4)))^2;

x0=[0;0;0;0];
A=[1 0 1 0;
-1 0 -1 0;
0 1 0 1;
0 -1 0 -1];
b=[1 -1 1 -1];
lb=[0 0 0 0];
[x,fval]=fmincon(@mmfun,x0,A,b,[],[],lb,[])