Analytic treatment of complete and incomplete geodesics in Taub-NUT space-times