TempExactFun = @(time)14.6*exp(-1.1*time);
timeMeasured = linspace(0,2,41);
TempMeasured = round(TempExactFun(timeMeasured));
CoefLinFit=polyfit(timeMeasured, TempMeasured,1);
time=0.7;
exact=TempExactFun(time)
TempLinFit=polyval(CoefLinFit,time)
timePlot=linspace(0,2,200);
TempExactPlot=TempExactFun(timePlot);
TempLinFitPlot=polyval(CoefLinFit,timePlot);
plot(timePlot,TempExactPlot,'--k',...
    timeMeasured,TempMeasured,'ok',...
    timePlot,TempLinFitPlot,'g')