diff --git a/src/gamma_inc.jl b/src/gamma_inc.jl
index 5a8b8f5c..14994207 100644
--- a/src/gamma_inc.jl
+++ b/src/gamma_inc.jl
@@ -425,7 +425,7 @@ function gamma_inc_taylor(a::Float64, x::Float64, ind::Integer)
     apn = a
 
     # compute and store larger terms in wk, to add from small to large
-    t = 1
+    t = 1.0
     i = 0
     while i < 20
         i += 1
@@ -445,8 +445,8 @@ function gamma_inc_taylor(a::Float64, x::Float64, ind::Integer)
     end
 
     # sum terms from small to large
-    for v ∈ @view wk[i:-1:1]
-        sm += v
+    for j ∈ i:(-1):1
+        sm += wk[j]
     end
     
     p = (rgammax(a, x) / a) * (1.0 + sm)