One of the ecient and powerful schemes to solve linear and nonlinear equationsis homotopy analysis method (HAM). In this work, we obtain the approximate solution ofa system of partial dierential equations (PDEs) by means of HAM. For this purpose, wedevelop the concept of HAM for a system of PDEs as a matrix form. Then, we prove theconvergence theorem and apply the proposed method to nd the approximate solution ofsome systems of PDEs. Also, we show the region of convergence by plotting the H-surface.