pub fn strong_error(exact: &[f64], numerical: &[f64]) -> f64
Strong error: max |exact(t_i) - numerical(t_i)|.