read recden; read NGCD; read PGCD; read AGCD; # This test hits the zero divisor y g := rpoly( x^3+y*x^2+1, [x,y], y^3 ); f := rpoly( x^4+y*x+1, [x,y], y^3 ); h := rpoly( x^5-y*x^3+y, [x,y], y^3 ); F := mulrpoly(f,g); H := mulrpoly(h,g); gcdrpoly(F,H); NGCD(F,H); # Don't try to split y^3-0 because the roots are not separable g := rpoly( x^3+(y-2)*x^2+1, [x,y], y^3 ); f := rpoly( x^2+(y+31)*x+1, [x,y], y^3 ); h := rpoly( x^3-(11-y)*x^3+31, [x,y], y^3 ); F := mulrpoly(f,g); H := mulrpoly(h,g); gcdrpoly(F,H); NGCD(F,H);