Improved Crank-Nicolson method for modelling diffusion of electromagnetic fields in hysteretic soft magnetic laminations