subgrouplist(bnrinit(bnfinit(x^3-x-1),10));
bnrL1(bnrinit(bnfinit(x^2+1),10,1),1);
setrand(5);bnrL1(bnrinit(bnfinit(y^2+6),1,1),0);
