oefEsts2 = nlinfit(time, log(conc), @(p,x)log(modelFun(p,x)), startingVals); line(xgrid, modelFun(coefEsts2, xgrid), 'Color',[0 .5 0]); legend({'Raw Data' 'Additive Errors Model' 'Multiplicative Errors Model'});