lorenz-cli 0.1.1

Plots lorenz attractor in terminal.
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
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49

//! # lorenz-cli
//!
//! Plots lorenz attractor in terminal.
//! 
//! # Exit
//! to terminate the plotter enter `q` or `Esc`

use std::{error::Error, process};

const SIGMA: f64 = 10.0;
const RHO: f64 = 28.0;
const BETA: f64 = 8.0 / 3.0;

pub fn generate_lorenz_points() -> Vec<(f64, f64)> {
    let dt = 0.01;
    let num_steps = 10000;

    let mut x = 5.0;
    let mut y = 18.0;
    let mut z = 10.0;
    let mut points = Vec::with_capacity(num_steps);

    for _ in 0..num_steps {
        let (dx, dy, dz) = lorenz(x, y, z, SIGMA, RHO, BETA).unwrap_or_else(|err| {
            eprintln!("{:?}", err);
            process::exit(1);
        });
        x += dx * dt;
        y += dy * dt;
        z += dz * dt;
        points.push((x, y));
    }
    points
}

pub fn lorenz(
    x: f64,
    y: f64,
    z: f64,
    sigma: f64,
    rho: f64,
    beta: f64,
) -> Result<(f64, f64, f64), Box<dyn Error>> {
    let dx = sigma * (y - x);
    let dy = x * (rho - z) - y;
    let dz = x * y - beta * z;

    Ok((dx, dy, dz))
}