Function levinson 

pub fn levinson(r: &[f64]) -> Result<Levinson>

Run Levinson–Durbin on autocovariances r[0..=K].