Trying about as high of a limit as I can

src/problem.rs

@@ -30,7 +30,7 @@ where
         }
     }
     pub fn solve(&mut self) -> Solution<S, D> {
-        let mut convergency = SimpleConvergency{ eps: 1e-12, max_iter: 1000 };
+        let mut convergency = SimpleConvergency{ eps: 1e-12, max_iter: 10000 };
         let mut times: Vec<f64> = vec![self.ode.t];
         let mut states: Vec<SVector<f64, D>> = vec![self.ode.y];
         let mut dense_coefficients: Vec<Vec<SVector<f64, D>>> = Vec::new();