Implementing hybrid explicit/implicit ODE stepper

[Brian Gough]

At Wed, 3 Jun 2009 09:31:40 -0500,
Rhys Ulerich wrote:
> What's the "right way" to shoehorn this into the gsl_odeiv_system
> struc?  L(u) as gsl_odeiv_system.function with N(u) passed within
> gsl_odeiv_system.params ?  The reverse?  Neither?

That's a good question, I'm not sure - maybe Tuomo has some
suggestions? 

My gut instict would be to create a new class of solver for this but
it would be better if they could all be in one framework of course.