Function stirling_accurate 

fn stirling_accurate(dx: f64, parity: f64, f_res: DoubleDouble) -> DoubleDouble