multivariate polynomial division, groebner bases; using sympy