这可能更像是一个数学问题，但我不知道如何在python中编程 我需要整合这种形式的功能：
(r1 + r2)**j * (r1 - r2)**k * (numpy.abs(r_1 - r_2))**m * numpy.exp(-Z*(r1+r2) / (k * r0)
关于r1和r2在0和无穷大之间。 j，k和m将会变化，但是在整合之前定义的整数，所以它们不应该太重要。 k和r0也是在进行积分之前定义的常量。
我正在尝试使用scipy.integrate.dblquad。然而问题是r_1和r_2是矢量，即（x，y，z），然后取r1-r2的模数。有没有办法告诉python这些是向量，它需要将它们整合为向量？
希望我已经解释了我的问题。谢谢