jivanu 1.0.0

Jivanu — microbiology engine for growth kinetics, metabolism, genetics, and epidemiology
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

//! Basic examples of jivanu microbiology formulas.

fn main() {
    // Growth: Monod kinetics
    let mu = jivanu::growth::monod_kinetics(5.0, 0.5, 1.0).unwrap();
    println!("Monod growth rate (S=5, mu_max=0.5, K_s=1): {mu:.3} /hr");

    // Growth: doubling time
    let td = jivanu::growth::doubling_time(mu).unwrap();
    println!("Doubling time: {td:.2} hr");

    // Metabolism: Michaelis-Menten
    let v = jivanu::metabolism::michaelis_menten(2.0, 10.0, 1.0).unwrap();
    println!("Michaelis-Menten rate (S=2, Vmax=10, Km=1): {v:.2}");

    // Genetics: GC content
    let gc = jivanu::genetics::gc_content("ATGCATGCATGC").unwrap();
    println!("GC content of ATGCATGCATGC: {:.1}%", gc * 100.0);

    // Epidemiology: R0 and herd immunity
    let r0 = jivanu::epidemiology::r_naught(0.5, 0.2).unwrap();
    let h = jivanu::epidemiology::herd_immunity_threshold(r0).unwrap();
    println!("R0 = {r0:.1}, herd immunity = {:.1}%", h * 100.0);

    // Hardy-Weinberg
    let (p2, pq2, q2) = jivanu::genetics::hardy_weinberg(0.6).unwrap();
    println!("Hardy-Weinberg (p=0.6): AA={p2:.2}, Aa={pq2:.2}, aa={q2:.2}");
}