[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...


More information about the R-help mailing list