| LIB "equising.lib";
int p=printlevel;
printlevel=1;
ring r = 0,(a,b,c,d,e,f,g,x,y),ds;
poly F = (x2+2xy+y2+x5)+ax+by+cx2+dxy+ey2+fx3+gx4;
list M = esStratum(F);
M[1];
→ _[1]=g
→ _[2]=f
→ _[3]=b
→ _[4]=a
→ _[5]=-4c+4d-4e+d2-4ce
printlevel=3; // displays additional information
esStratum(F,2); // es stratum over Q[a,b,c,d,e,f,g] / <a,b,c,d,e,f,g>^2
→ //
→ // Compute HN development
→ // ----------------------
→ // finished
→ //
→ // Blowup Step 1 completed
→ // Blowup Step 2 completed
→ // Blowup Step 3 completed
→ // 1 branch finished
→ //
→ // Elimination starts:
→ // -------------------
→ // finished
→ //
→ // output of 'esStratum' is list consisting of:
→ // _[1] = ideal defining equisingular stratum
→ // _[2] = 0
→ [1]:
→ _[1]=b
→ _[2]=a
→ _[3]=c-d+e
→ _[4]=g
→ _[5]=f
→ [2]:
→ 0
ideal I = f-fa,e+b;
qring q = std(I);
poly F = imap(r,F);
esStratum(F);
→ //
→ // Compute HN development
→ // ----------------------
→ // finished
→ //
→ // Blowup Step 1 completed
→ // Blowup Step 2 completed
→ // Blowup Step 3 completed
→ // 1 branch finished
→ //
→ // Elimination starts:
→ // -------------------
→ // finished
→ //
→ // output of 'esStratum' is list consisting of:
→ // _[1] = ideal defining equisingular stratum
→ // _[2] = 0
→ [1]:
→ _[1]=e
→ _[2]=a
→ _[3]=-4c+4d+d2
→ _[4]=g
→ [2]:
→ 0
printlevel=p;
|