function [ alpha ] = bis (f,a,b)
% bis(f,a,b) metodo di bisezione
m = (a + b)/2;
for k = 1:100
if f(m) == 0 alpha = m; else
    if sign(f(a)) == sign(f(m)) a = m; 
	    else b = m;
		end;
		m= (a + b)/2;
end;
alpha = m;
end
