| 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.
Copyright Sage, 2010. Some rights reserved under creative commons license. Content on this site is licensed under a Creative Commons Attribution Share Alike 3.0 license.