x=[-2,0,1,3];
y=[2.0782,2.0427,5.2494,25.1510];
PolynomialOrder=2;
[P,S]=polyfit(x,y,PolynomialOrder);
[Y,DELTA]=polyval(P,x,S);
errorbar(x,Y,DELTA,'rd-')
BestFitString=[];
for ii=1:PolynomialOrder+1
BestFitString=[BestFitString num2str(P(ii),'%+g') '*x^' num2str(PolynomialOrder-ii+1)];
end
title(['Best fit polynomial is : y=' BestFitString])
xlabel('x')
ylabel('y')
grid on
