numpy - Meijer G-function in Python and scipy -

i'm in need of meijer g function in scipy. read somewhere on internet due generality, meier g function not supported special function in scipy, should write according personal use case.

my problem have no experience whatsoever complex integration. latex forbidden here, here's i'm trying solve numerically:

(the first line being general case, second line case i'm trying compute), p(a), k, k2 given

as wikipedia states, there 3 ways l:

• l runs −i∞ +i∞ such poles of Γ(bj − s), j = 1, 2, ..., m, on right of path, while poles of Γ(1 − ak + s), k = 1, 2, ..., n, on left.
• l loop beginning , ending @ +∞, encircling poles of Γ(bj − s), j = 1, 2, ..., m, once in negative direction, not encircling pole of Γ(1 − ak + s), k = 1, 2, ..., n.
• l loop beginning , ending @ −∞ , encircling poles of Γ(1 − ak + s), k = 1, 2, ..., n, once in positive direction, not encircling pole of Γ(bj − s), j = 1, 2, ..., m.

how l , solve integral? way i'm used compute integrals on reals to

import numpy np myl = np.linspace(0, 1, 100) densityl = myl[1] - myl[0] myintegral = (f(myl)*densityl).sum() 

i'm not after efficiency, i'd prefer simple , slow working example can use understand methodology.

for complicated, think should avoid calculating integral yourself, if don't have experience complex integration, , use tested existing implementation.

meijer g-function implemented in mpmath , possibly in sympy.