? idealpow(nfinit(Pol([1,6,6^3,6])),x^6+1,7,1)==matid(3)

0

All earlier version for many years until very recently would return 1.  

Thanks,

Igor