1 | LIB "tst.lib"; tst_init(); |
---|
2 | LIB "resbinomial.lib"; |
---|
3 | ring r = 0,(x(1..4)),dp; |
---|
4 | list flag=identifyvar(); |
---|
5 | ideal J=x(1)^2-x(2)^2*x(3)^5, x(1)*x(3)^3+x(4)^6; |
---|
6 | list Lmb=1,list(0,0,0,0),list(0,0,0,0),list(0,0,0,0),iniD(4),iniD(4),list(0,0,0,0),-1; |
---|
7 | list L=data(J,2,4); |
---|
8 | list LL=determinecenter(L[1],L[2],2,4,0,0,Lmb,flag,0,-1); // Compute the first center |
---|
9 | LL[1]; // index of variables in the center |
---|
10 | LL[2]; // exponents of ideals J_4,J_3,J_2,J_1 |
---|
11 | LL[3]; // list of orders of J_4,J_3,J_2,J_1 |
---|
12 | LL[4]; // list of critical values |
---|
13 | LL[5]; // components of the resolution function t |
---|
14 | LL[6]; // list of D_4,D_3,D_2,D_1 |
---|
15 | LL[7]; // list of H_4,H_3,H_2,H_1 (exceptional divisors) |
---|
16 | LL[8]; // list of all exceptional divisors acumulated |
---|
17 | LL[9]; // auxiliary invariant |
---|
18 | LL[10]; // intvec pointing out the last step where the function t has dropped |
---|
19 | ring r= 0,(x(1..4)),dp; |
---|
20 | list flag=identifyvar(); |
---|
21 | ideal J=x(1)^3-x(2)^2*x(3)^5, x(1)*x(3)^3+x(4)^5; |
---|
22 | list Lmb=2,list(0,0,0,0),list(0,0,0,0),list(0,0,0,0),iniD(4),iniD(4),list(0,0,0,0),-1; |
---|
23 | list L2=data(J,2,4); |
---|
24 | list L3=determinecenter(L2[1],L2[2],2,4,0,0,Lmb,flag,0,-1); // Example with rational exponents in E-Coeff |
---|
25 | L3[1]; // index of variables in the center |
---|
26 | L3[2]; // exponents of ideals J_4,J_3,J_2,J_1 |
---|
27 | L3[3]; // list of orders of J_4,J_3,J_2,J_1 |
---|
28 | L3[4]; // list of critical values |
---|
29 | L3[5]; // components of the resolution function |
---|
30 | tst_status(1);$ |
---|