 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);
matrix S,IS=secondary_charp(L[1..size(L)],1);
==>
==> We need to find
==> 1 secondary invariant in degree 0
==> 0 secondary invariants in degree 1
==> 0 secondary invariants in degree 2
==> 2 secondary invariants in degree 3
==> 1 secondary invariant in degree 4
==> In degree 0 we have: 1
==>
==> Searching in degree 3 , we need to find 2 invariant(s)...
==> Looking for Power Products...
==> There are 2 irreducible secondary invariants in degree 3
==> We found all 2 irreducibles in degree 3
==>
==> Searching in degree 4 , we need to find 1 invariant(s)...
==> Looking for Power Products...
==> There are 1 irreducible secondary invariants in degree 4
==> We found all 1 irreducibles in degree 4
==>
==> We're done!
==>
print(S);
==> 1,xyz,x2zy2z,x3yxy3
print(IS);
==> xyz,x2zy2z,x3yxy3
