dea.rar

s=zeros(nx+ny,nx+ny); for i=1:nx s(i,i)=1; end for i=nx+1:ny+nx s(i,i)=-1; end l=zeros(m+nx+ny+1,1); f=zeros(m+nx+ny+1,l); for i=1:m+nx+ny f(i,1)=-0.000001; end f(m+nx+ny+1,1)=1; result=zeros(m,1); result1=zeros(m,1); for dmui=1:m dmux1=dmu(1:nx,dmui:dmui); dmux=[dmux1;zeros(ny,1)]; dmuy1=dmu(nx+1:ny+nx,dmui:dmui); dmuy=[zeros(nx,1);dmuy1]; a=[dmu s -dmux]; temp1=zeros(1,m); for i=1:m temp1(1,i)=1; end temp2=zeros(1,nx+ny+1); temp=[temp1,temp2]; a1=[a;temp]; b=dmuy; b1=[b;[1]]; [X,FVAL]=linprog(f,[],[],a,b,l,[]); for i=1:m+nx+ny+1 result(dmui)=FVAL; end [X,FVAL]=linprog(f,[],[],a1,b1,l,[]); result1(dmui)=FVAL; end table=[result1(1:m,1:1) result(1:m,1:1)]; save deajg table -ascii