Eigenvalues with numpy.linalg: normal modes of coupled oscillators

Afterwards you can turn a physical system into a matrix, solve it with numpy.linalg.solve, pick eigh or eig, and check eigenvectors by what they must satisfy.