| LIB "arnold.lib";
ring R = 0,(x,y,z),ds;
poly g = x^2*y^4+y^2*x^4+5*x^(10)+y^(11)+z^2;
poly phix = x+y^2+x^2+x*y+x^2*y+x*y^3;
poly phiy = y+y^2+2*x^2+x*y+y*x^2+y^2*x+x*y^4;
poly phiz = z+2x+x^2+y^4*x;
map phi = R,phix,phiy,phiz;
g = phi(g);
Poly F = makePoly(g);
NormalForm N= determineNormalForm(F);
N=determineExceptionalHypersurface(N);
exceptionalHypersurface(N);
==> a(4)*a(5)^2*a(6)-4*a(4)^2*a(6)^2
==>
|