Use np.linalg.eig or scipy.sparse.linalg.eigs for large systems.
Most dynamic systems (vibrations, chemical reactions) are described by ODEs. Runge-Kutta Methods (RK4) . Numerical Methods In Engineering With Python 3 Solutions