Solving an eigenvalue problem in Python.

eigenvalues.py

import numpy as np

matrix = np.array([[0., 0.16, 0., -9.81],
                   [-0.1250, -0.5, 212., 0.],
                   [0., -0.0082, -0.5, 0.],
                   [0., 0., 1., 0.]])

print("Matrix:\n" + str(matrix))

# Get the eigenvalues (eig returns tuple with eigenvalues and eigenvectors)
eigenvalues = np.linalg.eig(matrix)[0]

print("\nEigenvalues:")
for eigenvalue in eigenvalues:
    print(eigenvalue)