x=[6行2列的矩阵];
y=[6行25列矩阵];
n=size(x',1);
m=size(x,1);
s=size(y,1);
a=[-x' y'];b=zeros(n,1);
lb=zeros(m+s,1);ub=[];
options=optimset('maxiter',1e4,'tolfun',1e-6);
for i=1:n
    aeq=[x(:,i)' zeros(1,s)];beq=1;   
     f=[zeros(1,m) -y(:,i)'];
     w(:,i)=linprog(f,a,b,aeq,beq,lb,ub,[],options);
     E(i,i)=y(:,i)'*w(m+1:m+s,i);
     for k=1:n
          f=[zeros(1,m) y(:,k)'];
          aeq=[x(:,k)' zeros(1,s);E(i,i)*x(:,i)' -y(:,i)'];
          beq=[1;0];
      v=linprog(f,a,b,aeq,beq,lb,ub,[],options);
       E(i,k)=(y(:,k)'*v(m+1:m+s))/(x(:,k)'*v(1:m));
     end
end
该程序出现错误，请给予改正，您可设置出售帖，1-500币悬赏，请讲诚信