SMS dit: fortran subroutine for simpson's method of numerical integration /