I want to simulate the motion of moon around the earth accounting the solar perturbation. I read that Newton theorem for revolving orbits couldn't account for the apsidal precession of moon. Then I read about Hill's solution for restricted three body problem which gives quite good approximation for lunar precession. I found equations for circular restricted three body problem

http://ccar.colorado.edu/imd/2015/documents/CRTBP_Handout.pdf

but I couldn't find for elliptical restricted three body problem. I have been searching for the solution for many weeks but I am in a mess right now. Could anyone explain me about Hill's equation for elliptical restricted three body problem and suggest me how do I proceed with this problem?