Stability of a saddle node bifurcation under numerical approximations

In this paper, we show that the solution flows generated by a one-parameter family of ordinary differential equations are stable under their numerical approximations in a vicinity of a saddle node. Our result sharpens the one in [1] and the proof is adapted from the method of Sotomayor in [2., 3.].