I’m using ModelingToolkit v9.30.1:
[c3fe647b] AbstractAlgebra v0.42.0
[0f109fa4] BifurcationKit v0.3.6
[0c46a032] DifferentialEquations v7.13.0
[23fbe1c1] Latexify v0.16.4
[961ee093] ModelingToolkit v9.30.1
[8913a72c] NonlinearSolve v3.13.1
[91a5bcdd] Plots v1.40.5
[f27b6e38] Polynomials v4.0.11
[438e738f] PyCall v1.96.4
[9672c7b4] SteadyStateDiffEq v2.3.0
[24249f21] SymPy v2.1.1
[0c5d862f] Symbolics v5.35.1
Like I said in the original post, I’ve been using an older version before (v8.36.0) because of the requirements in FindSteadyStates.jl but updated to the latest version after seeing this post where an update solved the problem that t_nounits and D_nounits couldn’t be imported. In my case, the update didn’t solve the problems.