Analysis of the numerical scheme of the one-dimensional fractional Rayleigh-Stokes model arising in a heated generalized problem

In this paper, we present a well-organized method to estimate the one-dimensional fractional Rayleigh-Stokes model using the construction of orthogonal Gegenbauer polynomials (GBPs) and Lagrange square interpolation to estimate the time derivative. Therefore, we design an authentic and fast numerical calculation approach based on the elaborated convergence rate recovery method. The matrix of the derivative operation of an orthogonal GBP is gained by employing the characteristic of this type of polynomial. The privilege of the numerical method is the orthogonality of the GBP and operational matrices, which reduces time computation and increases speed. Eventually, we propose three numerical examples to check the validity and numerical studies to illustrate the precision and efficiency of the new approach.