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 Mc and valence Mv bands. Mv was insensitive to temperature and exhibited sub-linear field dependence, which is consistent with the previous research on Dirac electrons. By contrast, Mc 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 kBT ≃ EF , where EFis the Fermi energy measured from the band bottom and kBis the Boltzmann constant. At low temperatures kBT ≲ EF , M exhibited sub-linear field dependence, whereas M exhibited super-linear field dependence at high temperatures kBT ≳ EF . This qualitative change in the field dependence of M will play a significant role in the magnetization of Dirac electrons with small EF .