Revision history [back]

Here's a shot at it. You can call some maxima commands.

var('x,y')
maxima('load(grobner)')
maxima('poly_monomial_order:grlex')
F=[x*y^2-x , x-y^3]
ans=maxima('poly_pseudo_divide(x^2*y^2+x^3*y^2-y+1,[x*y^2-x,x-y^3],[x,y])')
print ans
(quo,rem,n,m)=ans
p=quo[0]*F[0]+quo[1]*F[1]+rem
print p
print p.expand()