Non-polynomial series solution of a coupled first order partial and ordinary differential equation (PDE-ODE) system for a discontinuous initial and boundary condition has been developed. Linear equation systems are constructed to calculate the constant coefficients of the series solution. Explicit expressions have been found to the solution of these linear equation systems. Different forms of the solution have been compared to the numerical solution of the PDE-ODE system and the rate of the convergence is also investigated. The studied first order PDE-ODE system describes an unsteady convection dominated heat transfer process induced by a buoyant plume entrainment.