Seol mar théacs é seo: fortran subroutine for simpson's method of numerical integration /