def _handle_positive_dimensional(polys, symbols, denominators):
    from sympy.polys.polytools import groebner
    _symbols = list(symbols)
    _symbols.sort(key=default_sort_key)
    basis = groebner(polys, _symbols, polys=True)
    new_system = []
    for poly_eq in basis:
        new_system.append(poly_eq.as_expr())
    result = [{}]
    result = substitution(new_system, symbols, result, [], denominators)
    return result
def _handle_positive_dimensional(polys, symbols, denominators):
    from sympy.polys.polytools import groebner
    _symbols = list(symbols)
    _symbols.sort(key=default_sort_key)
    basis = groebner(polys, _symbols, polys=True)
    new_system = []
    for poly_eq in basis:
        new_system.append(poly_eq.as_expr())
    result = [{}]
    result = substitution(new_system, symbols, result, [], denominators)
    return result