The hybrid function operational matrix together with the delay matrix are used for obtaining the solution of a linear time-varying system with delay in both the state and control. The method is based upon reducing the system to a set of algebraic equations.