Pātuhitia tēnei: fortran subroutine for simpson's method of numerical integration /