I am trying to solve a complex system of differential equations. The equations are stiff so I need to need to use a method which can handle both complex ODEs and stiffness switching. I have landed on using odeintw (link) which is able to solve my problem, however I also need to have event handling.