Anfonwch hwn fel neges destun: fortran subroutine for simpson's method of numerical integration /