# Solve an ODE with positivity-preserving property unconditionally

I have an ODE for a scalar function $u=u(t)$ of the form: $$\frac{du}{dt}=L(u).$$ Here the function $L=L(u)$ satisfies: $$L(0)=0, \quad L'(u)\le0.$$ Then it is easy to see that the solution $u=u(t)$ has the following property:

(i) If the initial value $u(0)\ge0$, then $u(t)\ge0$ for any $t>0$;

(ii) If the initial value $-M\le u(0)\le M$ with $M>0$, then $-M\le u(t)\le M$.

Question: I need to find a numerical method to solve this ODE and under arbitrary step size maintain the two properties or only one of two. I only know that the Euler backward method has these two properties. Does anyone of you know high order (e.g. third-order) methods to solve this problem?

Any links to literature or for further reading would be greatly appreciated.

• @DavidKetcheson The solution of an ODE has the property of monotonicity-preserving, i.e., if $u(t)$ and $v(t)$ are solutions to the same ODE and $u(0)\le v(0)$, then $u(t)\le v(t)$ for any $t>0$. Since $u\equiv0$ is a solution to this ODE, property (i) holds for both positivity-preserving and negativity-preserving with the aid of monotonicity-preserving property. – Michael Jun 16 '16 at 18:34