This work presents the application of the differential transform method (DTM) to the model of pollution for a system of three lakes interconnected by channels. Three input models (periodic, exponentially decaying, and linear) are solved to show that DTM can provide analytical solutions of pollution model in convergent series form. In addition, we present the posttreatment of the power series solutions with the Laplace-Padé resummation method as a useful strategy to extend the domain of convergence of the approximate solutions. The Fehlberg fourth-fifth order Runge-Kutta method with degree four interpolant (RKF45) numerical solution of the lakes system problem is used as a reference to compare with the analytical approximations showing the high accuracy of the results. The main advantage of the proposed technique is that it is based on a few straightforward steps and does not generate secular terms or depend of a perturbation parameter.