See https://github.com/sympy/sympy/issues/5031. We could implement the solution given there in `Equation()`