Abstract
This paper investigates the numerical solution of time‐fractional Volterra integrodifferential equations that arise in many real‐life problems. A higher‐order numerical scheme is proposed using the L1‐2‐3 method for the time‐fractional derivative, Simpson’s 1/3 rule approximates the integral part, and the spatial derivatives are approximated using the central finite difference method. The stability and convergence analysis of the proposed method is presented. The method’s efficiency is verified through illustrative examples, and the results are compared with those reported in other works.