scipy.special.rgamma

scipy.special.rgamma(x[, out]) = <ufunc 'rgamma'>

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