DataMuseum.dk

top - metrics - download

    Length: 1536 (0x600)
    Types: TextFile
    Names: »algflinda3«

└─⟦621cfb9a2⟧ Bits:30002817 RC8000 Dump tape fra HCØ.  Detaljer om "HC8000" projekt.
    └─⟦0364f57e3⟧ 
        └─⟦7e928b248⟧ »algbib« 
            └─⟦this⟧ 

TextFile

;kemlab5 1
clear flinda
gonewtext=set 45
permanent gonewtext.17
gonewtext=edit algflinda
l./array T(1:NK/,r/TD/X1,TD/,l-1,r/steps/steps,solvsteps/
l12,r//solvsteps:= 0; /,l21,r/ q:= 0;//,l1,
r/ begin/ p:= p + a(i,j)**2;/,l1,d2,r/ X(j):= q;//
l2,i/
for j:=1 step 1 until NK do X1(j):= 0;
H:= 0;
repsolve:
/,l2,r/ q/
   q/,r/;/;  p:= 0;/,l1,d2,i/
   for i:=1 step 1 until Na do p:= p + a(i,j)*a(i,NC);
   for k:=j+1 step 1 until NK do  p:= p - T(j,k)*TD(k);
   TD(j):= q*p; X1(j):= X1(j) + TD(j); h:= h + X1(j)**2
/,l2,i?
setposition(z,0,0); h:= h*DF;
for i:=1 step 1 until Na do begin
   inrec(z,NC); q:= 0;
   for j:=1 step 1 until NK do q:= q + z(j)*X1(j);
   q:= a(i,NC):= z(NC)-q;
   h:= h + q*q
end;
if H=0 or h<H then begin
   for j:=1 step 1 until NK do X(j):= X1(j);
   solvsteps:= solvsteps+1;
   H:= h; goto repsolve
end;
?,l./DFsteps/,r/dfsteps/dfsteps,solvsteps/,f
i gonewtext
▶EOF◀