ToolJuliaBeginner
Eigenvalues with LinearAlgebra: normal modes of coupled oscillators
- w2, V = eigen(Symmetric(K ./ m)) # eigenvalues are ω², in 1/s^2
- @printf("eigen(Symmetric(A)): %7.4f %7.4f %7.4f\n", eigvals(Symmetric(A))...)
- w2_co2, X = eigen(Symmetric(K1), M) # K a = ω² M a, ascending
- @printf("eigen(Symmetric(K1), M): %.4f %.4f %.4f\n", w2_co2...)
- w2 = eigvals(Symmetric(stiffness([0, k_bond, k_bond, 0])), Diagonal(masses_u .* u))