 Models of ferromagnetic hysteresis can be developed using thermodynamics. These models must satisfy certain conditions, including the second law of thermodynamics, Clausius-Duhem inequality, and Euclidean invariance. Additionally, the entropy production must be non-negative for all admissible thermodynamic processes. This implies that the entropy production must take on a specific form when describing different behaviors during loading and unloading. In the special case of a one-dimensional setting, a detailed model is derived for the magnetization function, which depends on a given susceptibility function. Starting from different initial magnetized states, hysteresis cycles can then be obtained by solving a nonlinear ordinary differential equation. Cycles with both large and small amplitudes can be generated from materials like soft iron. This article was authored by Claudio Giorgi and Angelo Moro.