Ver 0.37.0

- Whole new root finding
- from thiserror to anyhow
- Connect with Zenodo

README.md

@@ -256,7 +256,7 @@ The example code demonstrates how Peroxide can be used to simulate the Lorenz at
 use peroxide::fuga::*;
 
 fn main() -> Result<(), Box<dyn Error>> {
-    let rkf45 = RKF45::new(1e-4, 0.9, 1e-6, 1e-1, 100);
+    let rkf45 = RKF45::new(1e-4, 0.9, 1e-6, 1e-2, 100);
     let basic_ode_solver = BasicODESolver::new(rkf45);
     let (_, y_vec) = basic_ode_solver.solve(
         &Lorenz,

examples/lorenz_rkf45.rs

@@ -2,7 +2,7 @@ use peroxide::fuga::*;
 
 #[allow(unused_variables)]
 fn main() -> Result<(), Box<dyn Error>> {
-    let rkf45 = RKF45::new(1e-4, 0.9, 1e-6, 1e-1, 100);
+    let rkf45 = RKF45::new(1e-4, 0.9, 1e-6, 1e-2, 100);
     let basic_ode_solver = BasicODESolver::new(rkf45);
     let (_, y_vec) = basic_ode_solver.solve(
         &Lorenz,

src/util/plot.rs

@@ -21,8 +21,9 @@
 //!     let y1 = x.fmap(|t| t.powi(2));
 //!     let y2 = x.fmap(|t| t.powi(3));
 //!
+//!     let mut rng = SmallRng::seed_from_u64(42);
 //!     let normal = Normal(0f64, 0.1);
-//!     let eps = normal.sample(100);
+//!     let eps = normal.sample_with_rng(&mut rng, x.len());
 //!     let y3 = y2.add_v(&eps);
 //!
 //!     let mut plt = Plot2D::new();