implemented simple scipy odeint algorithm for testing purposes