eom 0.11.0

Configurable ODE/PDE solver
Documentation 
1
2
3
4
5
6
7
8
9
10
11
12

use eom::*;
use ndarray::*;

fn main() {
    let dt = 0.01;
    let eom = ode::Lorenz63::default();
    let mut teo = explicit::RK4::new(eom, dt);
    let mut x: Array1<f64> = arr1(&[1.0, 0.0, 0.0]);
    for _ in 0..100_000_000 {
        teo.iterate(&mut x);
    }
}