% investigate the global convergence just by running run_mueller :
diary assignment_one_data
A = zeros(8,3);
for i = 1:3
  for d = 3:10
     p = poly(rand(1,d));
     [root,iter,resi,fail] = run_mueller(p,1.0e-12,50);
     A(d-2,i) = iter(1);
  end;
end;
A
sum(A)/8  % average with column sum
% investigate the local convergence by running mueller :
B = zeros(8,3);
for i = 1:3
  for d = 3:10
     r = rand(1,d);
     p = poly(r);
     x0 = r(1) + randn*1.0e-2;
     x1 = r(1) + randn*1.0e-2;
     x2 = r(1) + randn*1.0e-2;
     [root,iter,resi,fail] = mueller(p,x0,x1,x2,1.0e-12,15);
     B(d-2,i) = iter;
  end;
end;
B
sum(B)/8  % average with column sum
diary off;