eom 0.11.0

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

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

fn main() {
    let dt = 0.01;
    let eom = ode::Roessler::default();
    let mut teo = explicit::RK4::new(eom, dt);
    let ts = adaptor::time_series(arr1(&[1.0, 0.0, 0.0]), &mut teo);
    let end_time = 50000;
    println!("time,x,y,z");
    for (t, v) in ts.take(end_time).enumerate() {
        println!("{},{},{},{}", t as f64 * dt, v[0], v[1], v[2]);
    }
}