pub fn boys_function(n: u32, x: f64) -> f64
Boys function F_n(x) for nuclear attraction integrals.
F_n(x) = ∫₀¹ t^{2n} exp(-x·t²) dt