Add a convenience function to increase the field and calculate magnetization and work
Add calculation of work as the integral of magnetization