function cost = SolvDep_2pathways_cost(par, data, constants) [CDcalc] = SolvDep_2pathways_simulator(par,constants,data); cost = (data(:,3)-CDcalc); if 1 == 0 figure subplot(1,2,1) plot(data(:,2),data(:,3),'o','MarkerSize',3) hold on plot(data(:,2),CDcalc) subplot(1,2,2) plot(data(:,2),cost) end end