OK, I have managed to do Gaussian fits. First, I define a delta phi search region, then identify the faintest magnitude within the search region, smooth out the noise, determine eclipse well and then fit a Gaussian curve to the well. Presenting the first O-C table. Looks like it’s linear. ANOVA period determined incorrectly?