不太明白楼主意思，对数字求导不成0了么？是不是说这一系列的数是将一组自变量值代入函数后得到的函数值？那要用数值方法求导数了

对，就是你说的那个意思，明白了，我用数值求导

SUBROUTINE Diff(N, DX, P, DP)
! COMPUTES THE FIRST DERIVATIVE DP OF P USING
! RICHARDSON EXTRAPOLATION FORMULA. THE DERIVATIVE AT THE
! BOUNDARIES ARE ASSUMED TO BE ZERO

     IMPLICIT NONE
      INTEGER :: N
      REAL*8 :: DX
      REAL*8 :: P(N), DP(N)
      INTEGER :: I

      DP(1)=0.0D0
      DP(2) = (P(3) - P(1))/(2.0D0*DX)
      DP(3) = (P(4) - P(2))/(2.0D0*DX)
      DO I = 4, N-3
          DP(I) = (P(I-2)-8.0D0*P(I-1)+8.0D0*P(I+1)-P(I+2))/(12.0D0*DX)
      END DO
      DP(N-2) = (P(N-1) - P(N-3))/(2.0D0*DX)
      DP(N-1) = (P(N) - P(N-2))/(2.0D0*DX)
      DP(N) = 0.0D0
      RETURN
      END SUBROUTINE DIFF