of interest may be the cyclotomic (http://hackage.haskell.org/packages/archive/cyclotomic/0.3.1...) package for haskell, which can perform exact computation on a decent amount of values.