Python SciPy特殊功能gammaincc和gammainccinv不匹配

如何解决Python SciPy特殊功能gammaincc和gammainccinv不匹配

我正在尝试从一个较高的不完全伽马函数生成随机数，并且我意识到SciPy中实现的反不完全伽马函数gammainccinv与gammaincc并不完全对应。以下代码可用于说明差异：

import numpy
from scipy.stats import poisson
from scipy.special import gammainc,gammaincc
from scipy.special import gammaincinv,gammainccinv
import matplotlib.pyplot as plt

lam = 10
xr = numpy.arange(0,20,0.01)
xri = numpy.arange(20)
y = numpy.arange(0,1,0.001)

plt.plot(xri,poisson.pmf(xri,lam),"*",label="Poisson")
plt.plot(xr-0.5,(gammaincc(xr+1e-6,lam)-gammaincc(xr,lam))/1e-6,label="Upper incomplete gamma")
x = gammainccinv(lam,y)
plt.plot(0.5*(x[1:]+x[:-1]),(y[:-1]-y[1:])/(x[1:]-x[:-1]),label="Inverse upper incomplete gamma\n(transformed)")
plt.legend(fontsize=8)
plt.show()

在Python中是否存在逆上不完全伽马函数的替代实现？