ToolJuliaBeginner
Eigenvalues with LinearAlgebra: normal modes of coupled oscillators
- residual = maximum(abs.(K ./ m * V - V * Diagonal(w2)))
- @printf("largest entry of (K/m) V - V Diagonal(w2): %.0e 1/s^2\n", residual)
- S = V * Diagonal(w2) * V' # K/m put back together from its modes
- M = Diagonal(m_co2) # u
- w2 = eigvals(Symmetric(stiffness([0, k_bond, k_bond, 0])), Diagonal(masses_u .* u))
- println("K1 * X ≈ M * X * Diagonal(w2_co2): ", K1 * X ≈ M * X * Diagonal(w2_co2))