x=[45681.64;38.9;845.17];
y=[33469.23;11169.48;337.86];
n=size(x',1);
m=size(x,1);
s=size(y,1);
epsilon=1.0e-10;
f=[zeros(1,n) -epsilon*ones(1,m+s) 1];
A=zeros(1,n+m+s+1);
b=0;
LB=zeros(n+m+1,1);
UB=[];
LB(n+m+s+1)=-inf;
for i=1:n
Aeq=[x eye(m) zeros(m,s) -x(:,i); y zeros(s,m) -eye(s) zeros(s,1)];
beq=[zeros(m,1) y(:,i)];
W(:,i)=linprog(f,A,b,Aeq,beq,LB,UB);
end
W
u=W(1:n,:)
s1=W(n+1:n+m,:)
s2=W(n+m+1:n+m+s,:)
p=W(n+m+s+1,:)
Optimization terminated.

W =

    1.0000
    0.0000
    0.0000
    0.0000
    0.0000
    0.0000
    0.0000
    1.0000

u =

    1.0000

s1 =

  1.0e-007 *

    0.4322
    0.0004
    0.0080

s2 =

  1.0e-006 *

    0.5397
    0.1801
    0.0054

p =

    1.000

请问哪个是结果？