> Any suggestions for solving a nonlinear system of ~1000 coupled ODEs with > complex coefficients?
Try your luck with the ODE solvers: ODE.jl, DASSL.jl and Sundials.jl. For sundials you will probably have to split the equations into real and imaginary part.
