Help debugging differential equation

[ChrisRackauckas]

@Ge0rges

@ChrisRackauckas I wonder if this is an acceptable solution to issue with LU factorization you highlighted?

yes, but instead of setting it to a small value, set it to a real zero.

zero_condition(u, t, integrator) = u[1] < 1e-12
function zero_out!(integrator)
        integrator.u[i] = 0
end