A new approximation scheme for first-order ordinary differential equations with non-negative right-hand sides