>>> from scipy import integrate
>>> x = integrate.quad(lambda x: x ** 2, 0, 1)
>>> x
(0.33333333333333337, 3.700743415417189e-15)
>>> x = integrate.quad(lambda x: x ** 2, 0, 3)
>>> x
(9.000000000000002, 9.992007221626411e-14)
>>> x = integrate.quad(lambda x: x ** 2, 0, 2)
>>> x
(2.666666666666667, 2.960594732333751e-14)
>>> 8 / 3 # (1 / 3) * x ** 3
2.6666666666666665
>>> from math import sin
>>> x = integrate.quad(sin, 0, 3.1415)
>>> x
(1.9999999957076557, 2.2204460444848536e-14)