title: How does odeint() from scypy python module work?
body: <p>I am a physics student interested in solving ODEs numerically. I usually write my own solvers in C using Runge–Kutta methods.</p>
  
  <p>I recently learned Python, and I used SciPy’s <code>odeint</code> function to solve ODEs. But I am worried about how the function algorithm works, because <strong>it not take a step size argument</strong>. So, how can I learn how it works? How can I know what is the precision of their results?</p>
  
  <p>I consulted <a href="https://docs.scipy.org/doc/scipy/reference/generated/scipy.integrate.ode.html#scipy.integrate.ode" rel="noreferrer">this documentation</a>,
  but it does not offer very much information, and I don’t really understand the optional arguments that they describe.</p>