bug fix to polygamma for order > 1, found by @andreasnoackjensen in #…

…7033
JuliaLang · Jun 3, 2014 · 4ba3b7b · 4ba3b7b · JeffBezanson · Jun 3, 2014
1 parent b80f270
commit 4ba3b7b
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/base/special/gamma.jl b/base/special/gamma.jl
@@ -308,7 +308,7 @@ function psifn(x::Float64, n::Int, kode::Int, m::Int)
     ans[1] = s - log(xq)
     return ans
 end
-polygamma(k::Int, x::Float64) = (2rem(k,2) - 1)*psifn(x, k, 1, 1)[1]/gamma(k + 1)
+polygamma(k::Int, x::Float64) = (2rem(k,2) - 1)*psifn(x, k, 1, 1)[1]*gamma(k + 1)
 polygamma(k::Int, x::Float32) = float32(polygamma(k, float64(x)))
 polygamma(k::Int, x::Real) = polygamma(k, float64(x))