symbolica 0.2.1

A blazing fast computer algebra system
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

use std::sync::Arc;

use symbolica::{
    domains::{
        factorized_rational_polynomial::FactorizedRationalPolynomial, integer::IntegerRing,
        rational_polynomial::RationalPolynomial,
    },
    parser::Token,
    state::{State, Workspace},
};

fn univariate() {
    let mut state = State::new();
    let workspace: Workspace = Workspace::new();

    let var_names = vec!["x".into(), "y".into()];
    let var_map = Arc::new(
        var_names
            .iter()
            .map(|n| state.get_or_insert_var(n).into())
            .collect(),
    );

    let field = IntegerRing::new();

    let rat: RationalPolynomial<_, u8> = Token::parse("1/((x+1)*(x+2)(x^3+2x+1))")
        .unwrap()
        .to_rational_polynomial(&workspace, &mut state, &field, &field, &var_map, &var_names)
        .unwrap();

    println!("Partial fraction {}:", rat.printer(&state));
    for x in rat.apart(0) {
        println!("\t{}", x.printer(&state));
    }
}

fn multivariate() {
    let mut state = State::new();
    let workspace: Workspace = Workspace::new();

    let var_names = vec!["x".into(), "y".into()];
    let var_map = Arc::new(
        var_names
            .iter()
            .map(|n| state.get_or_insert_var(n).into())
            .collect(),
    );

    let field = IntegerRing::new();

    let rat: FactorizedRationalPolynomial<_, u8> = Token::parse("1/((x+y)*(x^2+x*y+1)(x+1))")
        .unwrap()
        .to_factorized_rational_polynomial(
            &workspace, &mut state, &field, &field, &var_map, &var_names,
        )
        .unwrap();

    println!("Partial fraction {} in x:", rat.printer(&state));
    for x in rat.apart(0) {
        println!("\t{}", x.printer(&state));
    }
}

fn main() {
    univariate();
    multivariate();
}