sagbigrob.lib
LIB "sagbigrob.lib"; ring r=0,(x,y),Dp; ideal A=x2,xy; ideal I=x3y+x2,x4+x2y2,-x3y3-x2y2; poly f=x6y6+x2; poly g=x12+x6y6-x8+x5y; SGNF(f,I,A); ==> x2y2+x2 SGNF(g,I,A); ==> 0