Return to qr CVS log

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

Annotation of researchv10dc/cmd/matlab/exec/qr, revision 1.1

1.1     ! root        1: scale = s(m);
        !             2: sm = s(m)/scale; smm1 = s(m-1)/scale; emm1 = e(m-1)/scale;
        !             3: sl = s(l)/scale; el = e(l)/scale;
        !             4: b = ((smm1 + sm)*(smm1 - sm) + emm1**2)/2;
        !             5: c = (sm*emm1)**2;
        !             6: shift = sqrt(b**2+c); if b < 0, shift = -shift;
        !             7: shift = c/(b + shift)
        !             8: f = (sl + sm)*(sl-sm) - shift
        !             9: g = sl*el
        !            10: for k = l: m-1, exec('qrstep',ip)
        !            11: e(m-1) = f

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