I'm teaching a course on dynamical systems based on the book by Devaney. Today I tried to show the map f(θ)=2θ on the circle using Sage, where θ is given in radians. I could not find functions already defined to perform addition in R modulo 2π, so I tried to devise my own. I stumbled upon this problem:
atan(sin(1/5pi)/cos(1/5pi))
produces
arctan(4sin(1/5pi)/(sqrt(5) + 1)).
Is there a way to obtain 1/5*pi?
(By the way, _.simplify() does nothing in this case.)