def _pdf(self, x, df):
        r = asarray(df*1.0)
        Px = exp(gamln((r+1)/2)-gamln(r/2))
        Px /= sqrt(r*pi)*(1+(x**2)/r)**((r+1)/2)
        return Px
dataset/ETHPy150Open scipy/scipy/scipy/stats/_continuous_distns.py/t_gen._pdf/original