The Spectral Homotopy Analysis Method Extended to Systems of Partial Differential Equations