The formalism of SUSYQM (SUperSYmmetric Quantum Mechanics) is properly modified in such a way to be suitable for the description and the solution of a classical maximally superintegrable Hamiltonian System, the so-called Taub-Nut system, associated with the Hamiltonian:In full agreement with the results recently derived by A. Ballesteros et al. for the quantum case, we show that the classical Taub-Nut system shares a number of essential features with the Kepler system, that is just its Euclidean version arising in the limit η → 0, and for which a "SUSYQM" approach has been recently introduced by S. Kuru and J. Negro. In particular, for positive η and negative energy the motion is always periodic; it turns out that the period depends upon η and goes to the Euclidean value as η → 0. Moreover, the maximal superintegrability is preserved by the η-deformation, due to the existence of a larger symmetry group related to an η-deformed Runge-Lenz vector, which ensures that in R 3 closed orbits are again ellipses. In this context, a deformed version of the third Kepler's law is also recovered. The closing section is devoted to a discussion of the η < 0 case, where new and partly unexpected features arise.We consider the classical Hamiltonian in R N given by:where k and η are real parameters, q = (q 1 , . . . , q N ), p = (p 1 , . . . , p N ) ∈ R N are conjugate coordinates and momenta, and q 2 ≡ |q| 2 = ∑ N i=1 q 2 i . We recall that H η has been proven to be a maximally superintegrable Hamiltonian by making use of symmetry techniques [1]. This means that H η is endowed with the maximum possible number (2N − 1) of functionally independent constants of motion (including H η itself). In fact, besides the integrals of motion provided by the so(N) symmetry, H η is endowed with an η−deformed ND Laplace-Runge-Lenz vector R implying the existence of N additional constants of motion coming from the components of R, which are given by: