On numerical modelling of hysteresis in recording systems