In this study, we evaluated the dependence of magnetization of three-dimensional Dirac electrons in the quantum limit on the magnetic field and temperature. The magnetization was calculated by differentiating the free energy with respect to the magnetic field. The field and temperature dependence of the chemical potential were entirely considered under the canonical ensemble condition. The total magnetization M consisted of two contributions from the conduction and valence bands. was insensitive to temperature and exhibited sub-linear field dependence, which is consistent with the previous research on Dirac electrons. By contrast, was sensitive to both temperature and magnetic field, yielding a non-trivial contribution to the total M. As a result, the properties of total M changed at approximately , where is the Fermi energy measured from the band bottom and is the Boltzmann constant. At low temperatures , M exhibited sub-linear field dependence, whereas M exhibited super-linear field dependence at high temperatures . This qualitative change in the field dependence of M will play a significant role in the magnetization of Dirac electrons with small .