Return to four CVS log

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

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

1.1     ! root        1: n
        !             2: pi = 4*atan(1);
        !             3: i = sqrt(-1);
        !             4: w = exp(2*pi*i/n);
        !             5: F = <>;
        !             6: for k = 1:n, for j = 1:n, F(k,j) = w**((j-1)*(k-1));
        !             7: F = F/sqrt(n);
        !             8: alfa = r*pi;
        !             9: rho = exp(i*alfa);
        !            10: S = log(rho*F)/i - alfa*EYE;
        !            11: serr = norm(imag(S),1);
        !            12: S = real(S);
        !            13: serr = serr + norm(S-S',1)
        !            14: S = (S + S')/2;
        !            15: ferr = norm(F-exp(i*S),1)

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