Normal forms for real linear Hamiltonian systems with purely imaginary eigenvalues

This note gives a concise algorithm for computing a normal form for a real linear Hamiltonian differential equatin which has purely imaginary eigenvalues. This algorithm is then applied to the differential equation which comes from the quadratic terms of the Hamiltonian of the restricted three body problem at a Lagrange equilateral triangle equilibrium point.