ToolJuliaBeginner
Eigenvalues with LinearAlgebra: normal modes of coupled oscillators
- w2, V = eigen(Symmetric(K ./ m)) # eigenvalues are ω², in 1/s^2
- λ, W = eigen(A)
- @printf("eigen(A) : %7.4f %7.4f %7.4f, element type %s\n", λ..., eltype(λ))
- @printf("eigen(Symmetric(A)): %7.4f %7.4f %7.4f\n", eigvals(Symmetric(A))...)
- @printf("dot product of the first and last vector of eigen(A): %.3f\n", dot(W[:, 1], W[:, 3]))
- w2_co2, X = eigen(Symmetric(K1), M) # K a = ω² M a, ascending
and 1 more line in this tutorial