Abstract
This paper demonstrates an effective and powerful technique, namely fractional He-Laplace method (FHe-LM), to study a nonlinear coupled system of equations with time fractional derivative. The FHe-LM is designed on the basis of Laplace transform to elucidate the solution of nonlinear fractional Hirota-Satsuma coupled KdV and coupled mKdV system but the series coefficients are evaluated in an iterative process with the help of homotopy perturbation method manipulating He's polynomials. The fractional derivatives are considered in the Caputo sense. The obtained results confirm the suggested approach is extremely convenient and applicable to provide the solution of nonlinear models in the form of a convergent series, without any restriction. Also, graphical representation and the error estimate when compared with the exact solution are presented.