The solution of point reactor kinetics is approximated by the rational approximation of the matrix exponential function. Three rational approximating methods which are Padé approximation, Quadrature Rational Approximation Method and Chebyshev Rational Approximation Method have been studied. And Richardson extrapolation technique has been used to improve the accuracy. The calculation results demonstrated that when the large positive reactivity is inserted, Quadrature and Chebyshev Rational Approximation Method are with better effects and stability, while Padé approximation with extrapolation may face the stability problem.