|
|
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◀