Análisis Numérico - 21/05/2016 - Método de bisección

% Programa para calcular raices utilizando el algoritmo de bisección

a = input('Ingrese el límite inferior del intervalo: ');

b = input('Ingrese el límite superior del intervalo: ');

n = input('Ingrese la cantidad de iteraciones deseadas: ');

% Argumento s se utiliza para una entrada con caracteres

func = input('Ingrese la función: ', 's');

% inline se utiliza para declarar la variables de una función ingresada

f = inline(func);

for i = 1 : n
    c = (a+b)/2;

    e = abs((a-b)/2);

    A(i,:) = [i a b c f(c) e];

    if f(a)^f(c) < 0
        b = c;
    else
        a = c;
    end
end

disp('    n         a         b         c         f(c)      error');

disp(A);