ndarray-odeint

solve ODE with rust-ndarray

Example

use itertools::iterate;
use ndarray_odeint::lorenz63 as l63;

let dt = 0.01;
let p = l63::Parameter::default();
let l = |y| l63::f(p, y);
let ts = iterate(arr1(&[1.0, 0.0, 0.0]),
                 |y| ndarray_odeint::explicit::rk4(&l, dt, y.clone()));
let end_time = 10000;
for v in ts.take(end_time) {
  println!("{:?}", &v);
}

You can find complete code at src/bin/main.rs

This figure is plotted by matplotlib (see plot script).

Developement status

See issues

• Explicit Schemes
• Implicit Schemes

License

MIT-License, see LICENSE file.