Hello all,

I have a system of two differential equations that I am able to solve using `Rosenbrock23()`

.

I now would like to model this system with one modification, where at time `t`

I manually change the value of one of the differential equation results, causing a discontinuity.

I thought that simply passing `t`

to d_discontinuties would handle this, but I run out of iterations, or the solver produces very unusual results. I currently have it manually set up to solve from `0`

to `t-1`

then from `t`

to `f`

starting with values from `t-1`

, but that does not seem like an elegant solution and I get the following error:

`┌ Warning: dt(-1.6e-9) <= dtmin(1.862645149230957e-9) at t=1.059225e7. Aborting. There is either an error in your model specification or the true solution is unstable.`

Any tips for the proper way to handle this? My equations are stiff.