Random-bond Ising chain in a magnetic field

The magnetization $m$ as a function of magnetic field $h$ for various temperatures is numerically obtained for the one-dimensional random-bond Ising model with each bond energy either $+J$ or $-J$ with equal probability. For high temperatures ($\frac{\mathrm{kT}}{J}>~0.6$) the results obtained agree with those of Landau and Blume. We have obtained $m$ vs $h$ for low temperatures as a numerical approach to the zero-temperature limit.