pub fn boys_function(n: usize, t: f64) -> f64
Boys function $F_n(t) = \int_0^1 u^{2n} e^{-tu^2} du$.
Uses series expansion for small t and asymptotic form for large t.