Download this file
1 2 3 4 5 6 7 8 9 10
>>> result = quad(lambda x: expint(3, x), 0, Inf) >>> print result (0.33333333324560266, 2.8548934485373678e-09) >>> I3 = 1.0/3.0 >>> print I3 0.333333333333 >>> print I3 - result[0] 8.77306560731e-11