AstroLee的博客

Gauss-Legendre integraldef f(x): return x**4-2*x+1N = 3 #取3个样本点a,b = 0,2 #积分上下限x,w = roots_legendre(N)xp = x*(b-a)/2+(b+a)/2wp = w*(b-a)/2s = 0for i in range(N): s+=wp[i]*f(xp[i])print(s)