The classic model system for this is the Ising model, which I discussed in an earlier post. The difference here is that we vary the external field and see how the system's internal field reacts. The typical plot looks like this:
|Based on Wikipedia|
I decided to adapt my previous Ising script to try to demonstrate this effect, and I was surprised by my success: Hysteresis.py
Marika and I are packing things up to return to the States in a couple weeks, so I may miss posting.