1 | LIB "deform.lib"; |
---|
2 | ring R=32003,(x,y,z),ds; |
---|
3 | //---------------------------------------------------- |
---|
4 | // hypersurface case (from series T[p,q,r]): |
---|
5 | int p,q,r = 3,3,4; |
---|
6 | poly f = x^p+y^q+z^r+xyz; |
---|
7 | print(deform(f)); |
---|
8 | // the miniversal deformation of f=0 is the projection from the |
---|
9 | // miniversal total space to the miniversal base space: |
---|
10 | // { (A,B,C,D,E,F,G,H,x,y,z) | x3+y3+xyz+z4+A+Bx+Cxz+Dy+Eyz+Fz+Gz2+Hz3 =0 } |
---|
11 | // --> { (A,B,C,D,E,F,G,H) } |
---|
12 | //---------------------------------------------------- |
---|
13 | // complete intersection case (from series P[k,l]): |
---|
14 | int k,l =3,2; |
---|
15 | ideal j=xy,x^k+y^l+z2; |
---|
16 | print(deform(j)); |
---|
17 | list L=versal(j); // using default names |
---|
18 | def Px=L[1]; |
---|
19 | setring Px; |
---|
20 | show(Px); // show is a procedure from inout.lib |
---|
21 | listvar(ideal); |
---|
22 | // ___ Equations of miniversal base space ___: |
---|
23 | Js; |
---|
24 | // ___ Equations of miniversal total space ___: |
---|
25 | Fs; |
---|
26 | kill L; |
---|
27 | // the miniversal deformation of V(j) is the projection from the |
---|
28 | // miniversal total space to the miniversal base space: |
---|
29 | // { (A,B,C,D,E,F,x,y,z) | xy+A+Bz=0, y2+z2+x3+C+Dx+Ex2+Fy=0 } |
---|
30 | // --> { (A,B,C,D,E,F) } |
---|
31 | //---------------------------------------------------- |
---|
32 | // general case (cone over rational normal curve of degree 4): |
---|
33 | ring r1=0,(x,y,z,u,v),ds; |
---|
34 | matrix m[2][4]=x,y,z,u,y,z,u,v; |
---|
35 | ideal i=minor(m,2); // 2x2 minors of matrix m |
---|
36 | // Def_r will be the name of the miniversal base space with |
---|
37 | // parameters A(1),...,A(4) |
---|
38 | versal(i,0,"Def_r","A("); |
---|
39 | // the miniversal deformation of V(i) is the projection from the |
---|
40 | // miniversal total space to the miniversal base space: |
---|
41 | // { (A(1..4),x,y,z,u,v) | |
---|
42 | // -y^2+x*z+A(2)*x-A(3)*y=0, -y*z+x*u-A(1)*x-A(3)*z=0, |
---|
43 | // -y*u+x*v-A(3)*u-A(4)*z=0, -z^2+y*u-A(1)*y-A(2)*z=0, |
---|
44 | // -z*u+y*v-A(2)*u-A(4)*u=0, -u^2+z*v+A(1)*u-A(4)*v=0 } |
---|
45 | // --> { A(1..4) | |
---|
46 | // -A(1)*A(4) = A(3)*A(4) = -A(2)*A(4)-A(4)^2 = 0 } |
---|
47 | //---------------------------------------------------- |
---|
48 | LIB "tst.lib";tst_status(1);$ |
---|