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

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.