ndarray-odeint 0.7.0

solve ODE using rust-ndarray
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18


extern crate ndarray;
extern crate ndarray_odeint;

use ndarray::rcarr1;
use ndarray_odeint::*;

fn main() {
    let dt = 0.01;
    let eom = model::Roessler::default();
    let teo = explicit::rk4(eom, dt);
    let ts = time_series(rcarr1(&[1.0, 0.0, 0.0]), &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]);
    }
}