This contribution proposes a variational symplectic integrator aimed at linear systems issued from the least action principle. An internal quadratic finite-element interpolation of the state is performed at each time step. Then, the action is approximated by Simpson’s quadrature formula. The implemented scheme is implicit, symplectic, and conditionally stable. It is applied to the time integration of systems with quadratic Lagrangians. The example of the linearized double pendulum is treated. Our method is compared with Newmark’s variational integrator. The exact solution of the linearized double pendulum example is used for benchmarking. Simulation results illustrate the precision and convergence of the proposed integrator.