quad函數是SciPy庫中的一個函數,用于數值積分。它可以計算指定區間上函數的定積分。quad函數的語法如下: quad(fun, a, b[, limit=None[, args=()]], epsabs=1.49e-8, epsrel=1.0e-12) 其中:
import numpy as np
from scipy.integrate import quad
# 定義要積分的函數
def f(x):
return x**2 + 1
# 計算區間[0, 1]上的定積分
result, error = quad(f, 0, 1)
print("積分結果:", result)
print("誤差估計:", error)
輸出結果為:
積分結果: 1.3333333333333335
誤差估計: 1.0942262360578078e-08