function y=bisezione(f,a,b,k)
if f(a)*f(b)>=0
    error('L''intervallo non contiene uno zero');
end

for passo=1:k
    [a b]
    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
end
y=(a+b)/2;