Function rk4_step_fixed 

pub fn rk4_step_fixed<const N: usize, F>(
    state: &mut [f32; N],
    dt: f32,
    t: f32,
    derivative: F,
)
where
    F: Fn(&[f32; N], f32) -> [f32; N],