To quantify brightness of the stars
in maxima (MaxI in the phase 0.25, and MaxII in the phase 0.75)
and in minima (MinI for a primary minimum in the phase 0.0
and MinII for a secondary minimum in the phase 0.5), we performed
a weighted quadratic polynomial fit of the observed data
around these extremas in the phase interval ±0.1.