Tool Python Beginner 40 min
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.