The nonequilibrium relaxation method is applied to first-order transition systems. Using the mixed-phase initialization, one can avoid the influence of the meta-stable states and the hysteresis phenomenon in the simulation. The transition point is estimated accurately on very large lattices. To...