“…The geodesic of the generalized Taub-NUT metrics properly explains the motion of well-separated monopole-monopole interactions which provide non-trivial generalizations of the Kepler and harmonic oscillator systems [3,4,5,6,7,8,9,10,11]. The 5D Kepler systems interacting with su(2) monopoles or Yang-Coulomb monopoles have been investigated in [12,13,14,15,16,17,18,19,20,21,22]. Recently, there is a lot of international research interest in generalizing such monopole models to higher dimensions via different approaches [11,23,24,25,26,27].…”