% Rootfinding and Condition Number
% An Example
echo on
x = 1:8
a = poly(x)
X = roots(a)
ahat = a; ahat(9)=40325
Xhat = roots(ahat);
[X Xhat]
relouterror = (Xhat-X)./X
relinerror = 5/40320
relouterror/relinerror