[R] ODE's in R
Markku Karhunen
markku.karhunen at ktl.fi
Tue Jan 22 12:01:35 CET 2008
Thanks, Dr. Maechler.
> No, there's no such track.
> [ Matlab users coming to R may produce wrong R code
> by using 0:n-1 instead of 0:(n-1) ; but I don't assume this
> would be the case ]
>
>
Been there, done that!
> MK> We use just a simple discretisation written in a for loop
> MK> and a variable (i.e. user-fed) time step.
>
> I don't think you should use your own code instead of "professional"
> ODE solvers, such as the one in R package 'odesolve'....
>
>
We must look into that. The problem, maybe, is that in fact half of the
equations are, in fact, simple PDE's and I don't know, if you can put
them into odesolve.
> MK> Maybe, I'm too neurotic about this, but I guess I just want some comfort
> MK> after seeing a few particularly nasty orbits.
>
> As we know ``from Chaos theory'', there can be delicate
> inhereent and numerical problems in ODE solving..
>
But - to our best knowledge - they should not be any more acute in R,
than on any other platform...
BR,
Markku
More information about the R-help
mailing list