Advertisement
Guest User

Untitled

a guest
Jan 29th, 2020
103
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.72 KB | None | 0 0
  1. clear
  2. clc
  3. warning off
  4. A=[2,1;3,3;2,0];
  5. B=[10;24;8];
  6. Z=[-1;-1;-1];
  7. F=[300,200];
  8.  
  9.  
  10. [n,m]=size(Z);
  11. if n==(sum(Z)*(-1)) %+4 punkty
  12. AA=[A,eye(n)];
  13. end
  14.  
  15. CJ=[F,zeros(1,length(Z))];
  16. WB=size(F,2)+1:1:size(AA,2);
  17. CB=CJ(WB);
  18. FC=CB*B;
  19. WW=CB*AA-CJ;
  20. [W,KK]=min(WW);
  21.  
  22. while W<0
  23. for i=1:size(AA,1)
  24. if AA(i,KK)>0
  25. tmp(i)=B(i)/AA(i,KK);
  26. else
  27. tmp(i)=inf;
  28. end
  29. end
  30. [W,WK]=min(tmp);
  31. WB(WK)=KK;
  32. CB=CJ(WB);
  33. ER=AA(WK,KK);
  34. AA(WK,:)=AA(WK,:)/ER;
  35. B(WK)=B(WK)/ER;
  36. for i=1:size(AA,1)
  37. if i~=WK
  38. B(i)=B(i)-B(WK)*AA(i,KK);
  39. AA(i,:)=AA(i,:)-AA(WK,:)*AA(i,KK);
  40. end
  41. end
  42. FC=CB*B;
  43. WW=CB*AA-CJ;
  44. [W,KK]=min(WW);
  45. end
  46. FC
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement