Function dsc_model 

fn dsc_model(t: f64, p: &Fit, b0: f64) -> (f64, f64)

The authors’ dscModel.m: the C_minus compartment with its frequency negated ((2/3)·(−χ−)·γ·B0), the rest at the baseline decay.