Return to jacstep CVS log

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

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

1.1       root        1:         d = (a(q,q)-a(p,p))*0.5/a(p,q);
                      2:         t = 1/(abs(d)+sqrt(d*d+1));
                      3:         if d < 0, t = -t; end;
                      4:         c = 1/sqrt(1+t*t);  s = t*c;
                      5:         R = eye(n); r(p,p)=c; r(q,q)=c; r(p,q)=s; r(q,p)=-s;
                      6:         X = X*R;
                      7:         A = R'*A*R;

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