Control Theory & Systems
Shiva Karimi; Elyas Shivanian; Zahra Barikbin
Abstract
This paper presents a novel numerical method for solving variable-order fractional integro-differential equations using two-dimensional fractional-order Fibonacci wavelets. The proposed approach employs fractional-order Fibonacci wavelets together with their associated integral and derivative operational ...
Read More
This paper presents a novel numerical method for solving variable-order fractional integro-differential equations using two-dimensional fractional-order Fibonacci wavelets. The proposed approach employs fractional-order Fibonacci wavelets together with their associated integral and derivative operational matrices. First, new integral and derivative operational matrices are derived. These matrices, which exhibit improved accuracy in the numerical examples reported herein, are then employed to transform the governing equation into a system of algebraic equations. The collocation method is subsequently applied to solve this system and determine the unknown coefficients. Finally, error analysis, convergence results based on relevant theorems, and numerical examples are provided to demonstrate the accuracy, reliability, and efficiency of the proposed method.

Mathematics & Theoretical Foundations
Saad Qasim Abbas; Wasan Saad Ahmed
Abstract
This paper presents a systematic comparative study of two widely used numerical solvers --- HOFiD_bvp (high-order finite difference scheme) and bvp4c (collocation-based) --- for solving singular second-order ordinary differential equations (ODEs) with first-kind (regular) boundary singularities. Four ...
Read More
This paper presents a systematic comparative study of two widely used numerical solvers --- HOFiD_bvp (high-order finite difference scheme) and bvp4c (collocation-based) --- for solving singular second-order ordinary differential equations (ODEs) with first-kind (regular) boundary singularities. Four representative benchmark problems drawn from fluid dynamics, materials science, and radially symmetric diffusion models are used to evaluate solver performance across key metrics: maximum residual, maximum error, mesh point count, and ODE/BC function call counts. Results show that HOFiD_bvp consistently achieves lower residuals and errors with fewer function evaluations, making it computationally more efficient. Conversely, bvp4c demonstrates superior robustness for nonlinear singular problems and offers better adaptive mesh refinement capabilities. These findings provide practical guidance for selecting the appropriate numerical technique in applied science and engineering contexts, with implications for optimization of computational simulation workflows.
