Scipy中积分的计算

 一重积分的运算

形如f(x)={k(x**n)dx,首先积分运算需先引用scipy.intergrate,一重积分引入quad(),

用法：quad(func,a,b,args(n,k))，func是函数体，a是下限，b是上限，n是多少次方，k代表常系数。

# 计算3*x**0到3*x**14的积分结果
from scipy.integrate import quad
import matplotlib.pyplot as plt

def func(x, n, k):
    return k * x ** n

result = []
for i in range(15):
    data = quad(func, 0, 2, args=(i, 3))
    result.append(data[0])
# data[0] 代表只取第一个数值
print(result)

# 图像表示
plt.plot(result)
plt.show()

二重积分的运算

形如f(x)={{k(x**2)*(y**2)dxdy,二重积分引入dblquad()

用法：dblquad(func,a,b,gfun=c(常数），hfun=f(x)),func是函数体，c1表示下限，c2表示上限(c1,c2为写在积分式子前面的上下限），gfun为下限，hfun为上限，一般hfun为一个函数式子，需在前面提前定义。

# 计算{{3*(x**2)*(y**2)dxdy 的积分
from scipy.integrate import dblquad

def fun(x, y):
    return 3 * (x ** 2) * (y ** 2)

def y_area(x):
    return 1 - x ** 2

data = dblquad(fun, c1, c2, gfun=c3, hfun=y_area)
print(data[0])