function [b, c] = Factoriza(a, s0)
% Factoriza(a, s0) Factoriza una función polinomial con coeficientes
% a y valor inicial s0

    s = s0;
    while(1)
        s = s + Evalua_Polinomio(-s, a)/Evalua_Derivada(-s, a);
        if abs(Evalua_Polinomio(-s, a)) < 1e-6
            break;
        end
    end;
    
    fprintf('s = %f\n', s);
    
    [b, r] = Division_Sintetica(a, s);
    c = [s, 1];
end