diff --git a/Project.toml b/Project.toml
index a614e23..b79419d 100644
--- a/Project.toml
+++ b/Project.toml
@@ -1,7 +1,7 @@
 name = "OrbitalTrajectories"
 uuid = "2b613a20-8d2a-5290-b19f-e06f4bcc2e7d"
 authors = ["Dan Padilha"]
-version = "0.1.13"
+version = "0.1.14"
 
 [deps]
 DataStructures = "864edb3b-99cc-5e75-8d2d-829cb0a9cfe8"
diff --git a/src/dynamics/dynamical_models.jl b/src/dynamics/dynamical_models.jl
index 68e2ea8..f3eca7a 100644
--- a/src/dynamics/dynamical_models.jl
+++ b/src/dynamics/dynamical_models.jl
@@ -57,7 +57,7 @@ has_jacobian(X::Type{<:DiffEqBase.ODEFunction}) = !isnothing(fieldtype(X, :jac))
     # The State Transition Matrix (STM) ODE function is defined as follows, including the N^2 Jacobian equations +
     # the N first-order equations of motion. [Koon 2011]
     # NOTE: the Differential is defined element-wise and flattened to a list.
-    stm_eqs = simplify.(D.(ϕ) .~ A * ϕ)
+    stm_eqs = collect(simplify.(D.(ϕ) .~ A * ϕ))
 
     # Create the ODE system and generate its functions
     stm_ode = ODESystem(