biseccion

function raiz=biseccion(def,x0,x1,tol)
f=inline(def);
if f(x0)*f(x1)<0
    x=x0;
    c=1;
    while abs(f(x))>tol
        c=c+1;
        x=(x0+x1)/2;       
        if f(x0)*f(x)<0
            x1=x;
        else
            x0=x;
        end       
    end
    raiz=x;   
else
    raiz='La raiz no se puede calcular porque no cambia de signo';   
end
fprintf('\n\nNumero de iteraciones: %d',c);
end