Orbital Propagation Around Irregular Celestial Bodies Using The Harmonic Balance Method