% 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)