In this work, we present a modified generalized Mittag–Leffler function method (MGMLFM) and Laplace Adomian decomposition method (LADM) to get an analytic-approximate solution for nonlinear systems of partial differential equations (PDEs) of fractional-order in the Caputo derivative. We apply the MGMLFM and LADM on systems of nonlinear time-fractional PDEs. Precisely, we consider some important fractional-order nonlinear systems, namely Broer–Kaup (BK) and Burgers, which have found major significance because they arise in many physical applications such as shock wave, wave processes, vorticity transport, dispersal in porous media, and hydrodynamic turbulence. The analysis of these methods is implemented on the BK, Burgers systems and solutions have been offered in a simple formula. We show our results in figures and tables to demonstrate the efficiency and reliability of the used methods. Furthermore, our outcome converges rapidly to the given exact solutions.