function val = R(yExp, t, params)
 f2 = @(t,x) params(1)*t./2+params(2)*t-params(3);
 [~, yMod] = ode45(f2,t,0);
 val = yExp-yMod;
 val = sum(val);