|
|
1.1 ! root 1: <n,n> = size(A); ! 2: q1 = rand(n,1); ! 3: ort ! 4: alfa = <>; beta = <>; ! 5: q = q1/norm(q1); r = A*q(:,1); ! 6: for j = 1:n, exec('lanstep',0);
This archive runs on limited infrastructure. Preserving old code on modern bandwidth. Automated agents are requested to crawl responsibly.