pub fn effective_radius(r_a: f32, r_b: f32) -> f32
Effective radius for two spheres in contact 1/R_eff = 1/R_a + 1/R_b