ODE with complex entires

Question

I am trying to solve a system of ODEs with complex entries. From the GSL documentation it follows that it accepts only real entries. Is there any way to pass complex ones (anything more straightforward than distinguish real and imaginary parts)? In case it is impossible, could you please advise any other good library for this purpose?

Answer 1

You can also use odeint. It ships with boost 1.53 and allows you directly to use complex numbers.