Singular

LIB "curveInv.lib";
//////////////////////////////////////////////
// Mutltiplicity of the conductor of curves //
//////////////////////////////////////////////
ring R = 0,(x,y,z),ds;
// Example 1:
ideal I = x2-y4z,z3y2+xy2;
I = std(radical(I));
curveConductorMult(I);
==> 23
// Example 2:
ideal I = x*(y+z)^3 - y3, x2y2 + z5;
==> // ** redefining I (ideal I = x*(y+z)^3 - y3, x2y2 + z5;)
I = std(radical(I));
curveConductorMult(I);
==> 19
kill R;
////////////////////////////////////////////////////////
// Mutltiplicity of the conductor of Gorenstein curve //
////////////////////////////////////////////////////////
ring R = 0,(x,y),ds;
ideal I = xy;
// In such a case, the conductor multiplicity c satisfies: c = 2*delta
// Delta invariant:
curveDeltaInv(I);
==> 1
// Conductor Multiplicity:
curveConductorMult(I);
==> 2