diff --git a/README.md b/README.md
index d93a18e3f7..12570c35a3 100644
--- a/README.md
+++ b/README.md
@@ -36,7 +36,7 @@ lower it to a first order system, symbolically generate the Jacobian function
 for the numerical integrator, and solve it.
 
 ```julia
-using DifferentialEquations, ModelingToolkit
+using OrdinaryDiffEqDefault, ModelingToolkit
 using ModelingToolkit: t_nounits as t, D_nounits as D
 
 @parameters σ ρ β