pub fn double_pendulum( l1: f64, l2: f64, m1: f64, m2: f64, g: f64, ) -> impl Fn(f64, &[f64]) -> Vec<f64>