TempExactFun = @(time)14.6*exp(-1.1*time);
TempMeasured = [15     8     5     3     2];
timeMeasured = [0 .5 1 1.5 2];
CoefQuadFit=polyfit(timeMeasured, TempMeasured,2);
time=0.7;
exact=TempExactFun(time);
TempQuadFit=polyval(CoefQuadFit,time);
timePlot=linspace(0,2,200);
TempExactPlot=TempExactFun(timePlot);
TempQuadFitPlot=polyval(CoefQuadFit,timePlot);
plot(timePlot,TempExactPlot,'--k',...
    timeMeasured,TempMeasured,'ok',...
    timePlot,TempQuadFitPlot,'g')