y=rgamma(z) returns one divided by the gamma function of x.
scipy.special.psi
scipy.special.polygamma
Enter search terms or a module, class or function name.