Abstract
A powerful analytical approach, namely the fractional residual power series method (FRPS), is applied successfully in this work to solving a class of fractional stiff systems. The methodology of the FRPS method gets a Maclaurin expansion of the solution in rapidly convergent form and apparent sequences based on the Caputo sense without any restriction hypothesis. This approach is tested on a fractional stiff system with nonlinearity ranging. Meanwhile, stability and convergence study are presented in the domain of interest. Illustrative examples justify that the proposed method is analytically effective and convenient, and it can be implemented in a large number of engineering problems. A numerical comparison for the experimental data with another well-known method, the reproducing kernel method, is given. The graphical consequences illuminate the simplicity and reliability of the FRPS method in the determination of the RPS solutions consistently.