import numpy as np
import scipy.integrate as integrate

result = integrate.quad(lambda x: np.sin(x),
                        0, np.pi)
print "Result =", result[0]
print "Estimated error =", result[1]

# np.inf for infinity
result = integrate.quad(lambda x: np.exp(-x**2),
                        -np.inf, np.inf)
print "Result =", result[0]
print "Estimated error =", result[1]
print "sqrt(pi) =", np.sqrt(np.pi)