In this paper we develop an approach for obtaining the solutions to systems of linear retarded and neutral delay differential equations. Our analytical approach is based on the Laplace transform, inverse Laplace transform and the Cauchy residue theorem. The obtained solutions have the form of infinite non-harmonic Fourier series. The main advantage of the proposed approach is the closed-form of the solutions, which are capable of accurately evaluating the solution at any time. Moreover, it allows one to study the asymptotic behavior of the solutions. A remarkable discovery, which to the best of our knowledge has never been presented in the literature, is that there are some particular linear systems of both retarded and neutral delay differential equations for which the solution approaches asymptotically to a limit cycle. The well-known method of steps in many cases is unable to obtain the asymptotic behavior of the solution and would most likely fail to detect such cycles. Examples illustrating the Laplace transform method for linear systems of DDEs are presented and discussed. These examples are designed to facilitate a discussion on how the spectral properties of the matrices determine the manner in which one proceeds and how they impact the behavior of the solution. Comparisons with the exact solution provided by the method of steps are presented. Finally, we should mention that the solutions generated by the Laplace transform are, in most instances, extremely accurate even when the truncated series is limited to only a handful of terms and in many cases becomes more accurate as the independent variable increases.