pub fn solar_declination(day_of_year: f64) -> f64
Compute solar declination for a given day of year (1-365) Returns angle in degrees