Metodo De Secante

function p = LinBairstow(a,varargin)

% raices = LinBairstow(a,EPS=eps)
%
% Obtiene todas las raíces (reales y complejas) delpolinomio P(x):
%
% P(x) = a(1) x^n + a(2) x^n-1 + ... + a(n-1) x + a(n+1)

% 24 sept 2007
% M. Valenzuelaglobal RASTREO
RASTREO = 1;

if length(varargin)>=1
EPS = varargin{1};else
EPS = eps;
end

p = [];
n = length(a)-1;
r = 0; s = 0; % valores iniciales de ry s
if RASTREO
fprintf('r inicial:%f s inicial:%f\n\n',r,s);
end
while n>2
b = divSintetica(a,r,s);c = divSintetica(b(1:n),r,s);
A = [c(n-1) c(n-2) -b(n);c(n) c(n-1) -b(n+1)];
Delta = Montante(A);
r = r + Delta(1);s = s + Delta(2);
if RASTREO
display(A);
fprintf('Delta r=%f, Delta s=%f\n',Delta);
fprintf('r=%fs=%f\n\n',r,s);
pause
end
if sum(abs(Delta))<EPS
sq = sqrt(r^2+4*s);
p = [p; (r+sq)/2;(r-sq)/2];
a = b(1:n-1);
n = n - 2;
r = 0; s = 0; % valores iniciales de r y s
if RASTREOfprintf('r inicial:%f s inicial:%f\n\n',r,s)
end
end
end
if n==1
p = [p; -a(2)/a(1)];...