Function rod_moment_of_inertia_about_end 

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

Computes rod moment of inertia about one end using I = (1 / 3)mL².