Abstract
In this paper, we propose an algorithm to find approximate solutions of the proposed system of the fractional heat-like equations. The proposed algorithm, basically, illustrates how the two powerful algorithms, the Homotopy Perturbation Method (HPM) and the Sumudu Transform Method (STM), can be combined and used to get exact solutions of fractional partial differential equations. We also present some examples to illustrate the accuracy and the effectiveness of this algorithm.