russell_lab 1.15.0

Scientific laboratory for linear algebra and numerical mathematics
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
50
51
52
53
54

use plotpy::{Curve, Plot};
use russell_lab::math::{GOLDEN_RATIO, PI};
use russell_lab::*;

const OUT_DIR: &str = "/tmp/russell_lab/";

fn main() -> Result<(), StrError> {
    // F
    let xx = Vector::linspace(0.0, PI / 2.0, 101)?;
    let yy = xx.get_mapped(|x| math::elliptic_f(x, 0.5).unwrap());
    let mut curve = Curve::new();
    curve
        .set_line_color("#E9708E")
        .set_line_width(2.5)
        .draw(xx.as_data(), yy.as_data());
    let path = format!("{}/math_elliptic_functions_f.svg", OUT_DIR);
    let mut plot = Plot::new();
    plot.add(&curve)
        .grid_labels_legend("$\\phi$", "$F(\\phi,1/2)$")
        .set_figure_size_points(GOLDEN_RATIO * 280.0, 280.0)
        .save(&path)?;

    // E
    let xx = Vector::linspace(0.0, PI / 2.0, 101)?;
    let yy = xx.get_mapped(|x| math::elliptic_e(x, 0.5).unwrap());
    let mut curve = Curve::new();
    curve
        .set_line_color("#4C689C")
        .set_line_width(2.5)
        .draw(xx.as_data(), yy.as_data());
    let path = format!("{}/math_elliptic_functions_e.svg", OUT_DIR);
    let mut plot = Plot::new();
    plot.add(&curve)
        .grid_labels_legend("$\\phi$", "$E(\\phi,1/2)$")
        .set_figure_size_points(GOLDEN_RATIO * 280.0, 280.0)
        .save(&path)?;

    // Π
    let xx = Vector::linspace(0.0, 0.9 * PI / 2.0, 101)?;
    let yy = xx.get_mapped(|x| math::elliptic_pi(1.0, x, 0.5).unwrap());
    let mut curve = Curve::new();
    curve
        .set_line_color("#58B090")
        .set_line_width(2.5)
        .draw(xx.as_data(), yy.as_data());
    let path = format!("{}/math_elliptic_functions_pi.svg", OUT_DIR);
    let mut plot = Plot::new();
    plot.add(&curve)
        .grid_labels_legend("$\\phi$", "$\\Pi(1,\\phi,1/2)$")
        .set_figure_size_points(GOLDEN_RATIO * 280.0, 280.0)
        .save(&path)?;

    Ok(())
}