Function critical_damping_coefficient 

pub fn critical_damping_coefficient(
    mass: f64,
    spring_constant: f64,
) -> Option<f64>

Computes the critical damping coefficient using c_critical = 2 * sqrt(mk).