Takes observation data inputs of angular separation and position angle and plots a fitted orbit using a non-linear least squares method with the primary star at the origin.
Uses Matplotlib, NumPy, and SciPy

Calculates:
- Sky Projected Orbital Elements (a, e, ω)
- Apoapsis and periapsis positions and distances
- Prediction of position angle and angular separation at a given time (WIP)
- Past and future times of apoapsis and periapsis (WIP)
- Residuals

Github

Simulates a light curve of eclipsing binary stars
Uses Matplotlib and NumPy

Input Variable Units
- Stellar Masses (m1, m2): Solar Masses
- Stellar Radii (r1, r2): Solar Radii
- Luminosities (L1, L2): Solar Luminosities
- Orbital Period: Days

Assumptions
- Uniform Luminosity Distribution (No limb darkening)
- Circular Orbit (e = 0)
- Celestial bodies m1 and m2 is perfectly spherical with projected areas being perfect circular

Desmos Github