Thank you! I have tried the new solvers TRBDF2() , KenCarp4() , or QNDF() but I got the same issue reported in the first figure. I then increased the tolerance with soln = solve( prob, AutoVern7(Rodas5()), callback=delay, abstol=1e-20, reltol=1e-20, tstops=[t_in] ) and I got it through:
So the case is closed.
However, I think the problem is not so much the solver or the ODEs (since they work given certain values) but how to pick the correct parameters. This is the same issue as in my previous post: how to get the range of values for the parameters?
Ideally, I would need a program that will output something of this sort:
↓quantity →time
1 2 3 …
10 No No Yes …
20 No Yes Yes …
30 No Yes No …
…
Is it possible to get such a list?
Thanks