如下做实证遇到困难，自己把参数代了几次都错了，公式见附件。

	x1	x2	x3	y1	y2
DMU1	30	25	130	35	60
DMU2	60	40	150	43	80
DMU3	55	70	120	76	53
DMU4	40	30	70	52	42
DMU5	70	90	180	63	71

求教各位大虾，这个交叉效率DEA参数怎样带进去。谢谢各位。

clear
X=[];
Y=[];
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=[];
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);
Eii=Y(:,i)'*w(m+1:m+s,i);
for k=1:n;
f=[zeros(1,m) Y(:,k)'];
Aeq=[X(:,k)'zeros(1,s) Eii*X(:,i)'-Y(:,i)']
beq=[1 0];
v=LINPROG(f,A,b,Aeq,beq,LB,UB);
E(i,k)=(Y(:,k)'*v(m+1:m+s))/(X(:,k)'*v(1:m));
end
end
E
Mean(E)
[Y I]=sort(mean(E));
fliplr(I)