让
在哪里. 我想以数字方式积分:
在哪里,, 和.
一种天真的方法在达到峰值时会出现问题远离积分矩形。
有没有我可以使用的方法?
让
在哪里. 我想以数字方式积分:
在哪里,, 和.
一种天真的方法在达到峰值时会出现问题远离积分矩形。
有没有我可以使用的方法?
scipy.integrate
我刚刚尝试了双积分( )附带的积分算法,dblquad
它似乎对您的问题很有效
from __future__ import division, print_function
import numpy as np
from scipy.integrate import dblquad
a = 2
b = 2
c = 1
fun = lambda x, y: np.exp(-a/2*x**2 - b/2*y**2 + c*x*y)*x**3
inte2 = dblquad(fun, 1, 4, lambda y: -5, lambda y: 8)
print(inte2)
并返回
(0.5814009878561697, 8.012117719083942e-09)
第一个值是积分,第二个是对积分误差的估计。