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