After we made the computer source code with mathematical model of Muramatsu et al. that was expressed by the set of simultaneous first-order ordinary differential equations in evaporation-pyrolysis zone of cigarette, we simulated the distribution profiles of temperature and density of flue-cured tobacco. Those equations were solved numerically with the Runge-Kutta-Gill algorithm assuming step size of 0.025mm by Muramatsu at at,, but in this study the advanced algorithm of Runge-Kutta 4th Order assuming step size of 0.0005mm. The initial conditions and physical parameters of Muramatsu et at. were used for solving them. The calculated values corresponded well with results of Muramatsu et al., especially the gradient of the temperature profile increased with smoldering speed and the thickness of the evaporation-pyrolysis zone decreased with increasing of smoldering speed. On the other hand, the temperature gradient decreased with increasing of the effective thermal-conductivity value and the thickness of the evaporation-pyrolysis zone increased with the effective thermal-conductivity value. Key words : cigarette, simulation, smoldering