I have a system of second order differential equation to be solved numerically. I would like to set up events to halt integration to find the values of phi when r(phi)=2/3 . Here is my code
The code only works for the second event so it halts for r(phi)=1.63... etc
How do i stop this?
Thanks for any help.