function [c, r, s] = Itera_Bairstow(a, r0, s0)

N = length(a);
res =[r0, s0]';

for k=1:100
[b, c] = Division_Sintetica_Derivadas(a, res(1), res(2));

A = [c(3) c(4); c(2) c(3)];
d = [-b(2); -b(1)];
disp('b=');
disp(b');
disp('c=');
disp(c');

disp(A);
disp(d);

disp('res = ');
res = res - inv(A)*d;
disp(res);
if sqrt(norm(d)) < 0.0001 break; end
end;

r = res(1);
s = res(2);
c = b(3:N);