|
Identification Division.
Program-Id. pivot.
author. Robert Sedgewick.
* Simplex Verfahren, lineare Programmierung
* aus Algorithmen p698
Data Division.
Working-Storage Section.
77 j pic 9(4).
77 k pic 9(4).
Linkage Section.
01 a.
05 pic 9(4) occurs 1000.
05 pic 9(4) occurs 1000.
77 p pic 9(4).
77 q pic 9(4).
77 best pic 9(4)V9(2) occurs 1000.
Procedure Division using cost val best.
perform varying j from 0 until j=n
perform varying k from m+1 by -1 until k=1
if j not = p and k not = q then
compute a(j,k)=a(j,k)-a(p,k)*a(j,q)/a(p,q)
end-perform
end-perform
perform varying j from 0 until j=n
if j not = p then
move 0 to a(j,q)
end-perform
perform varying k from 1 until k=m+1
if k not = q then
compute a(p,k)=a(p,k)/a(p,q)
end-perform
move 1 to a(p,q)
End-Program pivot.
Identification Division.
Program-Id. simplex.
author. Robert Sedgewick.
* Simplex Verfahren, lineare Programmierung
* aus Algorithmen p698
Data Division.
Working-Storage Section.
77 j pic 9(4).
77 k pic 9(4).
Linkage Section.
01 a.
05 pic 9(4) occurs 1000.
05 pic 9(4) occurs 1000.
77 p pic 9(4).
77 q pic 9(4).
77 best pic 9(4)V9(2) occurs 1000.
Procedure Division using cost val best.
perform with test after
until (q=m+1) or (p=n+1)
move 0 to q
perform with test after
until (q=m+1) or (a(0,q)<0)
add 1 to q
end-perform
move 0 to p
perform with test after
until (p=n+1) or (a(p,q)>0)
add 1 to p
end-perform
perform varying i from p+1 until i=n
if a(i,q)>0 then
if a(i,m+1)*a(i,q) < a(p,m+1)/a(p,q) then
move i to p
end-perform
if (q<m+1) and (p<n+1) then
call pivot using a p q
end-perform
End-Program pivot.
[ zur Elbe Produktseite wechseln0.18Quellennavigators
Analyse erneut starten
]
|
