function [p,q,x]=itcomplex(A,x)
%
% [p,q,x]=itcomplex(A,x)
%
%     Copyright: Eduardo Casas (Octubre, 1999)
%
% Esta función realiza una iteración para generar una sucesión
% (p_k,q_k) convergente hacia los coeficientes del polinomio 
% t^2+pt+q cuyas raices son los dos valores propios complejos
% conjugados dominantes de la matriz A.
%
% Significacion de las variables de entrada:
%
%    A     Matriz nxn real cuyos valores propios dominantes son dos
%          numeros complejos conjugados
%
%    x     vector real nx1 en el que se inician las iteraciones
%
% Significacion de las variables de salida:
%
%    (p,q) Una iteración del método
%
%    x     Vector que debe utilizarse en entrada para la siguiente
%          iteración
%    
y=A*x;y=[y A*y];y=[y A*y(:,2)];y=[y A*y(:,3)];
[v,k]=max(abs(y(:,4)));
y=y/y(k,4);
M=y(k,2)^2-y(k,3)*y(k,1);
N=y(k,1)-y(k,3)*y(k,2);
M1=y(k,3)^2-y(k,2);
p=N/M;q=M1/M;
x=y(:,4);