θ is negated from example (no idea why), for 2024 eclipse, returns completely wrong result.
Instead of using this, compute the time of greatest eclipse, then use the time to find the point on the curve of max eclipse at that time.