Function rod_moment_of_inertia_about_center 

pub fn rod_moment_of_inertia_about_center(mass: f64, length: f64) -> Option<f64>

Computes rod moment of inertia about its center using I = (1 / 12)mL².