高手帮忙看下这个matlab程序哪有问题，谢谢

x=[60,70,85,106,35;156,200,157,263,105;50,180,100,86,30];    
y=[80,60,90,96,30;12,13,20,17,8;27,25,15,28,3;4,2,5,5,1];
n=size(x', 1);
m=size(x, 1);
s=size(y, 1);
epsilon=10^-10;
f=[zeros(1, n) -epsilon*ones(1, m+s) 1];
A=zeros(1, n+m+s+1);
b=0;
LB=zeros(n+m+s+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
lambda=w(1:n,:);
disp(lambda);
s_minus=w(n+1:n+m,:);
disp(s_minus);
s_plus=w(n+m+1:n+m+s,:);
disp(s_plus);
theta=w(n+m+s+1,:);
disp(theta)