Differential equation model for accommodation magnetization