import scipy as sp
import numpy as np
from sympy import *
import matplotlib.pyplot as plt
def fun(x):
return [np.sin(el)/el for el in x]
end = 100
trap = 10000
x = np.linspace(0.001, end, trap)
y = fun(x)
#分割近似求積分
print sum(y)*end/trap*2
print scipy.integrate.quad(lambda x:np.sin(x)/x, 0.001,np.inf)
print scipy.integrate.quad(lambda x:np.exp(-x), 0,np.inf)
z = symbols('z')
print sympy.integrate(sin(z)/z, (z, 0, oo))
plt.plot(x,y)
plt.show()
scipy sin(x)/x算出來(lái)積分誤差很大，即使近似值已經(jīng)很接近pi了,注意我乘2了
