Function fit_model 

pub fn fit_model(
    observations: &[&Observation],
    terms: &[Box<dyn Term>],
    fixed: &[bool],
    coefficients: &[f64],
    latitude: f64,
) -> Result<FitResult>