| euler {SoilR} | R Documentation |
This function can solve arbitrary first order ode systems with the euler
forward method and an adaptive time-step size control given a tolerance for
the deviation of a coarse and fine estimate of the change in y for the next
time step. It is an alternative to deSolve.lsoda.wrapper and
has the same interface. It is much slower than ode and should probably be
considered less capable in solving stiff ode systems. However it has one
main advantage, which consists in its simplicity. It is quite easy to see
what is going on inside it. Whenever you don't trust your implementation of
another (more efficient but probably also more complex) ode solver, just
compare the result to what this method computes.
euler(times, ydot, startValues)
times |
A row vector containing the points in time where the solution is sought. |
ydot |
The function of y and t that computes the derivative for a given point in time and a column vector y. |
startValues |
A column vector with the initial values. |