% Script Eg5_2
% Illustrates the function EllipsePerimeter

  clc
  disp('  tol           b = 1.0       b = 0.1     b = 0.01    b = 0.001')
  disp('------------------------------------------------------------------')
  k = 10;
  for tol = logspace(-1,-k,k)
     fprintf('%8.1e',tol)
     for b = [1  .1 .01 .001]
        [P,n] = EllipsePerimeter(1,b,tol);
        fprintf('   %10d',n)
     end
     fprintf('\n')
  end

 
  

