# source:git/emacs/test.sing@8d7503

spielwiese
Last change on this file since 8d7503 was 8d7503, checked in by Jens Schmidt <schmidt@…>, 25 years ago
Initial revision git-svn-id: file:///usr/local/Singular/svn/trunk@2394 2c84dea3-7e68-4137-9b89-c4e89433aadc
• Property mode set to `100644`
File size: 800 bytes
Line
1// define a new ring
2option(prot);
3ring r1 = 32003,(x,y,z),ds;
4r1;
5
6// define the polynomial
7int a,b,c,t=11,5,3,0;
8poly f = x^a+y^b+z^(3*c)+x^(c+2)*y^(c-1)+x^(c-1)*y^(c-1)*z3+
9         x^(c-2)*y^c*(y^2+t*x)^2;
10f;
11
12// calculate its Jacobi ideal ...
13ideal i=jacob(f);
14i;
15
16// ... and the standard basis of the Jacobi ideal
17ideal j=std(i);
18
19// calculate the Milnor number ...
20"The Milnor number of f(11,5,3) for t=0 is", vdim(j);
21
22// ... and the Tjurina number
23j=i+f;
24j=std(j);
25"The Tjurina number of f(11,5,3) for t=0 is", vdim(j);
26
27// repeat everything for t=1
28t=1;
29f=x^a+y^b+z^(3*c)+x^(c+2)*y^(c-1)+x^(c-1)*y^(c-1)*z3
30  +x^(c-2)*y^c*(y^2+t*x)^2;
31ideal i1=jacob(f);
32ideal j1=std(i1);
33"The Milnor number of f(11,5,3) for t=1 is", vdim(j1);
34"The Tjurina number of f(11,5,3) for t=1 is", vdim(std(j1+f));
Note: See TracBrowser for help on using the repository browser.