Numerical Solution of the Linear Fractional Delay Differential Equation Using Gauss–Hermite Quadrature

Fractional order differential equations often possess inherent symmetries that play a crucial role in governing their dynamics in a variety of scientific fields. In this work, we consider numerical solutions for fractional-order linear delay differential equations. The numerical solution is obtained via the Laplace transform technique. The quadrature approximation of the Bromwich integral provides the foundation for several commonly employed strategies for inverting the Laplace transform. The key factor for quadrature approximation is the contour deformation, and numerous contours have been proposed. However, the highly convergent trapezoidal rule has always been the most common quadrature rule. In this work, the Gauss–Hermite quadrature rule is used as a substitute for the trapezoidal rule. Plotting figures of absolute error and comparing results to other methods from the literature illustrate how effectively the suggested approach works. Functional analysis was used to examine the existence of the solution and the Ulam–Hyers (UH) stability of the considered equation.