Based upon the basic Properties of Legendre PoLynomials, the shifted Legendre delay operational matrix is introduced and then applied to the approximation of the solutions of linear delay systems with arbitrary time delay. An estimation of the parameter of the delay Systems is also reported in this paper.