SciPy - integration.quad() 方法

SciPy integrate.quad() 方法用于执行定积分任务。它通常被称为四边形(两点的集合)。

语法

以下是 SciPy integrate.quad() 方法的语法 −

quad(custom_func, upper_lim, lower_lim)

参数

此方法接受以下参数 −

  • custom_func:此参数设置用户定义的积分计算函数。
  • upper_lim:此参数用于设置上限。
  • lower_lim:此参数用于设置下限。

返回值

此方法返回数值整数或绝对误差作为结果。

示例 1

以下是说明 SciPy integrate.quad() 方法用法的基本示例。

from scipy import integration

# 定义函数
def fun(x):
    return x**2

# 执行集成
res, err = integration.quad(fun, 0, 1)

print("结果为 ", res)
print("估计误差为 ", err)

输出

上述代码产生以下输出 −

结果为 0.33333333333333337
估计误差为3.700743415417189e-15

示例 2

此程序定义函数 f(x) = sin(x) 在 0pi 区间内的定积分。

import scipy.integrate as sp
import numpy as np

# 定义函数
def fun(x):
    return np.sin(x)

# 执行积分
res, err = sp.quad(fun, 0, np.pi)

print("结果为 ", res)
print("估计误差为 ", err)

输出

上述代码产生以下输出 −

结果为 2.0
估计误差为 2.220446049250313e-14

示例 3

在这里,我们在用户定义函数中说明指数函数(e-x2),并使用 quad() 执行结果积分。这里,范围间隔设置在 -infinf 之间。

import scipy.integrate as sp
import numpy as np

# 定义函数
def fun(x):
    return np.exp(-x**2)

# 执行积分
res, err = sp.quad(fun, -np.inf, np.inf)

print("结果为 ", res)
print("估计误差为 ", err)

输出

上述代码产生以下输出 −

结果为 1.7724538509055159
估计误差为1.4202636780944923e-08

scipy_reference.html