A numerical ODE solver that preserves the fixed points and their stability

In this paper, we provide a one-step predictor–corrector method for numerically solving first-order differential initial-value problems with two fixed points. The method preserves the stability behaviour of the fixed points, which results in an efficient integrator for this kind of problem. Some numerical examples are provided to show the good performance of the method.