diffurch 0.0.3

Numerical solver for ordinary and delay differential equations
Documentation 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32

#![allow(incomplete_features)]
#![feature(generic_const_exprs)]
use std::f64::consts::PI;

use diffurch::*;

fn main() {
    let eq = equation!(|[x, dx]| [dx, -4. * x]);
    let range = 0. ..PI * 3.5;

    let mut t_x = Vec::new();
    let mut t_dx = Vec::new();

    Solver::new()
        .rk(&rk::RK98)
        .stepsize(1.)
        .on_step(
            event!(|t, [x, dx]| [(t as f32, x as f32), (t as f32, dx as f32)])
                .subdivide(10)
                .to_vecs([&mut t_x, &mut t_dx]),
        )
        .run(eq, [0., 1.], range.clone());

    use textplots::*;
    Chart::new(160, 80, range.start as f32, range.end as f32)
        .lineplot(&textplots::Shape::Lines(&t_x))
        .linecolorplot(
            &textplots::Shape::Lines(&t_dx),
            rgb::Rgb { r: 0, g: 0, b: 255 },
        )
        .display();
}