


%openf45, 1/2 2010, LE
x=[0.5 1.5 2.5 3 4 5]';
y=[3 1.5 1.5 1 1 0]';
plot(x,y,'ro')
pause
plot(x,y,x,y,'ro')
pause
A=[ones(size(x)) x x.^2 x.^3 x.^4 x.^5];
c=A\y;
xp=[0.5:0.01:5]';
Ap=[ones(size(xp)) xp xp.^2 xp.^3 xp.^4 xp.^5];
yp=Ap*c;
plot(xp,yp,x,y,'ro')

