For this problem the solution is trivial, but I'd like you to include code here that goes from the Kane object to a ode in first order form that can be integrated with odeint (from scipy). This will make the example much easier to extend for multiple bodies. So please use the dummify, lambdify, odeint method to get your simulation results instead of writing in out explicitly.