diff --git a/docs/Project.toml b/docs/Project.toml
index 1a028e65a..1eb2a6513 100644
--- a/docs/Project.toml
+++ b/docs/Project.toml
@@ -60,7 +60,7 @@ FiniteDiff = "2"
 ForwardDiff = "0.10, 1"
 IncompleteLU = "0.2"
 JLD2 = "0.4, 0.5.1"
-ModelingToolkit = "10"
+ModelingToolkit = "9, 10"
 NonlinearSolve = "3.15, 4"
 ODEProblemLibrary = "0.1"
 Optimization = "3, 4"
diff --git a/docs/src/examples/outer_solar_system.md b/docs/src/examples/outer_solar_system.md
index 207c9261f..6643f37be 100644
--- a/docs/src/examples/outer_solar_system.md
+++ b/docs/src/examples/outer_solar_system.md
@@ -71,7 +71,7 @@ Thus, $\dot{q}$ is defined by the masses. We only need to define $\dot{p}$, and
 ```@example outersolarsystem
 eqs = vec(@. D(D(u))) .~ .-ModelingToolkit.gradient(potential, vec(u)) ./
                          repeat(M, inner = 3)
-@mtkcompile sys = System(eqs, t)
+@mtkbuild sys = System(eqs, t)
 prob = ODEProblem(ss, [vec(u .=> pos); vec(D.(u) .=> vel)], tspan)
 sol = solve(prob, Tsit5());
 ```