How do I numerically solve an ODE in Python?
Consider
\\ddot{u}(
scipy.integrate() does ODE integration. Is that what you are looking for?