//deform_l.tst
//long  tests for deform.lib
//---------------------------
LIB "tst.lib";
tst_init();
LIB "deform.lib";
 printlevel =  2;
ring  r0  = 0,(x,y,z),ls;
 ideal i2 = y4-xy,yz-xz,z4-xz;   
 list L=versal(i2);
kill L;
//-----------------------------------------------------------------------------
ring  r1  = 0,(x,y,z,u),dp;
 ideal i5 = x3-y2,xz,xu,yz,yu,z3-u2; 
 list L=versal(i5);
 ideal i6 = z2-xz,zu-xu,u2-yu,yz-xu;
 L=versal(i6); kill L;
//------- cone over rational normal curve ofe degree d -------------------------
 int    d = 4;
ring r2   = 0,(x(1..d+1)),ls;
 ideal i7 = maxideal(1);
 matrix m[2][d] = i7[1..d],i7[2..d+1];
       i7 = minor(m,2);
       i7 = minbase(i7);
 list L=versal(i7); 
 def Px=L[1]; setring Px; Fs; print(Js); print(Rs);
 kill L,Px;
//------- L_d_d: d generic lines in A^d ---------------------------------------
        d = 4;
ring r3=0,(x(1..d)),ls;
 int n,m; ideal i8;
 for (n=1;n<d;n=n+1)
 { for (m=n+1;m<=d;m=m+1) {i8=i8,x(n)*x(m);}
 } 
        i8 = simplify(i8,10);
 list L=versal(i8,3,"myring");
 def Px=L[1]; setring Px; Fs;
 kill L,Px;
//------- two elliptic singularity  (t1=7 and 13,t2=3) ------------------------
ring r6   = 0,(x,y,z,u,v,w),ds;
 ideal i9 = uv-zw,u2-yw,zu-yv,yu-xw,xu+v2-uw,yz-xv,y2+v2-uw,xy+zv-yw,x2+z2-yu;
 list L=versal(i9);
 def Px=L[1]; setring Px; Fs;
 kill L,Px;
ring r7   = 0,(x,y,z,u,v),dp;
 ideal I1 = xy,xz,xu,xv,yz,yu,yv,zu,zv,uv;    
 list L=versal(I1);
 def Px=L[1]; setring Px; Fs;
 kill L,Px;

ring r8   = 0,(w(1..9)),dp;
 ideal I3 = w(8)^2-w(6)*w(9),
            w(7)*w(8)-w(5)*w(9),
            w(5)*w(8)-w(3)*w(9),
            w(4)*w(8)-w(2)*w(9),
            w(7)^2-w(4)*w(9),
            w(6)*w(7)-w(3)*w(9),
            w(5)*w(7)-w(2)*w(9),
            w(4)*w(7)-w(1)*w(9),
            w(3)*w(7)-w(2)*w(8),
            w(2)*w(7)-w(1)*w(8),
            w(5)*w(6)-w(3)*w(8),
            w(4)*w(6)-w(2)*w(8),
            w(5)^2-w(3)*w(7),
            w(4)*w(5)-w(2)*w(7),
            w(3)*w(5)-w(2)*w(6),
            w(2)*w(5)-w(1)*w(6),
            w(4)^2-w(1)*w(7),
            w(3)*w(4)-w(2)*w(5),
            w(2)*w(4)-w(1)*w(5),
            w(2)^2-w(1)*w(3),
            w(2)*w(6)*w(8)-w(3)^2*w(9),
            w(1)*w(6)*w(8)-w(2)*w(3)*w(9),
            w(2)*w(6)^2-w(3)^2*w(8),
            w(1)*w(6)^2-w(2)*w(3)*w(8);
list L= versal(I3);
def Px=L[1]; setring Px; Fs;
kill L,Px;
kill r0,r1,r2,r3,r6,r7,r8,m,n,d;
//////////////////////////////////////////////////////////////////////////////////
ring   Po = 0,(x,y),dp;
ideal  Io = std(x^4+y^3);
matrix Mo;
//============= rk 1 ======================================
module k =[x],[y];
module m(0)=k;
module m(1)=[x,y],[-y2,x3];
module m(2)=[x3,y],[-y2,x];
module m(3)=[x2,y],[-y2,x2];
module m(4)=[y,x,0],[-x2,0,-y],[0,-y,-x];
//============= rk 2 ======================================
module m(5)=[xy,-y2,x2],[-y2,-x3,-xy],[-x3,x2y,y2];
module m(6)=[xy,-yy,xx,0],[-yy,-xxx,-xy,0],[-xx,0,0,-y],[0,x2,y,x];
int i';
int ii=6;
for (i'=4;i'<ii;i'=i'+1)
{
  Mo=m(i');
  list L=mod_versal(Mo,Io);
  def Qx=L[2]; 
  print(Ms);
  print(Ls);
  print(Js); 
  kill L,Qx;
}
list L=mod_versal(m(2),Io,4,"qq","W","dp","oo");
def Px=L[1]; def Qx=L[2]; def So=L[3]; def Ox=L[4];
listvar(ring);
show(Qx);
tst_status(1);$