1      SUBROUTINE SPECTOP(Dr,N)
2      DIMENSION d1(0:32,0:32) , Dr(0:32,0:32) , x(0:32)
3      DO k = 0 , N
4         fctr2 = o
5         DO j = 0 , N
6            fctr = fctr1*fctr2
7            IF ( j.NE.k ) THEN
8               d1(k,j) = ck*fctr/(cj*(x(k)-x(j)))
9            ENDIF
10            fctr2 = -o*fctr2
11         ENDDO
12         DO j = 0 , N
13            Dr(k,j) = d1(N-k,N-j)
14         ENDDO
15      ENDDO
16      END
17