1 | initial version |
A good source of information is this question: http://ask.sagemath.org/question/8714/ode_solver-unable-to-convert-to-float/
Your best bet is likely either to convert to real differential equations or use Scipy's scipy.integrate.ode zvode
which can be done within Sage.