Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- clear
- clc
- warning off
- A=[2,1;3,3;2,0];
- B=[10;24;8];
- Z=[-1;-1;-1];
- F=[300,200];
- [n,m]=size(Z);
- if n==(sum(Z)*(-1)) %+4 punkty
- AA=[A,eye(n)];
- end
- CJ=[F,zeros(1,length(Z))];
- WB=size(F,2)+1:1:size(AA,2);
- CB=CJ(WB);
- FC=CB*B;
- WW=CB*AA-CJ;
- [W,KK]=min(WW);
- while W<0
- for i=1:size(AA,1)
- if AA(i,KK)>0
- tmp(i)=B(i)/AA(i,KK);
- else
- tmp(i)=inf;
- end
- end
- [W,WK]=min(tmp);
- WB(WK)=KK;
- CB=CJ(WB);
- ER=AA(WK,KK);
- AA(WK,:)=AA(WK,:)/ER;
- B(WK)=B(WK)/ER;
- for i=1:size(AA,1)
- if i~=WK
- B(i)=B(i)-B(WK)*AA(i,KK);
- AA(i,:)=AA(i,:)-AA(WK,:)*AA(i,KK);
- end
- end
- FC=CB*B;
- WW=CB*AA-CJ;
- [W,KK]=min(WW);
- end
- FC
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement