Qi = nfinit(x^2+1);
idealintersect(Qi,0,1)
Q = nfinit(y);
idealaddtoone(Q,[1,[;]])
idealaddtoone(Q,[1,[;],Mat(1/2)])
ideallog(Q,2,idealstar(Q,4,1))
idealstar(Qi,matdiagonal([6,12]))
idealintersect(Qi, 1/2,1/2)

