np.linalg.eig

Вот цитата из рассыли [1]:

eig() is for nonsymmetric matrices and eigh() is for symmetric
(or hermitian matrices).

The former most likely will return complex eigen values. The latter
always returns real eigen values.

eig() is slower and numerically less reliable (it calls LAPACK's DGEEV
and ZGEEV). Even for symmetric matrices it can produce complex
eigen pairs due to numerical round-off.

linalg.eig -- функция вычисляющая собственные векторы и собственные значения матрицы

Пример:

np.linalg.eig(np.matrix('-1 -6; 2 6'))

Формат возвращаемых значений

Например для матрицы:

np.linalg.eig(np.matrix('-1 -6; 2 6'))

получим: