diff --git a/README.md b/README.md
index 6f4711f..d7d4d8a 100644
--- a/README.md
+++ b/README.md
@@ -31,6 +31,7 @@ the documentation, which contains the unreleased features.
 
 ```julia
 using ModelingToolkit, MethodOfLines, DifferentialEquations, ModelOrderReduction
+using SciMLBase: symbolic_discretize
 
 # firstly construct a ModelingToolkit.PDESystem for the FitzHugh-Nagumo model
 @variables x t v(..) w(..)
diff --git a/docs/Project.toml b/docs/Project.toml
index c7a4d8f..9796956 100644
--- a/docs/Project.toml
+++ b/docs/Project.toml
@@ -7,6 +7,7 @@ MethodOfLines = "94925ecb-adb7-4558-8ed8-f975c56a0bf4"
 ModelOrderReduction = "207801d6-6cee-43a9-ad0c-f0c64933efa0"
 ModelingToolkit = "961ee093-0014-501f-94e3-6117800e7a78"
 Plots = "91a5bcdd-55d7-5caf-9e0b-520d859cae80"
+SciMLBase = "0bca4576-84f4-4d90-8ffe-ffa030f20462"
 
 [compat]
 DifferentialEquations = "8"
@@ -16,3 +17,4 @@ MethodOfLines = "0.11"
 ModelOrderReduction = "0.1"
 ModelingToolkit = "11"
 Plots = "1.40"
+SciMLBase = "3"
diff --git a/docs/src/tutorials/deim_FitzHugh-Nagumo.md b/docs/src/tutorials/deim_FitzHugh-Nagumo.md
index 93e75e2..af11ecb 100644
--- a/docs/src/tutorials/deim_FitzHugh-Nagumo.md
+++ b/docs/src/tutorials/deim_FitzHugh-Nagumo.md
@@ -73,6 +73,7 @@ Next, we apply finite difference discretization using
 
 ```@example deim_FitzHugh_Nagumo
 using MethodOfLines
+using SciMLBase: symbolic_discretize
 N = 15 # equidistant discretization intervals
 dx = (L - 0.0) / N
 dxs = [x => dx]
diff --git a/test/Project.toml b/test/Project.toml
index 6a2d831..9a517fd 100644
--- a/test/Project.toml
+++ b/test/Project.toml
@@ -4,6 +4,7 @@ MethodOfLines = "94925ecb-adb7-4558-8ed8-f975c56a0bf4"
 ModelingToolkit = "961ee093-0014-501f-94e3-6117800e7a78"
 OrdinaryDiffEq = "1dea7af3-3e70-54e6-95c3-0bf5283fa5ed"
 SafeTestsets = "1bc83da4-3b8d-516f-aca4-4fe02f6d838f"
+SciMLBase = "0bca4576-84f4-4d90-8ffe-ffa030f20462"
 SciMLTesting = "09d9d899-5365-40a9-917a-5f67fddea283"
 Test = "8dfed614-e22c-5e08-85e1-65c5234f0b40"
 
@@ -13,4 +14,5 @@ MethodOfLines = "0.11"
 ModelingToolkit = "11"
 OrdinaryDiffEq = "7"
 SafeTestsets = "0.1"
+SciMLBase = "3"
 SciMLTesting = "2.4"
diff --git a/test/deim.jl b/test/deim.jl
index 7868a31..436cb3c 100644
--- a/test/deim.jl
+++ b/test/deim.jl
@@ -1,5 +1,6 @@
 using Test, ModelOrderReduction
 using ModelingToolkit, MethodOfLines, OrdinaryDiffEq
+using SciMLBase: symbolic_discretize
 
 # construct an ModelingToolkit.ODESystem with non-empty field substitutions
 @variables x t v(..) w(..)