Abstract
In this study, solutions of system of time-space fractional partial differential equations (FPDEs) are obtained by utilizing the Shehu transform iterative method. The utility of the technique is shown by getting numerical solutions to a large number of system of FPDEs.