The spatial discretization form of the space-dependent reactor kinetics equation is a first-order simultaneous ordinary differential equation in time. Conventional numerical methods of the space-dependent kinetics equation, i.e., the generalized Runge-Kutta method, the implicit method (backward Euler method), and the Theta method, are based on the time difference approximation. However, the present study adopts the analytical solution of the space-dependent kinetics equation expressed by the matrix exponential and no time difference approximation is used. In this context, our present approach is classified as an explicit method in which no iteration calculation on space and energy is necessary. The Krylov subspace method is used to evaluate the matrix exponential observed in the solution of the spatially discretized space-dependent kinetics equation. The Krylov subspace method is implemented into a space-dependent kinetics solver. In order to examine the effectiveness of the Krylov subspace method, the TWIGL benchmark problem is analyzed as a verification calculation. The calculation results show the effectiveness of the present method especially in the step reactivity perturbation.