fieker-DuValspielwiese
Last change
on this file since 3fec5d0 was
3a823d,
checked in by Hans Schoenemann <hannes@…>, 9 years ago
|
fixes by pfister: tr #634, tr #682
|
-
Property mode set to
100644
|
File size:
669 bytes
|
Line | |
---|
1 | LIB "tst.lib"; |
---|
2 | tst_init(); |
---|
3 | |
---|
4 | // raducal in alga.extensions |
---|
5 | LIB "primdec.lib"; |
---|
6 | ring R=(0,a),(x,y,z),dp; |
---|
7 | minpoly=a2+a+1; |
---|
8 | |
---|
9 | poly L=x; |
---|
10 | poly u1=1; |
---|
11 | poly u2=1; |
---|
12 | poly v1=y; |
---|
13 | poly v2=z; |
---|
14 | poly g1a=(a+2)/(-3)*(u1*v1*a-u2*v2); |
---|
15 | poly g2a=(a+2)/(-3)*(u1*v1-u2*v2); |
---|
16 | poly g1=L*g1a; |
---|
17 | poly g2=L*g2a; |
---|
18 | poly v3=g1a+(a+1)*g2a; |
---|
19 | poly f1=(u1*u2*L^3+v1*v2*v3)/2; |
---|
20 | poly f2=(-u1*u2*L^3+v1*v2*v3)/2; |
---|
21 | |
---|
22 | poly F=g1^3-f1^2; |
---|
23 | ideal I=jacob(F); |
---|
24 | I=std(I); |
---|
25 | hilb(I); |
---|
26 | |
---|
27 | ideal I1=radical(I,"KL"); |
---|
28 | I1=std(I1); |
---|
29 | hilb(I1); |
---|
30 | |
---|
31 | ideal I2=radical(I); |
---|
32 | I2=std(I2); |
---|
33 | hilb(I2); |
---|
34 | |
---|
35 | ideal I3=radicalEHV(I); |
---|
36 | I3=std(I3); |
---|
37 | hilb(I3); |
---|
38 | |
---|
39 | // the Hilbert functions of the radicals must coincide: |
---|
40 | hilb(I1,1); |
---|
41 | hilb(I2,1); |
---|
42 | hilb(I3,1); |
---|
43 | |
---|
44 | tst_status(1);$ |
---|
Note: See
TracBrowser
for help on using the repository browser.