Numerical stability analysis around L4 in the post-Newtonian restricted three-body problem

The post-Newtonian equations of motion of the restricted, circular three-body problem are integrated numerically. Using this integration the stable region around the triangular libration point L4 for a relativistic system of mass ratio earth-moon is calculated. The set of stable points is compared with the results of McKenzie and Szebehely found for the Newtonian case. The stable region shows a very complex structure which shows similarities to the one found in the Newtonian problem.