diff --git a/Project.toml b/Project.toml
index 14995750..e0e20083 100644
--- a/Project.toml
+++ b/Project.toml
@@ -1,7 +1,7 @@
 name = "HarmonicBalance"
 uuid = "e13b9ff6-59c3-11ec-14b1-f3d2cc6c135e"
 authors = ["Jan Kosata <kosataj@phys.ethz.ch>", "Javier del Pino <jdelpino@phys.ethz.ch>", "Orjan Ameye <orjan.ameye@hotmail.com>"]
-version = "0.10.7"
+version = "0.10.8"
 
 [deps]
 BijectiveHilbert = "91e7fc40-53cd-4118-bd19-d7fcd1de2a54"
@@ -48,7 +48,7 @@ HomotopyContinuation = "2.9"
 JET = "0.9.9"
 JLD2 = "0.4.48, 0.5"
 Latexify = "0.16"
-ModelingToolkit = "9.34"
+ModelingToolkit = "9.34 - 9.40"
 NonlinearSolve = "3.14"
 OrderedCollections = "1.6"
 OrdinaryDiffEq = "v6.89"
@@ -56,7 +56,7 @@ Peaks = "0.5"
 Plots = "1.39"
 PrecompileTools = "1.2"
 ProgressMeter = "1.7.2"
-SteadyStateDiffEq = "2.3"
+SteadyStateDiffEq = "2.3.2"
 SymbolicUtils = "3.5"
 Symbolics = "6.4"
 julia = "1.10.0"
diff --git a/ext/ModelingToolkitExt.jl b/ext/ModelingToolkitExt.jl
index 932deb2d..1866f7ed 100644
--- a/ext/ModelingToolkitExt.jl
+++ b/ext/ModelingToolkitExt.jl
@@ -14,7 +14,8 @@ using ModelingToolkit:
     varmap_to_vars,
     parameters,
     @parameters,
-    @mtkbuild
+    @mtkbuild,
+    @independent_variables
 
 swapsides(eq::Equation) = Equation(eq.rhs, eq.lhs)
 
@@ -30,7 +31,7 @@ function ModelingToolkit.ODESystem(eom::HarmonicEquation)
         eom = rearrange_standard(eom)
     end
 
-    slow_time = (@variables T; T)
+    slow_time = (@independent_variables T; T)
     par_names = declare_parameter.(eom.parameters)
     vars = get_variables(eom)