I realize this is an old post, but this might be helpful to others:
One reason for EIG to return complex values with very small imaginary part, could be that A is very close to, but not exactly, symmetric. In this case,
[U, D] = eig( ( A + A')/2 );
will make EIG treat the input as symmetric, and always return real eigenvalues.