> There is a rk routine in the contrib packages written by Jaime > Villate (I only have used it for 2-d vector fields). There are also the (undocumented) functions of package diffeq (share/numeric/diffeq.mac), which were written in 1980: runge1, runge2 and rungen. -- Mario Rodriguez Riotorto www.biomates.net