|
DataMuseum.dkPresents historical artifacts from the history of: RC4000/8000/9000 |
This is an automatic "excavation" of a thematic subset of
See our Wiki for more about RC4000/8000/9000 Excavated with: AutoArchaeologist - Free & Open Source Software. |
top - metrics - download
Length: 768 (0x300) Types: TextFile Names: »extinvert«
└─⟦621cfb9a2⟧ Bits:30002817 RC8000 Dump tape fra HCØ. Detaljer om "HC8000" projekt. └─⟦0364f57e3⟧ └─⟦58ca399f1⟧ »extbib« └─⟦this⟧
invert=algol list.yes index.no external procedure invert(n,a,b); value n; integer n; array a,b; begin integer i,j; integer array p(1:n); array bj(1:n); comment the procedure inverts the nonsingular matrix a(1:n,1:n) and stores the result in b; decompose(a,p,0); for j:=1 step 1 until n do begin for i:=1 step 1 until n do bj(i):= 0; bj(j):= 1; solve(a,p,0,bj); for i:=1 step 1 until n do b(i,j):= bj(i) end end; end ▶EOF◀