On the Construction of Exact Numerical Schemes for Linear Delay Models

Exact numerical schemes have previously been obtained for some linear retarded delay differential equations and systems. Those schemes were derived from explicit expressions of the exact solutions, and were expressed in the form of perturbed difference systems, involving the values at previous delay intervals. In this work, we propose to directly obtain expressions of the same type for the fundamental solutions of linear delay differential equations, by considering vector equations with vector components corresponding to delay-lagged values at previous intervals. From these expressions for the fundamental solutions, exact numerical schemes for arbitrary initial functions can be proposed, and they may also facilitate obtaining explicit exact solutions. We apply this approach to obtain an exact numerical scheme for the first order linear neutral equation x ′ ( t ) − γ x ′ ( t − τ ) = α x ( t ) + β x ( t − τ ) , with the general initial condition x ( t ) = φ ( t ) for − τ ≤ t ≤ 0 . The resulting expression reduces to those previously published for the corresponding retarded equations when γ = 0 .