Entanglement entropy encodes important features of strongly interacting quantum many-body systems and gauge theories, but its analytical study is still limited to systems with high level of symmetry. This motivates the search for efficient techniques to investigate this quantity numerically, through Monte Carlo calculations on the lattice. In this contribution, we discuss the computation of the entropic c-function by means of an algorithm based on Jarzynski's equality, which is an exact theorem from non-equilibrium statistical mechanics. After presenting benchmark results for the Ising model in two dimensions, where our algorithm successfully reproduces the analytical predictions from conformal field theory, we discuss its generalization to the three-dimensional Ising model, for which we were able to extract universal terms beyond the area law. Finally we point out some future generalizations of this calculation.