diff --git a/Project.toml b/Project.toml
index ae4e8a6..b151eff 100644
--- a/Project.toml
+++ b/Project.toml
@@ -1,7 +1,7 @@
 name = "DifferentiableStateSpaceModels"
 uuid = "beacd9db-9e5e-4956-9b09-459a4b2028df"
 authors = ["Jesse Perla <jesseperla@gmail.com> and contributors"]
-version = "0.5.1"
+version = "0.5.2"
 
 [deps]
 ChainRulesCore = "d360d2e6-b24c-11e9-a2a3-2a2ae2dbcce4"
diff --git a/src/generate_perturbation.jl b/src/generate_perturbation.jl
index 99797ff..b6d5860 100644
--- a/src/generate_perturbation.jl
+++ b/src/generate_perturbation.jl
@@ -28,7 +28,7 @@ function verify_steady_state(m, p_d, p_f; atol = 1e-8, args...)
     p_d_symbols = collect(Symbol.(keys(p_d)))
     p = order_vector_by_symbols(merge(p_d, p_f), m.mod.m.p_symbols)
 
-    w = hcat(sol.y, sol.x) # get a vector for the proposed steadystate
+    w = vcat(sol.y, sol.x) # get a vector for the proposed steadystate
     H = Vector{Float64}(undef, length(w)) # allocate it, but leave undef to make sure we can see if it goes to 0 or not
     m.mod.m.H̄!(H, w, p)  # evaluate in place
     return (norm(H) < atol)