pub type ODEFunc = Box<dyn Fn(f64, &State) -> State>;
A function that computes the derivative du/dt given time t and state u.
pub struct ODEFunc(/* private fields */);