 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 S,IS=secondary_and_irreducibles_no_molien(L[1..2],intvec(1,2),1);
==>
==> We need to find 4 secondary invariants.
==>
==> In degree 0 we have: 1
==>
==> Searching in degree 3
==> Looking for Power Products...
==> Looking for irreducible secondary invariants in degree 3
==> We found irreducible sec. inv. number 1 in degree 3
==> We found irreducible sec. inv. number 2 in degree 3
==>
==> Searching in degree 4
==> Looking for Power Products...
==> Looking for irreducible secondary invariants in degree 4
==> We found irreducible sec. inv. number 1 in degree 4
==>
==>
==> We're done!
==>
print(S);
==> 1,xyz,x2zy2z,x3yxy3
print(IS);
==> xyz,x2zy2z,x3yxy3
