Gaussian quadratures Calculation of the nodes and weights of the Gaussian quadratures
Simpson's method One-dimensional integration
Adaptive integration using a 61-point Gauss-Kronrod quadrature One-dimensional integration