astrodynamics 0.1.0

Numerical astrodynamics engine for orbit propagation, force models, and flight-dynamics primitives
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

pub struct DP54Tableau {
    pub c: [f64; 7],
    pub a: Vec<Vec<f64>>,
    pub b5: [f64; 7],
    pub b4: [f64; 7],
}

impl Default for DP54Tableau {
    fn default() -> Self {
        // Dormand-Prince 5(4) coefficients (DOPRI5)
        Self {
            c: [0.0, 1.0/5.0, 3.0/10.0, 4.0/5.0, 8.0/9.0, 1.0, 1.0],
            a: vec![
                vec![],
                vec![1.0/5.0],
                vec![3.0/40.0, 9.0/40.0],
                vec![44.0/45.0, -56.0/15.0, 32.0/9.0],
                vec![19372.0/6561.0, -25360.0/2187.0, 64448.0/6561.0, -212.0/729.0],
                vec![9017.0/3168.0, -355.0/33.0, 46732.0/5247.0, 49.0/176.0, -5103.0/18656.0],
                vec![35.0/384.0, 0.0, 500.0/1113.0, 125.0/192.0, -2187.0/6784.0, 11.0/84.0],
            ],
            b5: [35.0/384.0, 0.0, 500.0/1113.0, 125.0/192.0, -2187.0/6784.0, 11.0/84.0, 0.0],
            b4: [5179.0/57600.0, 0.0, 7571.0/16695.0, 393.0/640.0, -92097.0/339200.0, 187.0/2100.0, 1.0/40.0],
        }
    }
}