pysimp.f

Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file pysimp.f

 
C*********************************************************************
 
C...PYSIMP
C...Simpson formula for an integral.
 
      FUNCTION pysimp(Y,X0,X1,N)
 
C...Double precision and integer declarations.
      IMPLICIT DOUBLE PRECISION(a-h, o-z)
      IMPLICIT INTEGER(i-n)
      INTEGER pyk,pychge,pycomp
 
C...Local variables.
      DOUBLE PRECISION y,x0,x1,h,s
      dimension y(0:n)
 
      s=0d0
      h=(x1-x0)/n
      DO 100 i=0,n-2,2
        s=s+y(i)+4d0*y(i+1)+y(i+2)
  100 CONTINUE
      pysimp=s*h/3d0
 
      RETURN
      END