pub type SystemFunction = fn(time: f64, state: &State) -> DState;
Function that receives as input the time and state of the system and returns the differential of the state.