First of all, there is a big difference between symbolic and numeric solvers. Sympy is a symbolic library, and while it has some facilities to do numerics, if numerics is all that you will be doing (as you imply) you should use numpy/scipy which are created for that.
For your case check out http://docs.scipy.org/doc/scipy/reference/generated/scipy.integrate.odeint.html However, do not forget that you need to rewrite your second order equation as a vector of two first order equations: u' = v v' = u'' = f(v, u) -- You received this message because you are subscribed to the Google Groups "sympy" group. To unsubscribe from this group and stop receiving emails from it, send an email to sympy+unsubscr...@googlegroups.com. To post to this group, send email to sympy@googlegroups.com. Visit this group at http://groups.google.com/group/sympy?hl=en-US. For more options, visit https://groups.google.com/groups/opt_out.