Product Type Two-Point Gauss-Legendre-Simpson's
Integration [D1] (Algorithm A438)
CACM December, 1972
Boland, W. R.
