Daylight length follows a sine curve almost perfectly β almost. This explores the extra wobble caused by Earth's tilted, elliptical orbit.
Start here β these are the sources that inspired this exploration.
A plain sinusoidal model gets you close to real daylight-length data, but not exact. The gap comes from two effects layered on top of each other: Earth's axis is tilted, and its orbit is elliptical, so it moves faster in January than in July. This combination, called the equation of time, is why the earliest sunset of the year in the northern hemisphere falls well before the winter solstice, and the latest sunrise falls well after it. Fit a plain sine curve to daylight-length data for one location, then look at where it misses, and see if you can explain the gap using the equation of time rather than just latitude.