a method for simulating a concentration of impurities within a semiconductor device while the semiconductor device is being oxidized includes setting a transition region in a part of the semiconductor device which is oxidized for a unit time period , and then calculating the concentration of impurities within the transition region by solving a diffusion equation using an impurity diffusion coefficient peculiar to the transition region , thus obtaining a redistribution of impurities within the transition region .