We study the realization problem for linear time-invariant systems described by higher-order differential equations, which are J-lossless with constant J. Our approach is based on the factorization of a constant matrix obtained in a straightforward way from the storage function of the system. State equations are obtained directly from this factorization and from the original system representation.