SciPy has a numerical integration package, scipy.integrate, which has no equivalent in NumPy. The quad() function can integrate a one-variable function between two points. These points can be at infinity. The function uses the simplest numerical integration method: the trapezoid rule.