[R] More accurate ODE solver?

[Tianchan Niu]

Dear All,

I was using rk4 and lsoda to solve a ODE system. However, both of them gave
bad accurate solutions, especially compared with Matlab solver ODE45. For
example, ODE45 gave solutions that can go to a stable level (about 1.6) when
time goes to infinity, however, the solutions from lsoda are decreasing to
very very small (about 1e-130) numbers. 

Does R have more accurate ODE solvers as ODE45 in Matlab?

Thank you,
Jane