[Rcpp-devel] C++ code for ordinary differential equations

[Douglas Bates]

I am working with pharmacokinetic/pharmacodynamic models, some of
which can only be represented as the solutions of ordinary
differential equations.  I am looking for ODE sofware written in C++.
The deSolve package for R has lots of wonderful functions but I may
want to take things apart a bit and I really don't have the patience
for wading through Fortran code to try to discover what it is doing.

I can, of course, google likely phrases but if anyone has experience
with such solvers and can make any recommendations I would appreciate
it (off-list responses are probably better than on-list)