clear,clc,clf
x=[0 1 3 4 6 7 9 10]';
y=[0 1.2 2 2.4 2.5 2.2 1.4 0.3]';
plot(x,y,'*')
hold on
A=[ones(8,1) x x.^2];
c=A\y
xlin=0:0.1:10;
ylin=c(1)+c(2)*xlin+c(3)*xlin.^2;
plot(xlin,ylin)