typo

Author: oscardssmith

src/gamma.jl

@@ -9,12 +9,12 @@ function gamma(_x::Float64)
     end
     if x > 11.5
         w = inv(x)
-        s = (
+        coefs = (1.0, 
             8.333333333333331800504e-2, 3.472222222230075327854e-3, -2.681327161876304418288e-3, -2.294719747873185405699e-4,
             7.840334842744753003862e-4, 6.989332260623193171870e-5, -5.950237554056330156018e-4, -2.363848809501759061727e-5,
             7.147391378143610789273e-4
         )
-        w = muladd(w, evalpoly(w, s), 1.0)
+        w = evalpoly(w, coefs)
         # avoid overflow
         v = x ^ muladd(0.5, x, -0.25)
         res = SQ2PI(T) * v * (v / exp(x)) * w