Novel derivative operational matrix in Caputo sense with applications

The main objective of this study is to present a computationally efficient numerical method for solving fractional-order differential equations with initial conditions. The proposed method is based on the newly developed generalized derivative operational matrix and generalized integral operational matrix derived from Laguerre polynomials, which belong to the class of orthogonal polynomials. Through the utilization of these operational matrices, the fractional-order problems can be transformed into a system of Sylvester-type matrix equations. This system is easily solvable using any computational software, thereby providing a practical framework for solving such equations. The results obtained are compared against various benchmarks, including an existing exact solution, Podlubny numerical techniques, analytical and numerical solvers, and reported solutions from stochastic techniques employing hybrid approaches. This comparative analysis serves to validate the accuracy of our proposed design scheme.