Patched conics give you reasonable starting values. If you have accurate models of Sun, Earth, and Moon motion in your solution, integrating the state vector at 0.1 second intervals will work very nicely. My celestial body models are harmonic equations with about 200 terms each, but that's because I use the same models for interplanetary orbits.