scirs2-special 0.4.2

Special functions module for SciRS2 (scirs2-special)
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
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95

use scirs2_special::*;

#[allow(dead_code)]
fn main() {
    // Test j0 values
    println!("J0 function values:");
    println!("j0(0.5) = {}", j0(0.5));
    println!("j0(1.0) = {}", j0(1.0));
    println!("j0(5.0) = {}", j0(5.0));
    println!("j0(10.0) = {}", j0(10.0));
    println!();

    // Test j1 values
    println!("J1 function values:");
    println!("j1(0.5) = {}", j1(0.5));
    println!("j1(1.0) = {}", j1(1.0));
    println!("j1(5.0) = {}", j1(5.0));
    println!("j1(10.0) = {}", j1(10.0));
    println!();

    // Test i0 values
    println!("I0 function values:");
    println!("i0(0.5) = {}", i0(0.5));
    println!("i0(1.0) = {}", i0(1.0));
    println!("i0(5.0) = {}", i0(5.0));
    println!();

    // Test i1 values
    println!("I1 function values:");
    println!("i1(0.5) = {}", i1(0.5));
    println!("i1(1.0) = {}", i1(1.0));
    println!("i1(5.0) = {}", i1(5.0));
    println!();

    // Test y0 values
    println!("Y0 function values:");
    println!("y0(0.5) = {}", y0(0.5));
    println!("y0(1.0) = {}", y0(1.0));
    println!("y0(5.0) = {}", y0(5.0));
    println!("y0(10.0) = {}", y0(10.0));
    println!();

    // Test jn values for n=2
    println!("Jn function values for n=2:");
    println!("jn(2, 0.5) = {}", jn(2, 0.5));
    println!("jn(2, 1.0) = {}", jn(2, 1.0));
    println!("jn(2, 5.0) = {}", jn(2, 5.0));
    println!();

    // Test iv values for v=2
    println!("Iv function values for v=2:");
    println!("iv(2.0, 0.5) = {}", iv(2.0, 0.5));
    println!("iv(2.0, 1.0) = {}", iv(2.0, 1.0));
    println!("iv(2.0, 5.0) = {}", iv(2.0, 5.0));
    println!();

    // Test gamma function values
    println!("Gamma function values:");
    println!("gamma(0.5) = {}", gamma(0.5));
    println!("gamma(1.0) = {}", gamma(1.0));
    println!("gamma(5.0) = {}", gamma(5.0));
    println!();

    // Test gammaln function values
    println!("Gammaln function values:");
    println!("gammaln(0.5) = {}", gammaln(0.5));
    println!("gammaln(1.0) = {}", gammaln(1.0));
    println!("gammaln(5.0) = {}", gammaln(5.0));
    println!();

    // Test beta function values
    println!("Beta function values:");
    println!("beta(0.5, 0.5) = {}", beta(0.5, 0.5));
    println!("beta(1.0, 1.0) = {}", beta(1.0, 1.0));
    println!("beta(5.0, 2.0) = {}", beta(5.0, 2.0));
    println!();

    // Test betaln function values
    println!("Betaln function values:");
    println!("betaln(0.5, 0.5) = {}", betaln(0.5, 0.5));
    println!("betaln(1.0, 1.0) = {}", betaln(1.0, 1.0));
    println!("betaln(5.0, 2.0) = {}", betaln(5.0, 2.0));
    println!();

    // Test Wright omega values
    println!("Wright omega values:");
    println!(
        "wright_omega_real(0.0, 1e-10) = {}",
        wright_omega_real(0.0, 1e-10).expect("Operation failed")
    );
    println!(
        "wright_omega_real(1.0, 1e-10) = {}",
        wright_omega_real(1.0, 1e-10).expect("Operation failed")
    );
}