Function cap_solid_angle 

pub fn cap_solid_angle(alpha: f64) -> f64

Solid angle of a spherical cap with half-angle alpha (radians).

Formula: Ω = 2π(1 − cos α). Returns a value in [0, 4π].