  It has been previously shown that moons of extrasolar planets may be
detectable with the Kepler Mission, for moon masses above ~0.2 Earth masses
Kipping et al. 2009c. Transit timing effects have been formerly identified as a
potent tool to this end, exploiting the dynamics of the system. In this work,
we explore the simulation of transit light curves of a planet plus a single
moon including not only the transit timing effects but also the light curve
signal of the moon itself. We introduce our new algorithm, LUNA, which produces
transit light curves for both bodies, analytically accounting for shadow
overlaps, stellar limb darkening and planet-moon dynamical motion. By building
the dynamics into the core of LUNA, the routine automatically accounts for
transit timing/duration variations and ingress/egress asymmetries for not only
the planet, but also the moon. We then generate some artificial data for two
feasibly detectable hypothetical systems of interest: a i) prograde and ii)
retrograde Earth-like moon around a habitable-zone Neptune for a M-dwarf
system. We fit the hypothetical systems using LUNA and demonstrate the
feasibility of detecting these cases with Kepler photometry.
