Abstract
In this paper, we broaden the utilization of a beautiful computational scheme, residual power series method (RPSM), to attain the fractional power series solutions of nonhomogeneous and homogeneous nonlinear time-fractional systems of partial differential equations. This paper considers the fractional derivatives of Caputo-type. The approximate solutions of given systems of equations are calculated through the utilization of the provided initial conditions. This iterative scheme generates the fast convergent series solutions with conveniently determinable components. The implementation of this numerical scheme clearly exhibits its effectiveness, reliability and easiness regarding the procedure of the solution, as well as its better approximation. The repercussions of the fractional order of Caputo derivatives on solutions are depicted through graphical presentations for various particular cases.