2012MNRAS.420.1630P

In this paper we describe an algorithm and deduce the related mathematical formulae that allow the computation of observed fluxes in stellar and planetary systems with an arbitrary number of bodies that are part of a transit or occultation event. The method presented does not have any limits or constraints on the geometry and can be applied for almost all of the available limb-darkening models as well. As a demonstration, we apply this scheme to gather information about the orbital inclinations from multiple transiting planetary systems in cases in which mutual transits occur. We also show here that these mutual events constrain the inclinations unambiguously, yielding a complete picture for the whole system.