function y=bisecc(f,a,b,M)
%Metodo de biseccion. Se parte del intervalo [a,b] y en M pasos se busca la
%raiz de la ecuacion f(x)=x
for i=1:M
    c=(a+b)/2;
    if f(c)==0
        y=c;
    elseif f(c)*f(a)>0
        a=c;
        c=(b+c)/2;
    else
        b=c;
        c=(a+c)/2;
    end
end
y=c;