Function hertz_contact_force 

pub fn hertz_contact_force(
    eff_young: f32,
    eff_radius: f32,
    penetration: f32,
) -> f32

Hertz contact force: F = (4/3) * E_eff * sqrt(R_eff) * delta^(3/2)