We present the complete set of analytical solutions of the geodesic equation
in Taub-NUT space-times in terms of the Weierstrass elliptic function. We
systematically study the underlying polynomials and characterize the motion of
test particles by its zeros. Since the presence of the "Misner string" in the
Taub-NUT metric has led to different interpretations, we consider these in
terms of the geodesics of the space-time. In particular, we address the
geodesic incompleteness at the horizons discussed by Misner and Taub, and the
analytic extension of Miller, Kruskal and Godfrey, and compare with the
Reissner-Nordstr\"om space-time.Comment: 22 pages, 14 figures, accepted for publication in PR