ODE solvers with complex data

Hi again,

My first task with GSL is to convert some Matlab code to C. The code uses
the Matlab ODE solvers, and so I'm planning on using the GSL odeiv routines.
However, the Matlab code uses complex data, but it seems the GSL odeiv
routines only use double-type data. The GSL odeiv system "function" routine
specifies input and output as type double. I figured that since GSL supports
complex data types, there would be matching odeiv functions for complex
types if it were possible to use them.

I'm new to ODE solvers, so I could be missing something. Does anyone have
any ideas about this?

Cheers,
Michael