I want to expand cos(2πn) for n=5, 17, 257,65537 to radicals, for example: cos(2π5)=14√5−14
Is this a possible task?
x = var('x')
test = (x+1)^2 # test
test.expand() # works :)
term = cos(2*pi/5)
term.expand() # works :)
term = cos(2*pi/17)
term.expand() # works not :(
PS: I tried
cos(2*pi/17).trig_simplify()
cos(2*pi/17).trig_expand()
as well.