Function sun_position 

pub fn sun_position(at: OffsetDateTime) -> Equatorial

Apparent geocentric solar position (equatorial, epoch of date). Feeds twilight and crate::moon_phase_angle internally.

Meeus ch. 25 low-accuracy method: ~0.01° (36″), inside the crate’s ≤1′ contract (pinned against AstroPy get_sun).

use skymath::sun_position;
use time::macros::datetime;

// Meeus example 25.a.
let sun = sun_position(datetime!(1992-10-13 00:00 UTC));
assert!((sun.ra().degrees() - 198.38083).abs() < 5e-3);
assert!((sun.dec().degrees() + 7.78507).abs() < 5e-3);