subroutine ssinti (n,wsave) real wsave(1), dt, fk, pi data pi / 3.141592653 5897932384 6264338327 950 e0 / c if (n .le. 1) return np1 = n+1 ns2 = n/2 dt = pi/float(np1) fk = 0.e0 do 101 k=1,ns2 fk = fk+1.e0 wsave(k) = 2.e0*sin(fk*dt) 101 continue c call srffti (np1,wsave(ns2+1)) c return end