function y=bisezione2(f,a,b,epsilon);
passo=0;
while (b-a)/2 >= epsilon
    c=f((a+b)/2);
    if c==0
        break
    end
    if f(a)*c<0
        b=(a+b)/2;
    else
        a=(a+b)/2;
    end
    passo=passo+1;
end
passo
y=(a+b)/2;