Equation not in exact domain. Try converting to rational

Raise code

                result = substitution(
                    polys_expr, symbols, exclude=denominators)
                return result

        # positive dimensional system
        res = _handle_positive_dimensional(polys, symbols, denominators)
        if res is EmptySet and any(not p.domain.is_Exact for p in polys):
            raise NotImplementedError("Equation not in exact domain. Try converting to rational")
        else:
            return res

    else:
        # If all the equations are not polynomial.
        # Use `substitution` method for the system
        result = substitution(