The identification of mean motion resonances in exoplanetary systems or in the Solar System might be cumbersome when several planets and large number of smaller bodies are to be considered. Based on the geometrical meaning of the resonance variable, an efficient method is introduced and described here, by which mean motion resonances can be easily find without any a priori knowledge on them. The efficiency of this method is clearly demonstrated by using known exoplanets engaged in mean motion resonances, and also some members of different families of asteroids and Kuiper-belt objects being in mean motion resonances with Jupiter and Neptune respectively.

In the spatial circular restricted three-body problem librational motion around the Lagrangian points $L_4$ and $L_5$ exists up to high inclinations of the massless object. We report the results of numerical investigations on the stability of this librational motion for systematic variations in the inclination of the Trojan and the mass-ratio $\mu$ of the two massive bodies. We show that stable motion prevails for inclinations below 60 degrees and mass-ratios $\mu < 0.04$ for typical integration times up to $10^6$ revolution periods. At even higher mass-ratios -- beyond the critical mass-ratio for the planar case -- stable orbits were found to exist for up to $10^7$ periods at moderate inclinations. We extracted the librational frequencies on a grid in the parameter space from the Fourier spectra and traced their variation. Several resonances between the short and long period librational frequency as well as the vertical frequency lie inside the investigated region. The application of the Lyapunov characteristic indicator and spectral number methods also reveals the chaotic regions. This simple model is equally applicable in the Solar system for low mass ratios for Trojans of the planets, as well as to Trojan-type exoplanets in binary star systems at high mass-ratios.

The stability of the Lagrangian point $L_4$ is investigated in the elliptic restricted three-body problem by using Floquet's theory. Stable and unstable domains are determined in the parameter plane of the mass parameter and the eccentricity by computing the characteristic exponents. Frequencies of motion around $L_4$ have been determined both in the stable and unstable domains and fitting functions for the frequencies are derived depending on the mass parameter and the eccentricity. Resonances between the frequencies are studied in the whole parameter plane. It is shown that the 1:1 resonances are not restricted only to single curves but extend to the whole unstable domain. In the unstable domains longer escape times of the test particle from the neighbourhood of $L_4$ are related to certain resonances, but changing the parameters the same resonances may lead to faster escape.