Getting the symbolic solution of system of equations using ModelingToolKit

Yes, I run “using Groebner”, and I can see it’s loaded by:

using Pkg
Pkg.status()
Status `~/.julia/environments/v1.11/Project.toml`
  [6e4b80f9] BenchmarkTools v1.5.0
  [a93c6f00] DataFrames v1.7.0
⌃ [82cc6244] DataInterpolations v6.5.2
⌃ [0c46a032] DifferentialEquations v7.14.0
  [0b43b601] Groebner v0.8.2
⌃ [7073ff75] IJulia v1.25.0
  [a98d9a8b] Interpolations v0.15.1
⌃ [961ee093] ModelingToolkit v9.49.0
  [16a59e39] ModelingToolkitStandardLibrary v2.17.0
⌅ [8913a72c] NonlinearSolve v3.15.1
⌃ [1dea7af3] OrdinaryDiffEq v6.89.0
⌃ [f0f68f2c] PlotlyJS v0.18.14
⌃ [91a5bcdd] Plots v1.40.8
  [24249f21] SymPy v2.2.0
  [2efcf032] SymbolicIndexingInterface v0.3.34

The full error stack is:

┌ Info: Assuming (R2*RisoN + R2*RisoP + R3*RisoN + R3*RisoP + RisoN*RisoP) != 0
└ @ Symbolics /home/ricardo/.julia/packages/Symbolics/6CYZh/src/solver/preprocess.jl:53
┌ Info: Assuming (R1*R2*RisoN + R1*R2*RisoP + R1*R3*RisoN + R1*R3*RisoP + R1*RisoN*RisoP + R2*RisoN*RisoP + R3*RisoN*RisoP) != 0
└ @ Symbolics /home/ricardo/.julia/packages/Symbolics/6CYZh/src/solver/preprocess.jl:53

"Groebner bases engine is required. Execute `using Groebner` to enable this functionality."

Stacktrace:
 [1] solve_multivar(eqs::Vector{Num}, vars::Tuple{Num, Num}; dropmultiplicity::Bool, warns::Bool)
   @ Symbolics ~/.julia/packages/Symbolics/6CYZh/src/solver/main.jl:328
 [2] symbolic_solve(expr::Vector{Equation}, x::Tuple{Num, Num}; dropmultiplicity::Bool, warns::Bool)
   @ Symbolics ~/.julia/packages/Symbolics/6CYZh/src/solver/main.jl:204
 [3] symbolic_solve(expr::Vector{Equation}, x::Tuple{Num, Num})
   @ Symbolics ~/.julia/packages/Symbolics/6CYZh/src/solver/main.jl:145
 [4] top-level scope
   @ ~/jupyter/isolation_sensing/jl_notebook_cell_df34fa98e69747e1a8f8a730347b8e2f_W1sZmlsZQ==.jl:1