function [eigv, eigd, eigvalues, theta1, theta2, theta1D] = eigSym2(C)
if nargin < 1
    C = [0.9025, 0.125;0.125, 1.0001];
end
[eigv, eigd] = eig(C);
eigv(:,1) = eigv(:,1) / norm(eigv(:,1));
eigv(:,2) = eigv(:,2) / norm(eigv(:,2));
if (eigd(1,1) < eigd(2,2))
    tmp = eigv;
    eigv(:,1) = eigv(:,2);
    eigv(:,2) = tmp(:,1);
    
    tmp = eigd;
    eigd(1,1) = eigd(2,2);
    eigd(2,2) = tmp(1,1);
end
eigvalues(1) = eigd(1,1);
eigvalues(2) = eigd(2,2);
theta1 = atan(eigv(2,1)/ eigv(1,1));
theta2 = theta1 + pi/2;
theta1D = theta1 * 180 / pi;