dynamics that the size of the regular tori in phase space decays exponentially with the degree of freedom of the system (Lichtenberg & Liebermann 1983; Altmann & Kantz 2007).

Our primary goal is to avoid the escape and to keep the test particle in a finite volume of the configuration space with minimal energy consumption. To do this we use the basic idea of chaos control where the main argument is that small changes along a chaotic orbit should be efficient enough to reach a desired performance of the system according to some criteria. The focus is on the expression “*small changes*” that can be applied with success, since the chaos is naturally sensitive to small changes and therefore tiny perturbations are suitable for control and manipulation of the chaotic process (Boccaletti & Pucacco 2000). In this study, a slightly modified chaos control method is used which is suitable for our purposes. We construct the stable manifold of the chaotic saddle numerically for a given state of the phase space, called initial condition map (ICM) (Utku, Hagen & Palmer 2015) and then try to adjust the test particle’s velocity (and consequently the position) to be as close as possible to 