Function cos_accurate 

fn cos_accurate(
    y: DoubleDouble,
    msin_k: DoubleDouble,
    cos_k: DoubleDouble,
) -> f64