bisect.m

various matlab programs to slove various problems

function [xvect,xdif,fx,nit]=bisect(a,b,toll,nmax,fun)
err=toll+1; 
nit=0; 
xvect=[]; 
fx=[]; 
xdif=[];
while (nit < nmax & err > toll)
    nit=nit+1; 
    c=(a+b)/2; 
    x=c; 
    fc=eval(fun); 
    xvect=[xvect;x]; 
    fx=[fx;fc]; 
    x=a; 
    if (fc*eval(fun) > 0)
      a=c; 
    else 
      b=c; 
    end; 
    err=abs(b-a); 
    xdif=[xdif;err];
end
return