 LIB "finvar.lib";
ring R=3,(x,y,z),dp;
matrix A[3][3]=0,1,0,1,0,0,0,0,1;
list L=primary_invariants(A,intvec(1,1,0));
// In that example, there are no secondary invariants
// in degree 1 or 2.
matrix IS=irred_secondary_no_molien(L[1..2],intvec(1,2),1);
==>
==> Searching irred. sec. inv. in degree 3
==> We have 4 candidates for irred. secondaries
==> We found irr. sec. inv. number 1 in degree 3
==> We found irr. sec. inv. number 2 in degree 3
==> Searching irred. sec. inv. in degree 4
==> We have 1 candidates for irred. secondaries
==> We found irr. sec. inv. number 3 in degree 4
==> Searching irred. sec. inv. in degree 5
==> Searching irred. sec. inv. in degree 6
==> Searching irred. sec. inv. in degree 7
==> Searching irred. sec. inv. in degree 8
==> Searching irred. sec. inv. in degree 9
==> Searching irred. sec. inv. in degree 10
==> Searching irred. sec. inv. in degree 11
==> Searching irred. sec. inv. in degree 12
==> Searching irred. sec. inv. in degree 13
print(IS);
==> x2zy2z,xyz,x3yxy3
