diff --git a/plots/solutions1.png b/plots/solutions1.png
index bbea945..ee6368b 100644
Binary files a/plots/solutions1.png and b/plots/solutions1.png differ
diff --git a/plots/solutions2.png b/plots/solutions2.png
index 07aeb3e..68e3e23 100644
Binary files a/plots/solutions2.png and b/plots/solutions2.png differ
diff --git a/plots/solutions3.png b/plots/solutions3.png
index 619ab47..fd66206 100644
Binary files a/plots/solutions3.png and b/plots/solutions3.png differ
diff --git a/solve.jl b/solve.jl
index c72eec1..650bd8e 100644
--- a/solve.jl
+++ b/solve.jl
@@ -20,7 +20,7 @@ function solve(F, (G, roots) = start_system(F), maxsteps=10000)
   #  F=homogenize(F)
   H=homotopy(F,G)
   solutions = []
-  steps = 0
+  step_array = []
 
   @time Threads.@threads for r in roots
     t = 1.0
@@ -37,9 +37,10 @@ function solve(F, (G, roots) = start_system(F), maxsteps=10000)
       steps += 1
     end
     push!(solutions, x0)
+    push!(step_array, steps)
   end
 
-  return (solutions, steps)
+  return (solutions, step_array)
 end
 
 # Input polynomial system