subroutine scosqi (n,wsave) real wsave(1), dt, fk, pih data pih / 1.570796326 7948966192 3132169163 975 e0 / c dt = pih/float(n) fk = 0.e0 do 101 k=1,n fk = fk+1.e0 wsave(k) = cos(fk*dt) 101 continue c call srffti (n,wsave(n+1)) c return end