blob: 4bebed016782973257988944e6b4dd2fbf400791 [file] [log] [blame]
SUBROUTINE SPECTOP(Dr,N)
DIMENSION d1(0:32,0:32) , Dr(0:32,0:32) , x(0:32)
DO k = 0 , N
fctr2 = o
DO j = 0 , N
fctr = fctr1*fctr2
IF ( j.NE.k ) THEN
d1(k,j) = ck*fctr/(cj*(x(k)-x(j)))
ENDIF
fctr2 = -o*fctr2
ENDDO
DO j = 0 , N
Dr(k,j) = d1(N-k,N-j)
ENDDO
ENDDO
END