code: def eclipseprob(self):
          
          #TODO: incorporate eccentricity/omega for exact calculation?
          s = self.stars
          return ((s['radius_1'] + s['radius_2'])*RSUN /
                  (semimajor(s['P'],s['mass_1'] + s['mass_2'])*AU))
docs: Array of eclipse probabilities.