CONAN.occultquad_pya.OccultQuadPy.occultquad
============================================

.. py:method:: CONAN.occultquad_pya.OccultQuadPy.occultquad(z0, u1, u2, p, nz=None)

   Occultquad routine

   Python implementation of code provided by Mandel and Agol (2002).
   Please cite Mandel & Agol (2002) if you make use of this routine.
   The original authors are not responsible for the python
   re-implmentation.

   Produces same result as the fortran code but ~30 times slower.

   :param z0: Impact parameter in units of stellar radius
   :type z0: array
   :param u1: Linear limb-darkening coefficient (gamma_1 in paper)
   :type u1: int
   :param u2: Quadratic limb-darkening coefficient (gamma_2 in paper)
   :type u2: int
   :param p: Occulting star (or planet) size in units of stellar radius
   :type p: int
   :param nz: Number of values in z0. Determined here from z0, parameter
              introduced to maintain API.
   :type nz: int, optional

   :returns: * **muo1** (*array*) -- Fraction of flux at each z0 (impact parameter) for a limb-darkened source
             * **mu0** (*array*) -- Fraction of flux at each z0 (impact parameter) for a uniform source