fix typos

Author: oscardssmith

src/gamma.jl

@@ -17,41 +17,41 @@ function gamma(_x::Float64)
         w = muladd(w, evalpoly(w, s), 1.0)
         # avoid overflow
         v = x ^ muladd(0.5, x, -0.25)
-
         res = SQ2PI(T) * v * (v / exp(x)) * w
-	
-    	if _x < 0
-	    return π / (res * s)
-	else
-	    return res
-	end
+
+        if _x < 0
+            return π / (res * s)
+        else
+            return res
+        end
     end
     P = (
-	1.000000000000000000009e0, 8.378004301573126728826e-1, 3.629515436640239168939e-1, 1.113062816019361559013e-1,
-	2.385363243461108252554e-2, 4.092666828394035500949e-3, 4.542931960608009155600e-4, 4.212760487471622013093e-5
+        1.000000000000000000009e0, 8.378004301573126728826e-1, 3.629515436640239168939e-1, 1.113062816019361559013e-1,
+        2.385363243461108252554e-2, 4.092666828394035500949e-3, 4.542931960608009155600e-4, 4.212760487471622013093e-5
     )
     Q = (
-	9.999999999999999999908e-1, 4.150160950588455434583e-1, -2.243510905670329164562e-1, -4.633887671244534213831e-2,
-	2.773706565840072979165e-2, -7.955933682494738320586e-4, -1.237799246653152231188e-3, 2.346584059160635244282e-4,
-	-1.397148517476170440917e-5
+        9.999999999999999999908e-1, 4.150160950588455434583e-1, -2.243510905670329164562e-1, -4.633887671244534213831e-2,
+        2.773706565840072979165e-2, -7.955933682494738320586e-4, -1.237799246653152231188e-3, 2.346584059160635244282e-4,
+        -1.397148517476170440917e-5
     )
 
     z = 1.0
     while x >= 3.0
-	x -= 1.0
-	z *= x
+        x -= 1.0
+        z *= x
     end
     while x < 2.0
-	z /= x
-	x += 1.0
+        z /= x
+        x += 1.0
     end
 
     x -= 2.0
     p = evalpoly(x, P)
     q = evalpoly(x, Q)
-    retrun _x < 0 ? π * q / (s * z * p) : z * p / q
+    return _x < 0 ? π * q / (s * z * p) : z * p / q
 end
 
+
 function gamma(_x::Float32)
     x = Float64(_x)
     if _x < 0