 LIB "finvar.lib";
ring r= 0, (a,b,c,d,e,f),dp;
matrix A1[6][6] = 0,0,0,1,0,0,0,0,1,0,0,0,0,1,0,0,0,0,1,0,0,0,0,0,0,0,0,0,0,1,0,0,0,0,1,0;
matrix A2[6][6] = 0,1,0,0,0,0,1,0,0,0,0,0,0,0,0,0,0,1,0,0,0,0,1,0,0,0,0,1,0,0,0,0,1,0,0,0;
list L = primary_invariants(A1,A2);
matrix IS = irred_secondary_char0(L[1],L[2],L[3],0);
IS;
==> IS[1,1]=ab+cd+ef
==> IS[1,2]=a2d+ad2+b2e+be2+c2f+cf2
==> IS[1,3]=ac2+b2d+a2e+ce2+d2f+bf2
==> IS[1,4]=b2c+bc2+d2e+de2+a2f+af2
==> IS[1,5]=a2c+bd2+c2e+ae2+b2f+df2
==> IS[1,6]=a2cd+abd2+abe2+b2ef+c2ef+cdf2
==> IS[1,7]=ab2d+ac2d+a2be+cd2f+ce2f+bef2
==> IS[1,8]=a2bd+acd2+ab2e+c2df+be2f+cef2
==> IS[1,9]=a3d+ad3+b3e+be3+c3f+cf3
