% a quick demo of Horner's method and its effects

clear all
close all

% first a comparison of ways to compute

x = 0.988:.0001:1.012;
y1 = x.^7-7*x.^6+21*x.^5-35*x.^4+35*x.^3-21*x.^2+7*x-1;
plot(x,y1,'r--')
pause

% because this polynomial can be nicely factored ...

y2 = (x-1).^7;
plot(x,y2)
pause

% now let's see Horner's method

y3=-1+x.*(7+x.*(-21+x.*(35+x.*(-35+x.*(21+x.*(-7+x))))));
plot(x,y3)
pause

% now for the head-to-head comparison

plot(x,y1,'r--',x,y3)
