%
% Written by:
% --
% John L. Weatherwax                2009-04-21
%
% email: wax@alum.mit.edu
%
% Please send comments and especially bug reports to the
% above email address.
%
% Problem on EPage 63; Solutions on EPage 161
%
%-----

close all; clc; clear all; 

addpath('../../BookCode');

epsilons = [ 0.1, 0.01, 0.001 ];
for ei=1:length(epsilons)
  e = epsilons(ei);
  coefs = [ -(13+e), (57+8*e), -(95+17*e), 50+10*e ];
  disp(e)
  bairstow( coefs, 4, 1.e-5 )
  
  roots( [1, coefs] )
end