Equivalent Systems for Differential Equations with Polynomially Distributed Delay
Roland Pulch
TL;DR
This paper tackles delay differential equations with polynomially distributed delays by deriving an exact equivalent system that uses only two discrete delays, enabling standard DDE techniques to be applied; it also introduces a Gaussian-quadrature discretisation to approximate the distributed delay as a DDE with multiple discrete delays. The authors establish the auxiliary-function framework $x_i(t)=\int_a^b y(t-\tau) \tau^i d\tau$ and show how to replace the integral with $\sum_{i=0}^n \alpha_i x_i(t)$, leading to a coupled system whose stability is analyzed; a scaling $\tilde{t}=t/b$ is proposed to mitigate error amplification. Generalisations to systems and to multiple distributed delays are provided, along with a quadrature-based alternative that preserves the distributed-delay dynamics in approximate form. Numerical experiments on a nonlinear test model (SIR-type) compare both approaches, showing qualitative agreement and a convergence of quadrature-based solutions to the exact distributed-delay solution as the number of nodes grows, with time-step accuracy playing a crucial role in the observed error growth.
Abstract
We consider delay differential equations with a polynomially distributed delay. We derive an equivalent system of delay differential equations, which includes just two discrete delays. The stability of the equivalent system and its stationary solutions are investigated. Alternatively, a Gaussian quadrature generates a discretisation of the integral, which describes the distributed delay in the original delay differential equation. This technique yields an approximate differential equation with multiple discrete delays. We present results of numerical computations, where initial value problems of the differential equations are solved. Therein, the two approaches are compared.
