Exemplo n.º 1
0
def eigen_calc(M):
    x = Symbol('x')
    M = Matrix(M)
    I = Matrix([[x, 0, 0], [0, x, 0], [0, 0, x]])
    MI = M - I
    char_equation = MI.det()
    a = Poly(char_equation, x)
    eigen_values = np.roots(a.coeffs())
    return eigen_values