SciPy - integrate.romberg() 方法

SciPy integrate.romberg() 方法用於計算數值積分。該方法返回區間 (a,b) 上的函式積分。

語法

以下為 SciPy integrate.romberg() 方法的語法 −

romberg(func, a, b)

引數

此函式接受以下引數 −

• func：此引數用於積分效能。
• a：此引數表示積分的下限。
• b：此引數表示積分的上限。

返回值

此方法以浮點值作為結果。

示例 1

以下是對 SciPy integrate.romberg() 方法的基本演示，該方法顯示了區間 0 和 1。

在此，我們採用函式 f(x) = x² 來計算積分。

from scipy import integrate

def polynomial(x):
    return x**2

res = integrate.romberg(polynomial, 0, 1)
print("The result of integrating x^2 from 0 to 1:", res)

輸出

以上程式碼生成以下輸出 −

The result of integrating x^2 from 0 to 1: 0.3333333333333333

示例 2

在此，我們使用 cos() 函式執行區間 0 和 pi/2 上的積分任務。

import numpy as np
import scipy.integrate as sp
def cosine(x):
    return np.cos(x)

res = sp.romberg(cosine, 0, np.pi/2)
print("The result of integrating cos(x) from 0 to π/2:", res)

輸出

以上程式碼生成以下輸出 −

The result of integrating cos(x) from 0 to π/2: 0.63212055882857