Hi, I have a question for mathematica programming. I want to repeatedly use NDSolve to solve a system of high order ODEs, which depend on a few parameters. I wanna write some kind of loop structure to scan through certain ranges of the parameters. But in most cases NDSolve will crash at some point of the independent variable (say r0), either due to numertical issues or possible singularity. I want to extract the point where NDSolve crashes, and feed this value (r0) to some other operation i want to carry out, but in an automatic way (within the loop structure). Is there anyway I can do this? Thank you.