Simpson's Rule

 

In either Fortran or C, write a function subroutine of the form

C:        float simpint(float f[],float r[],float dr[],int N)

Fortran:  real function simpint(f,r,dr,N)
          integer N
          real f(N),r(N),dr(N)

which returns the result of Simpson's rule for the same integral

using the same inputs as your routine in (1.1).