Python 和 scipy 中的 Meijer G 函数
33
0
我需要一个 scipy 中的 Meijer G 函数。我在网上读到过，由于 Meier G 函数过于通用，所以 scipy 并不支持将其作为特殊函数，但大家应该...
我需要一个 scipy 中的 Meijer G 函数。我在网上读到过，由于 Meier G 函数过于通用，所以 scipy 不支持将其作为特殊函数，但每个人都应该根据自己的使用情况编写一些函数。
我的问题是，我对复杂积分毫无经验。由于这里禁止使用 LaTeX，因此我尝试用数字方法解决以下问题：
（第一行是一般情况，第二行是我要计算的情况），其中 p(a)、k、k2 给定
正如 维基百科所述 ，有三种方法可以获取 L :
- L 从 −i∞ 延伸到 +i∞，使得 Γ(bj − s)，j = 1, 2, ..., m 的所有极点均位于路径的右侧，而 Γ(1 − ak + s)，k = 1, 2, ..., n 的所有极点均位于路径的左侧。
- L 是一个始于并止于 +∞ 的环，围绕 Γ(bj − s), j = 1, 2, ..., m 的所有极点，在负方向上恰好一次，但不围绕 Γ(1 − ak + s), k = 1, 2, ..., n 的任何极点。
- L 是一个始于并止于 −∞ 的环，它绕 Γ(1 − ak + s), k = 1, 2, ..., n 的所有极点在正方向上恰好绕一次，但不绕 Γ(bj − s), j = 1, 2, ..., m 的任何极点。
我如何得到 L 并解决积分？我习惯计算实数积分的方法是
import numpy as np
myL = np.linspace(0, 1, 100)
densityL = myL[1] - myL[0]
myIntegral = (F(myL)*densityL).sum()
我并不太追求效率，我更喜欢一个简单而缓慢的工作示例，以便于我理解方法论。
