一、SciPy求函数积分
1、已知函数型求积分
一般情况
from scipy import integrate
import numpy as np
def f(x):
return x + 1
v , err =integrate.quad(f, -1, 1)
print(format(v, ".2f")) # 保留两位小数
遇到未知的参数的函数
from scipy import integrate
import numpy as np
def f(x, a, b):
return a * x + b
v , err =integrate.quad(f, -1, 1, args=(-1,1))
print(format(v, ".2f")) # 保留两位小数
遇到有断点的函数
from scipy import integrate
import numpy as np
def f(x):
return 1 / np.sqrt(abs(x))
v , err =integrate