First time here? Check out the FAQ!
answered 0 years ago
Like this:
mymod = 21 x = polygen(Zmod(mymod)) f = x^4 + x^2 + 1 print( f.roots(multiplicities=False) )