public void ValidateDensityLn(double a, double b, double x)
        {
            var    n        = new InverseGamma(a, b);
            double expected = Math.Log(Math.Pow(b, a) * Math.Pow(x, -a - 1.0) * Math.Exp(-b / x) / SpecialFunctions.Gamma(a));

            Assert.AreEqual(expected, n.DensityLn(x));
            Assert.AreEqual(expected, InverseGamma.PDFLn(a, b, x));
        }