the dynamical evolution of the comet within the last 10$^{5}$ yr. JFCs are highly chaotic objects, whose dynamic evolution must be studied in terms of statistics [@levison1994]. With this in mind we analysed the evolution of the original object, i.e., comet 21P, by considering the nominal values of its orbital parameters as they are defined in JPL-HORIZONS (orbital solution JPL K182/3). In the analysis, 200 clones were generated following the covariance matrix of its orbital parameters[^6] . We performed the integrations with the numerical package MERCURY [@chambers1999], using the integration algorithm Bulirsch-Stoer with a time-step of 8 d and we included the Sun, all planets and Pluto in the simulation. In addition, we also included non-gravitational forces. The results of the simulations are displayed in Figure \[dynamic\].

We find that the orbits of all the clones in the simulation are very compact for a period of $\sim$1650 yr. After that period, the orbits started to scatter, which was provoked by a close encounter with Jupiter, at mean distance of 0.1 au. Due to the chaotic nature of JFCs, a comparison of results from different authors who applied different methods (e.g., integration algorithm used, the number of clones, how their clones 