在SciPy中,可以使用scipy.integrate.quad
函數來實現自適應積分算法。該函數可以通過指定參數epsabs
和epsrel
來控制積分的絕對誤差和相對誤差的精度,從而實現自適應積分。具體的用法如下:
import scipy.integrate as spi
def integrand(x):
return x**2
result, error = spi.quad(integrand, 0, 1, epsabs=1.0e-6, epsrel=1.0e-6)
print(result)
在上面的例子中,quad
函數會自適應地對integrand
函數在區間[0, 1]上進行積分,并返回積分結果和估計的誤差。用戶可以根據需要調整epsabs
和epsrel
的值來控制積分的精度。