Representation of a solution of the Cauchy problem for an oscillating system with pure delay

We obtain a representation of a solution of the Cauchy problem for a linear inhomogeneous differential equation with constant coefficients and pure delay. We use special matrix functions called a delayed matrix sine and a delayed matrix cosine. They have the form of matrix polynomials of degree dependent on the value of delay.