scipy.special.rgamma

scipy.special.rgamma(x[, out])

y=rgamma(z) returns one divided by the gamma function of x.