ode_solver : unable to convert to float
I am trying to solve a system of 38 first-order differential equations using ode_solver(). Among many problems, I have a problem with complex numbers. The functions in question are complex functions (equations include conjugate(function)), and initial contitions are also complex numbers.
The error I get is following:
TypeError: Unable to convert 0.0353553390593274 - 0.0353553390593274*I
to float; use abs() or real_part() as desired
0.0353553390593274 - 0.0353553390593274*I is the value of one of the inital conditions. I am rather new to sage, so I wonder if ode_solver() can deal with complex functions and complex numbers as such?
Thank you!
If you have a specific example, that would be really helpful. If it's quite long, use http://pastebin.com or something to link to.
It's REALY long, and it contains lots of abbreviations and definitions within definitions... I doubt it would be very useful. I can describe the system: I have 39 functions (y[0], ..., y[38]) and 39 dependent first-order diferential equations. Functions are defined in the complex ring, the inital values are complex numbers, etc. RHS of equations are complicated combinations of y's. So i write them all down, define range and initial conditions and that's it. In Mathematica I did it with NDSolve, there were absolutly no problems (I have just typed the equations) and the whole job is done in about 2 seconds. Now I am trying to do the same job with Sage. But it turns out it's not that simple :(
Can you maybe post a much smaller test case that illustrates the problem?