use super::RK;
pub struct Fehlberg45 {}
impl RK for Fehlberg45 {
fn order() -> u8 {
5
}
fn stages() -> usize {
6
}
fn a_coeffs() -> &'static [f64] {
&[
1.0 / 4.0,
3.0 / 32.0,
9.0 / 32.0,
1932.0 / 2197.0,
-7200.0 / 2197.0,
7296.0 / 2197.0,
439.0 / 216.0,
-8.0,
3680.0 / 513.0,
-845.0 / 4104.0,
-8.0 / 27.0,
2.0,
-3544.0 / 2565.0,
1859.0 / 4104.0,
-11.0 / 40.0,
]
}
fn b_coeffs() -> &'static [f64] {
&[
16.0 / 135.0,
0.0,
6656.0 / 12825.0,
28561.0 / 56430.0,
-9.0 / 50.0,
2.0 / 55.0,
25.0 / 216.0,
0.0,
1408.0 / 2565.0,
2197.0 / 4104.0,
-1.0 / 5.0,
0.0,
]
}
}