In this paper we study the dynamics of a massless particle around the L1,2 libration points of the Earth–Moon system in a full Solar System gravitational model. The study is based on the analysis of the quasi-periodic solutions around the two collinear equilibrium points. For the analysis and computation of the quasi-periodic orbits, a new iterative algorithm is introduced which is a combination of a multiple shooting method with a refined Fourier analysis of the orbits computed with the multiple shooting. Using as initial seeds for the algorithm the libration point orbits of Circular Restricted Three Body Problem, determined by Lindstedt-Poincare methods, the procedure is able to refine them in the Solar System force-field model for large time-spans, that cover most of the relevant Sun–Earth–Moon periods.