default(realprecision,38);
setrand(1); K=bnfinit(x^2+31);bnrL1(bnrinit(K,1,1))
setrand(1); K=bnfinit(x^2-3);
setrand(1); bnrL1(bnrinit(K,5,1),,5)
setrand(1); bnrL1(bnrinit(K,[5,[1,0]],1),,5)
setrand(1); bnrL1(bnrinit(K,[5,[1,1]],1),,5)
setrand(1); K=bnfinit(x^2-168);\\#1601
setrand(1); bnrL1(bnrinit(K,[6,[1,1]],1),,5)
