pub fn calc_curvature(control_points: &[(f64, f64)], t: f64) -> f64
Calculate curvature at parameter t