%curvefit09.m

x=[1	2	3	4	5	6	7	8	9	10];
y=[2615 1943	1494	1087	765	538	484	290	226	204];
plot(x,y,'.k','markersize',20)

n=3;
a=polyfit(x,y,n)
x0=4.5;
y0=polyval(a,x0)
x1=11;
y1=polyval(a,x1)
x2=1:0.1:10;
hold on
y2=polyval(a,x2);
plot(x2,y2,'-r','markersize',20)
xlabel('x')