ToolPythonBeginner Eigenvalues with numpy.linalg: normal modes of coupled oscillators D = K1 / np.sqrt(np.outer(m_co2, m_co2)) # K_ij / sqrt(m_i m_j)D = stiffness([0, k_bond, k_bond, 0]) / np.sqrt(np.outer(mass, mass))