diff --git a/docs/Project.toml b/docs/Project.toml
index 075f32e31..aa6ae1a41 100644
--- a/docs/Project.toml
+++ b/docs/Project.toml
@@ -12,3 +12,6 @@ OrdinaryDiffEq = "1dea7af3-3e70-54e6-95c3-0bf5283fa5ed"
 Plots = "91a5bcdd-55d7-5caf-9e0b-520d859cae80"
 ProbNumDiffEq = "bf3e78b0-7d74-48a5-b855-9609533b56a5"
 Statistics = "10745b16-79ce-11e8-11f9-7d13ad32a3b2"
+
+[compat]
+ModelingToolkit = "9"
diff --git a/docs/src/tutorials/dae.md b/docs/src/tutorials/dae.md
index 544df1d0c..dcbe7df6a 100644
--- a/docs/src/tutorials/dae.md
+++ b/docs/src/tutorials/dae.md
@@ -115,7 +115,7 @@ simplified_sys = structural_simplify(dae_index_lowering(traced_sys))
 
 Let's build two different ODE problems, and check how well we can solve each:
 ```@example dae
-prob_index3 = ODEProblem(traced_sys, Pair[], tspan)
+prob_index3 = ODEProblem(complete(traced_sys), Pair[], tspan)
 prob_index1 = ODEProblem(simplified_sys, Pair[], tspan)
 
 sol3 = solve(prob_index3, EK1())