pub fn calculate_new_step_size<F: IntegrateFloat>(
h_current: F,
error_norm: F,
error_order: usize,
safety: F,
) -> Fpub fn calculate_new_step_size<F: IntegrateFloat>(
h_current: F,
error_norm: F,
error_order: usize,
safety: F,
) -> F