private static double PositiveLGamma(double v0) { double vg = v0 + 12.644565; // v0 + g - 0.5 double res = Math.Log(GammaRationalFunc(v0)) - vg; res += (v0 - 0.5) * Math.Log(vg); return res; }