Various numerical methods are studied to solve the point kinetics equations in nuclear dynamics. However the point kinetics equation has stiffness problem. So it is difficult to obtain accurate results. In this paper, we proposed a method based on the Taylor-Lie series method under the First order hold (FOH) assumption. To evaluate the performance of presented method, the procedure is tested using multi input signals. By applying the method to the equation, we can resolve the stiffness problem.