The Intel Math Kernel Library has some subroutine to solve coupled ODE's using Gear's Method?
I solved this once ago using DIVPAG subroutine on IMSL Library, but now i'm trying to solve this without use IMSL or any non-free library.
MKL does not have any routines, using Gear's method. There was some investigation into support for ODE's done awhile ago and this was made available via our whatif program. You can take a look at the archive of that project, but you should note that it is no longer actively supported.