function bisect(f,a,b,tol,n)
% Bisection method for solving the nonlinear 
%equation  f(x)=0.
a0=a;
b0=b; 
iter=0;
u=feval(f,a);
v=feval(f,b);
c=(a+b)*0.5;
err=abs(b-a)*0.5;
disp('_____________________________________________________________________')
disp(' iter    a            b           c           f(c)      |b-a|/2  ')
disp('_____________________________________________________________________')
fprintf('\n')
if (u*v<=0)
   while (err>tol)&(iter<=n)
      w=feval(f,c);
      fprintf('%2.0f  %10.4f  %10.4f  %12.6f  %10.6f  %10.6f\n',iter,a,b,c,w,err)
      if (w*u<0)
         b=c;v=w;
      end
      if (w*u>0)
         a=c;u=w;
      end
      iter=iter+1;
      c=(a+b)*0.5;
      err=abs(b-a)*0.5;
   end 
   if (iter>n)
      disp('  Method failed to converge')
   end   
else
   disp('   The method cannot be applied f(a)f(b)>0')  
end
% Plot f(x) in the interval [a,b].
fplot(f, [a0 b0])
xlabel('x');ylabel('f(x)');
grid