Out of scientific curiosity, I am interested in Lagrangian dynamical systems that can be expressed in a "linear" manner. By this, I mean that their Lagrangian can be expressed, linearly, as
L = 1/2 (dq/dt)^T.M.(dq/dt) - 1/2 q^T.K.q^T,
where 'q' is a vector in R^n; 'M' and 'K' are symmetric, positive-definite, n-dimensional square matrices with constant coefficients.
One example of such a system is the "linear" double pendulum which is simply a double pendulum where one assumes small angles during motion (see for example https://math24.net/double-pendulum.html).
Usually, it is easy to obtain the analytic solutions for this kind of system (as for the linear double pendulum). So my question is if there exist any other examples such as this one, with true technical/scientific interest? And particularly, is there a similar example where the analytic solution is not obvious at all (requiring the use of numerical methods)?