不完整的gamma函数不接受复杂的值输入

我正在计算 upper incomplete gamma function 定义类似于 this post . 如果我使用

from scipy.special import gamma,gammainc
from numpy import linspace

a = 0
z = (2+3j)*np.linspace(0,10)
gamma(a)*(1-gammainc(a,z))

哪里 z 是一个复杂的向量,我得到一个错误

TypeError: ufunc 'gammainc' not supported for the input types, and the inputs could not be safely coerced to any supported types according to the casting rule ''safe''

是否有其他函数来进行计算?当我试着用Wolframalpha做这个的时候,似乎没有问题。 Gamma 功能。