diff --git a/test/sympy.jl b/test/sympy.jl
index d768a18a1..030ff6d77 100644
--- a/test/sympy.jl
+++ b/test/sympy.jl
@@ -21,6 +21,8 @@ expr = x * p + (x^2 - 1 + y) * (p + 2t)
 sexpr = symbolics_to_sympy(expr)
 sp = symbolics_to_sympy(p)
 
+@test SymPy.simplify(symbolics_to_sympy(Symbolics.solve_for(expr, p))) == SymPy.simplify(SymPy.solve(sexpr, sp)[1])
+ 
 symbolics_sol = SymPy.simplify(symbolics_to_sympy(Symbolics.symbolic_linear_solve(expr, p)))
 sympy_sols = SymPy.solve(SymPy.expand(sexpr), sp)
 @test !isempty(sympy_sols) && isequal(symbolics_sol, sympy_sols[1])
@@ -73,4 +75,4 @@ sol_vars = Symbolics.get_variables(sol_ode)
 const_sym = only(filter(v -> startswith(string(Symbolics.nameof(v)), "C"), sol_vars))
 expected_sol = C1 * exp(2 * x)
 canonical_sol_ode = Symbolics.substitute(sol_ode, Dict(const_sym => C1))
-@test isequal(canonical_sol_ode, expected_sol)
+@test isequal(canonical_sol_ode, expected_sol)
\ No newline at end of file