debug in calling odeint, and add a test code