# import ... # result = ... # START CUT import numpy as np import scipy.integrate as integrate result = integrate.quad(lambda x: np.exp(-x), 0, np.inf) # END CUT print "Result =", result[0] print "Estimated error =", result[1]