Non-equilibrium evolution at absolute temperature T and approach to equilibrium of statistical systems in long-time (t) approximations, using both hierarchies and functional integrals, are reviewed. A classical non-relativistic particle in one spatial dimension, subject to a potential and a heat bath (hb), is described by the non-equilibrium reversible Liouville distribution (W) and equation, with a suitable initial condition. The Boltzmann equilibrium distribution Weq generates orthogonal (Hermite) polynomials Hn in momenta. Suitable moments Wn of W (using the Hn’s) yield a non-equilibrium three-term hierarchy (different from the standard Bogoliubov–Born–Green–Kirkwood–Yvon one), solved through operator continued fractions. After a long-t approximation, the Wn’s yield irreversibly approach to equilibrium. The approach is extended (without hb) to: (i) a non-equilibrium system of N classical non-relativistic particles interacting through repulsive short range potentials and (ii) a classical ϕ4 field theory (without hb). The extension to one non-relativistic quantum particle (with hb) employs the non-equilibrium Wigner function (WQ): difficulties related to non-positivity of WQ are bypassed so as to formulate approximately approach to equilibrium. A non-equilibrium quantum anharmonic oscillator is analyzed differently, through functional integral methods. The latter allows an extension to relativistic quantum ϕ4 field theory (a meson gas off-equilibrium, without hb), facing ultraviolet divergences and renormalization. Genuine simplifications of quantum ϕ4 theory at high T and large distances and long t occur; then, through a new argument for the field-theoretic case, the theory can be approximated by a classical ϕ4 one, yielding an approach to equilibrium.