The systems of reaction-diffusion equations coupled with moving boundaries defined by Stefan condition have been widely used to describe the dynamics of spreading population. There are several numerical difficulties to efficiently handle such systems. Firstly extremely small time steps are usually demanded due to the stiffness of the system. Secondly it is always difficult to efficiently and accurately handle the moving boundaries. In this talk, to overcome these difficulties, we first transform the one-dimensional problem with a moving boundary into a system with a fixed computational domain, and then introduce four different temporal schemes: Runge-Kutta, Crank-Nicolson, implicit integration factor (IIF) and Krylov IIF for handling such stiff systems. Numerical examples are examined to illustrate the efficiency, accuracy and consistency for different approaches, and it can be shown that Krylov IIF is superior to other three approaches in terms of stability and efficiency by direct comparison. I will also briefly review the integration factor methods and their applications in viscous fluid flows with moving interfaces.