Hurin's profile - activity

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 :(

Thank you for your answer, I had a feeling something like that will be the case :) Can you advise me with the best way to do that conversion? Do I have to do it "by hand" (separating each function f = fRe + fIm, and each value like y0 = y0Re + y0Im) , or is there some routine that can do it?

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!