Return to jacobi CVS log

Up to [Research Unix] / researchv10dc / cmd / matlab / exec

Annotation of researchv10dc/cmd/matlab/exec/jacobi, revision 1.1.1.1

1.1       root        1: <n, n> = size(A); 
                      2: X = eye(n);
                      3: anorm = norm(A,'fro');
                      4: cnt = 1;
                      5: while cnt > 0,...
                      6:   cnt = 0;...
                      7:   for p = 1:n-1,...
                      8:     for q = p+1:n,...
                      9:       if anorm + abs(a(p,q)) > anorm,...
                     10:         cnt = cnt + 1;...
                     11:         exec('jacstep');...
                     12:       end,...
                     13:     end,...
                     14:   end,...
                     15:   display(rat(A)),...
                     16: end

This archive runs on limited infrastructure. Preserving old code on modern bandwidth. Automated agents are requested to crawl responsibly.