PolyGamma[r]= e^-(EulerGamma r)