anchor.diff:  
         IMPLICIT none
         INTEGER NMAX
  -      PARAMETER (NMAX=50000)
  +      PARAMETER (NMAX=100000)
         REAL*8 x(NMAX),f(NMAX),fint(NMAX),dx
         CHARACTER*30 ifile,ofile
         LOGICAL ex
  
positive.code:      1 '  -h/-?    Online help (this screen)'
        write(*,*)
        write(*,*)
       1 'polyfit numerically integrates a time series using trapezoidal',
       2 ' integration'
        write(*,*)
        write(*,*)
  
negative.code:       fint(1) = 0.0d0
        do 103 i = 2,ntot
            dx = x(i) - x(i-1)
            fint(i) = fint(i-1) + dx*0.5d0*(f(i)+f(i-1))
    103 continue
  
  