The escape dynamics around the triangular Lagrangian point L5 in the real Sun-Earth-Moon-Spacecraft system is investigated. Appearance of the finite time chaotic behaviour suggests that widely used methods and concepts of dynamical system theory can be useful in constructing a desired mission design. Existing chaos control methods are modified in such a way that we are able to protect a test particle from escape. We introduce initial condition maps in order to have a suitable numerical method to describe the motion in high dimensional phase space. Results show that the structure of initial condition maps can be split into two well-defined domains. One of these two parts has a regular contiguous shape and is responsible for long time escape; it is a long-lived island. The other one shows a filamentary fractal structure in initial condition maps. The short time escape is governed by this object. This study focuses on a low-cost method which successfully transfers a reference trajectory between these two regions using an appropriate continuous control force. A comparison of the Earth-Moon transfer is also presented to show the efficiency of our method.

The size distribution of the stability region around the Lagrangian point L4 is investigated in the elliptic restricted three-body problem as the function of the mass parameter and the orbital eccentricity of the primaries. It is shown that there are minimum zones in the size distribution of the stability regions, and these zones are connected with secondary resonances between the frequencies of librational motions around L4. The results can be applied to hypothetical Trojan planets for predicting values of the mass parameter and the eccentricity for which such objects can be expected or their existence is less probable.

Several integration schemes exits to solve the equations of motion of the $N$-body problem. The Lie-integration method is based on the idea to solve ordinary differential equations with Lie-series. In the 1980s this method was applied for the $N$-body problem by giving the recurrence formula for the calculation of the Lie-terms. The aim of this works is to present the recurrence formulae for the linearized equations of motion of $N$-body systems. We prove a lemma which greatly simplifies the derivation of the recurrence formulae for the linearized equations if the recurrence formulae for the equations of motions are known. The Lie-integrator is compared with other well-known methods. The optimal step size and order of the Lie-integrator are calculated. It is shown that a fine-tuned Lie-integrator can be 30%-40% faster than other integration methods.

The dynamical structure of the phase space of the Pluto--Charon system is studied in the model of the spatial circular restricted three-body problem by using numerical methods. With the newly discovered two small satellites S/2005 P1 and S/2005 P2, the Pluto--Charon system can be considered as the first known binary system in which celestial bodies move in P-type orbits. It is shown that the two satellites are in the stable region of the phase space and their origin by capture is unlikely. Also the large mass parameter allows to regard the satellites as a model of a new class of exoplanets orbiting around stellar binary systems.