1pub fn sl9(x: f64) -> f64 {
2 let pi = core::f64::consts::PI;
3 let pi2 = pi*pi;
4 let pi4 = pi2*pi2;
5 let pi6 = pi4*pi2;
6 let pi8 = pi4*pi4;
7 let x2 = x*x;
8 x*(pi8/9450.0 + x2*(-1.0/5670.0*pi6 + x2*(pi4/10800.0 + x2*(-1.0/30240.0*pi2 + (pi/80640.0 - 1.0/725760.0*x)*x))))
9}