We consider the class of real second-order linear dynamical systems that admit real diagonal forms with the same eigenvalues and partial multiplicities. The nonzero leading coefficient is allowed to be singular, and the associated quadratic matrix polynomial is assumed to be regular. We present a method and algorithm for converting any such n-dimensional system into a set of n mutually independent second-, first-, and zeroth-order equations. The solutions of these two systems are related by a real, time-dependent, and nonlinear n-dimensional transformation. Explicit formulas for computing the 2 n× 2 n real and time-invariant equivalence transformation that enables this conversion are provided. This paper constitutes a complete solution to the problem of diagonalizing a second-order linear system while preserving its associated Jordan canonical form.