【正文】
5]。 N=[3。 3 10 0 0 1 300。運行結(jié)果如下: A=[9 4 1 0 0 360。 for i=1:mA if i~=outb A(i,:)=A(i,:)A(outb,:)*A(i,inb)。 % 出基變量下標(biāo) end end % 以下更新N for i=1:mA1 if i==outb N(i)=inb。sita(i)temp temp=sita(i)。 end end temp=inf。 % 進基變量的下標(biāo) end end sita=zeros(1,mA1)。 for i=1:nA1 if A(mA,i)temp temp=A(mA,i)。 break。)。A(1:mA1,i)=0 % 問題有無界解 disp(39。 end val=A(mA,nA)。 sol=zeros(1,nA1)。while flag kk=kk+1。 kk=0。 . A*x=b。4]。 0]。運行結(jié)果如下: A=[1 1 1 0 4。 for i=1:mA if i~=outb A(i,:)=A(i,:)A(outb,:)*A(i,inb)。 % 出基變量下標(biāo) end end % 以下更新N for i=1:mA1 if i==outb N(i)=inb。sita(i)temp temp=sita(i)。 end end temp=inf。 % 進基變量的下標(biāo) end end sita=zeros(1,mA1)。 for i=1:nA1 if A(mA,i)temp temp=A(mA,i)。 break。)。A(1:mA1,i)=0 % 問題有無界解 disp(39。 end val=A(mA,nA)。 sol=zeros(1,nA1)。while flag kk=kk+1。 kk=0。 . A*x=b。0],0,[2。運行結(jié)果如下:[x,fx,k]=bowell([2 2。*A*x+b39。endx=x0。 d02=d10。 x10=x0+a01*d01。*g10/(d1039。 g10=A*x02+b。 x02=x01+a02*d02。*g02/(d0239。 g02=A*x01+b。*A*d01)。 a01=d0139。%