Return to jacobi CVS log

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

researchv10 Dan Cross

<n, n> = size(A); 
X = eye(n);
anorm = norm(A,'fro');
cnt = 1;
while cnt > 0,...
  cnt = 0;...
  for p = 1:n-1,...
    for q = p+1:n,...
      if anorm + abs(a(p,q)) > anorm,...
        cnt = cnt + 1;...
        exec('jacstep');...
      end,...
    end,...
  end,...
  display(rat(A)),...
end

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