Increase precision in ref. system conversions

[ronisbr]

We can increase the precision in rotations that are about one axis. See, for example:

julia> JD_TT = 2.460115315972222e6
2.460115315972222e6

julia> JD_UTC = jd_tt_to_utc(JD_TT)
2.460115315171481e6

julia> D = r_eci_to_eci(TEME(), JD_UTC, TOD(), JD_UTC)
3×3 StaticArrays.SMatrix{3, 3, Float64, 9} with indices SOneTo(3)×SOneTo(3):
  1.0           3.8393e-5     4.15777e-19
 -3.8393e-5     1.0          -3.91043e-17
 -1.25473e-19  -3.91016e-17   1.0

julia> D[:]
9-element StaticArrays.SVector{9, Float64} with indices SOneTo(9):
  0.9999999992629907
 -3.8392950862527036e-5
 -1.2547321052914035e-19
  3.839295086252501e-5
  0.9999999992629907
 -3.910159364194184e-17
  4.157774012182647e-19
 -3.910427424736158e-17
  0.9999999999999999

Maybe we need a version of angle_to_rot with only one angle argument.