Abstract:We introduce a novel, time-efficient adaptive Runge-Kutta computational scheme tailored for systematically solving linear and nonlinear Volterra-type Integro-Differential Equations (VTIDEs). This scheme is particularly effective for equations featuring a specific class of convolution memory, K ( t , s ) = a κ ( t - s ) within convolution integrals, where a , κ ∈ R . Such equations frequently arise in fields like viscoelasticity, population dynamics, epidemiology, control systems, financial mathematics, and neuroscience, where current states are influenced by historical data, especially in dynamical systems within physics. Our proposed scheme demonstrates remarkable computational efficiency, achieving a cost-effective computation time of O ( N t ) for N t iterations, a significant improvement over existing schemes for VTIDE, which typically exhibit complexities of O ( N t log ⁡ N t ) and O ( N t ( log ⁡ N t ) ) 2 . The key feature of our scheme lies in its ability to elegantly and efficiently handle convoluted integrals of varying complexity such as those including nonlinearities and derivatives. Specifically, our approach eliminates the need for backstage integration, allowing the main Runge-Kutta scheme to operate efficiently with only a few data points at each time step. This not only saves time but also prevents potential data storage issues while rendering numerical implementation straightforward and tractable. Moreover, this article introduces flexible discretizations for integration and differentiation, accommodating unequally spaced data points. We implement a fast RK45-Fehlberg algorithm (TE-RK45-Fehlberg), which we use to solve three examples of nonlinear VTIDE presented in this study. We further demonstrate the effectiveness of our approach by using an example, where we employ a novel parabolic step size generator function to efficiently generate step sizes for a computationally fast RK4 scheme (Flex-RK4). For all examples tackled using TE-RK45-Fehlberg and Flex-RK4, we generate plots that compare numerical data with exact solutions, providing insight into the accuracy and reliability of our approach. We also show how the scheme provides straightforward ways of transforming it into implicit forms and obtaining stability curves.

Efficiently and easily integrating differential equations with JiTCODE, JiTCDDE, and JiTCSDE

$D$-Module Techniques for Solving Differential Equations in the Context of Feynman Integrals

D-module techniques for solving differential equations in the context of Feynman integrals

DelayDiffEq: Generating Delay Differential Equation Solvers via Recursive Embedding of Ordinary Differential Equation Solvers

Applying stiff integrators for ODEs and DDEs to problems with distributed delays

Collocation in Conjunction with the Local Variational Iteration Method

Rapid variable-step computation of dynamic convolutions and Volterra-type integro-differential equations: RK45 Fehlberg, RK4

Highly Accurate and Efficient Time Integration Methods with Unconditional Stability and Flexible Numerical Dissipation

Some robust integrators for large time dynamics

Symbolic Computation of Jacobi Elliptic Function Solutions to Nonlinear Differential-Difference Equations

Numerical integration of ordinary differential equations with rapidly oscillatory factors

Numerical simulation of differential-algebraic equations with embedded global optimization criteria

Exponential Integrators for Phase-Field Equations using Pseudo-spectral Methods: A Python Implementation

New efficient numerical methods for some systems of linear ordinary differential equations

Numerical Derivative-based Flexible Integration Algorithm for Power Electronic Systems Simulation Considering Nonlinear Components

A tutorial on solving ordinary differential equations using Python and hybrid physics-informed neural network

A Family of Conditionally Explicit Methods for Second-Order ODEs and DAEs: Application in Multibody Dynamics

Differentiable Programming for Differential Equations: A Review

An Adaptive Parallel Integrator of Ordinary Differential Equations System for Space Experiment Simulation

Dual stepsize explicit numerical integration method and applications