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
- Residuals
- Work in Progress:
-- Orbital Period and Elements
-- Past and future times of apoapsis and periapsis
-- Prediction of position angle and angular separation at a given time

Github

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

Outputs
- Light Curve and Orbital Diagram
- Eclipse Luminosity and Duration
- Points of Contact
- Impact Parameter
- L1 Point (Roche Lobe Overflow)

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

Desmos Github