clear;
close all;
clc;

a=5;
b=3;
c=0.1;

f1 = @(t,x) a*t./2+b*t-c;
t=linspace(0,10,100);
[~, yExp] = ode45(f1,t,0); % Generate experimental data

paramsGuess=[4.9;2.9;0.1];
optFun = @(params) R(yExp,t,params);
[x,fval] = fminsearch(optFun,paramsGuess) % optimize