Changeset c91fefd in git for Tst/Short/equising_s.tst


Ignore:
Timestamp:
Nov 8, 2005, 7:17:25 PM (18 years ago)
Author:
Hans Schönemann <hannes@…>
Branches:
(u'spielwiese', 'fe61d9c35bf7c61f2b6cbf1b56e25e2f08d536cc')
Children:
120c517824adf0e2f616cd70afef52a65d2683a7
Parents:
09171bdbde9e017eae7fe839e8512982e2f9d70f
Message:
*lossen: new lib, more examples


git-svn-id: file:///usr/local/Singular/svn/trunk@8794 2c84dea3-7e68-4137-9b89-c4e89433aadc
File:
1 edited

Legend:

Unmodified
Added
Removed
  • Tst/Short/equising_s.tst

    r09171bd rc91fefd  
    135135reduce(F,std(M[1][1]));
    136136kill LL,Px,r;
     137
    137138/////////////////////////////////////////////////
    138 
     139//              MEGA 05 Example
     140ring R = 0, (x,y), ls;
     141poly f = (y4-x4)^2 - x10;
     142list L = versal(f);        //compute semiuniversal deformation
     143def Px=L[1]; setring Px;
     144poly F = Fs[1,1];         
     145list M=esStratum(F);       //compute the equisingularity stratum
     146def ESSring = M[1]; setring ESSring;
     147option(redSB);
     148ES = std(ES);
     149size(ES);    //-> 42
     150ES[1];       //-> 8*A(1)+8*A(22)+A(1)^3
     151ES[34];      //-> 8*A(40)-A(1)^2+A(1)*A(22)
     152poly F = reduce(imap(Px,F),ES);  //A(1),A(22) both appear in F
     153poly g = subst(F, A(22), -A(1)-(1/8)*A(1)^3);
     154for (int ii=1; ii<=44; ii++){ g = subst(g,A(ii),random(1,100)); }
     155setring R;
     156milnor(f);                //-> 57
     157milnor(imap(ESSring,g));  //-> 57
     158kill R,L,Px,M,ESSring,ii;
    139159
    140160
Note: See TracChangeset for help on using the changeset viewer.