Function bohr_orbit_radius 

pub fn bohr_orbit_radius(n: u32) -> Option<f64>

Computes the hydrogen-like Bohr orbit radius for Z = 1 using r_n = a0 * n^2.

Returns None when n == 0 or when the computed result is not finite.