how to find an ideal which all variables' power is one

i find a function which can do for characteristic 0
but ideal is not my expected

LIB "finvar.lib";
ring R=0,(x,y,z),dp;
matrix A[3][3]=0,1,0,0,0,1,1,0,0;
matrix REY,M=reynolds_molien(A);
matrix P=primary_char0(REY,M);
print(P);

ring R=0,(x,y,z),dp;
matrix A[3][3]=0,1,0,0,0,1,1,0,0;
matrix P=primary_charp_without_random(A,1);
print(P);

ring R=(x,y,z),dp;
matrix A[3][3]=0,1,0,0,0,1,1,0,0;
matrix P=primary_charp_without_random(A,1);
print(P);

LIB "finvar.lib";
ring R=1,(x,y,z),dp;
matrix A[3][3]=0,1,0,0,0,1,1,0,0;
list L=group_reynolds(A);
list l=primary_charp_no_molien_random(L[1],1);
print(l[1]);