General nonstandard finite-difference schemes for differential equations with three fixed-points

We construct general nonstandard finite-difference (NSFD) schemes that provide explicit methods for first-order differential equations with three fixed-points y′(t)=±y(y−α)(y−1) where 0≤α≤10≤α≤1. For y≥0y≥0, these methods, regardless of the step-size chosen, are stable with respect to the monotonicity of solutions and are elementary stable. That is, they preserve the critical properties of the original differential equation such as the positivity of the solutions, the stability behavior of all fixed-points, and the monotonicity of solutions within each subinterval (0,α)(0,α), (α,1)(α,1), and (1,∞)(1,∞).