function u=root2(p,q)
%
% u=root2(p,q)
%
%     Copyright: Eduardo Casas (Octubre, 1999)
%
% Esta función calcula las raices del polinomio x^2+px+q
%
if q==0,
    u=[-p;0];
else
    d=sqrt(p^2-4*q);
    if p<0,
        c=2*q/(d-p);
    else
        c=-2*q/(d+p);
    end
    u=[q/c;c];
end