Kurkcu, Omur Kivanc2025-07-102025-07-1020250894-33701099-1204https://doi.org/10.1002/jnm.70061https://hdl.handle.net/20.500.13091/10143This study is concerned with treating the fractional integral equations with functional kernels and variable delays, introducing a hyperaccurate semi-analytical method based on the Stieltjes-Wigert polynomials, matrix expansions, and the Laplace transform. After analytically converting the terms in the governing equation into the matrix expansions of the Stieltjes-Wigert polynomials type at the collocation points, the method gathers these matrices into a unique matrix equation and then readily solves it by an elimination technique. The residual improvement technique is also introduced to correct the obtained solutions. The residual error bound analysis is theoretically proved via algebraical properties and the mean value theorem for fractional integral calculus, respectively. Six model equations are treated via the method, which runs on a devised computer program. Based on the outcomes, the method is straightforward to treat model equations and to encode its mainframe on a mathematical software.eninfo:eu-repo/semantics/closedAccessError Bound AnalysisIntegral EquationsMatrix-Collocation MethodRiemann-Liouville Fractional Integral OperatorStieltjes-Wigert PolynomialsArticle10.1002/jnm.700612-s2.0-105007425523