diff --git a/src/Numerics/Distributions/Dirichlet.cs b/src/Numerics/Distributions/Dirichlet.cs index 57eaf8b8..295d5fec 100644 --- a/src/Numerics/Distributions/Dirichlet.cs +++ b/src/Numerics/Distributions/Dirichlet.cs @@ -227,8 +227,8 @@ namespace MathNet.Numerics.Distributions { get { - var num = _alpha.Sum(t => (t - 1)*SpecialFunctions.DiGamma(t)); - return SpecialFunctions.GammaLn(AlphaSum) + ((AlphaSum - Dimension)*SpecialFunctions.DiGamma(AlphaSum)) - num; + var num = _alpha.Sum(t => (t - 1)*SpecialFunctions.DiGamma(t) - SpecialFunctions.GammaLn(t)); + return -SpecialFunctions.GammaLn(AlphaSum) + ((AlphaSum - Dimension)*SpecialFunctions.DiGamma(AlphaSum)) - num; } }