% Metodo da bissecçao
% Resolução de ex. da aula teorico-pratica
% Carlos Balsa - Métodos Numéricos 2007/08
%
a=1; b=3; tol=0.5e-2;
fa=a^2-4*sin(a);
k=0;
erro=b-a;

while ((erro)>tol) 
    k=k+1;
    m=(a+b)/2;
    fm=m^2-4*sin(m);
    if sign(fa)==sign(fm)
        a=m;
    else
        b=m;
    end
    erro=b-a;
end

k         % nº total de iterações
x=(a+b)/2 % raiz aproximada
erro      % erro absoluto