  Probability distributions of the magnetic work are computed for the 2D Ising
model by means of Monte Carlo simulations. The system is first prepared at
equilibrium for three temperatures below, at and above the critical point. A
magnetic field is then applied and grown linearly at different rates.
Probability distributions of the work are stored and free energy differences
computed using the Jarzynski equality. Consistency is checked and the dynamics
of the system is analyzed. Free energies and dissipated works are reproduced
with simple models. The critical exponent $\delta$ is estimated in an usual
manner.
