1use num::Complex;
2use poly::{coefficients, Polynomial};
3
4fn main() {
5 let a = Polynomial::new(coefficients![1f32, 2.0, 3.0, 0.0]);
6 let b = Polynomial::new(coefficients![1f32, 0.0, 1.0]);
7 let (q, r) = a.div_rem(&b);
8 println!(
9 "({0}) / ({1}) = ({1}) * ({2}) + {3}",
10 a.to_display("ω"),
11 b.to_display("ω"),
12 q.to_display("ω"),
13 r.to_display("ω")
14 );
15
16 let x = Complex::new(0f32, 1.0);
17 let e = a.eval(x);
18 println!("{} = {} for x = {}", a.to_display("x"), e, x);
19
20 let d = a.eval_der(1f32, 2);
21 println!("({})'' = {} for z = {}", a.to_display("z"), d, 1f32);
22}