计算双整数

Question

我想这样计算二重积分

如何用python编写代码？ 谢谢

def J(t, f, a, b, J, H, N, T):

   J = (M(t, f, a, b, J, H, N, T) * np.sin(t))/(4*np.pi)

   return J

Answer 1

x*y**2在x从0到2，y从0到1的盒子上的二重积分。你知道吗

from scipy import integrate
f = lambda y, x: x*y**2
print(f)
integrate.dblquad(f, 0, 2, lambda x: 0, lambda x: 1)

资料来源： https://docs.scipy.org/doc/scipy/reference/generated/scipy.integrate.dblquad.html

评论后编辑： 这两个lambda函数用于scipy docy = gfun(x)..hfun(x)中定义的y范围。如果你需要一个函数，比如说p和变量t,f，它应该给你同样的结果：

from scipy import integrate
p = lambda f, t: t*f**2
integrate.dblquad(p, 0, 2, lambda x: 0, lambda x: 1)

阅读更多关于https://mathinsight.org/double_integral_examples的二重积分

在另一条评论后编辑2：
如果您的函数包含多于f,t，那么在使用dtdf进行积分时，在开始双重积分之前，所有其他函数都必须有一个值：
示例：

from scipy import integrate
a=1
b=2
J=3
H=4
N=5
T=6
p = lambda f, t: t*f**2+(a+b+J+H+N+T)
integrate.dblquad(p, 0, 2, lambda x: 0, lambda x: 1)