pub fn deflection_angle(mass: f64, impact_param: f64, c: f64, G: f64) -> f64
Deflection angle for a light ray passing a mass M at impact parameter b. alpha = 4GM / (b * c^2)