function y=bisezione(effe,a,b,k)
%ritorna i valori di k iterate successive di bisezione per @effe su [a,b]
y=zeros(k,1);
if effe(a)*effe(b)>0
   error 'intervallo non valido' 
end

for i=1:k
    t=(a+b)/2;
    if effe(t)==0
        y(i:end)=t;
        break;
    elseif effe(t)*effe(a)>0
        a=t;
    else
        b=t;
    end
    y(i)=(a+b)/2;
end
