Our study proposes a MATLAB simulation to analyze the hysteresis loss in single-phase transformers under a frequency of 50 Hz, calculating the variables and the relationship between the magnetic flux density (B) and the magnetic field intensity (H) in the theoretical equation to design used in numerous conditions, including consideration of parameter of the ferromagnetic silicon steel core material, sizes, power rated, and amplitude of alternating current (AC) sinusoid that affects the B-H curve, and saturations point in the time domain. Comparing the hysteresis loss characteristics between the model and empirical experiments and the experiment test results are given in this paper.