In this paper, a method to construct an analytic-numerical solution for homogeneous parabolic coupled systems with homogeneous boundary conditions of the type ut = Auxx, A1u(o,t) + B1ux(o,t) = 0, A2u(1,t) + B2ux(1,t) = 0, ...