title: scipy.integrate.solve_ivp unclear how to solve implicit ODE of form 0=F(t, y(t), y'(t))
body: <p>Currently, I do use <a href="https://jmodelica.org/assimulo/tutorial_imp.html" rel="nofollow noreferrer">assimulos solver suite</a> to solve an implicit differential equation of the form 0 = F(t, y(t), y'(t)). I would like to use a solver which comes with a native scipy installation and had a look at the documentation for <a href="https://jmodelica.org/assimulo/tutorial_imp.html" rel="nofollow noreferrer">scipy.integrate.solve_ivp</a>. The documentation says that the solver solves a a system of the form:</p>
  
  <blockquote>
    <p>Solve an initial value problem for a system of ODEs.
    This function numerically integrates a system of ordinary differential equations given an initial value <code>dy / dt = f(t, y)</code>    <code>y(t0) = y0</code></p>
  </blockquote>
  
  <p>I do miss the ability to solve a system of the form <code>dy / dt = f(t, y, y')</code>. Is this possible with <code>scipy.integrate.solve_ivp</code> ?</p>