| LIB "normal.lib";
ring r=0,(x,y),dp;
poly f=(x2+y^2-1)^3 +27x2y2;
ideal I=f,jacob(f);
I=std(I);
list qr=minAssGTZ(I);
size(qr);
==> 6
// each component of the singular locus either describes a cusp or a pair
// of conjugated nodes:
deltaLoc(f,qr[1]);
==> [1]:
==> 1
==> [2]:
==> 2
==> [3]:
==> 1
deltaLoc(f,qr[2]);
==> [1]:
==> 1
==> [2]:
==> 2
==> [3]:
==> 1
deltaLoc(f,qr[3]);
==> [1]:
==> 1
==> [2]:
==> 2
==> [3]:
==> 1
deltaLoc(f,qr[4]);
==> [1]:
==> 1
==> [2]:
==> 2
==> [3]:
==> 1
deltaLoc(f,qr[5]);
==> [1]:
==> 2
==> [2]:
==> 2
==> [3]:
==> 4
deltaLoc(f,qr[6]);
==> [1]:
==> 2
==> [2]:
==> 2
==> [3]:
==> 4
|