add RandomTamedEM() solver for RODEs

[rmsrosa]

This implements the RODE equivalent of the "Tamed Euler" SDE solver proposed in Strong convergence of an explicit numerical method for SDEs with nonglobally Lipschitz continuous coefficients by Martin Hutzenthaler, Arnulf Jentzen, Peter E. Kloeden.

This scheme, for a RODE $dX_t/dt = f(t, X_t, Y_t)$, with noise $Y_t$, is an explicit scheme consisting of the time steps
$$X_{t_{j+1}}^n = X_{t_j}^n + \frac{\Delta t f(t_j, X_{t_j}, Y_{t_j})}{1 + \Delta t ||f(t_j, X_{t_j}, Y_{t_j})||}$$