How can I find the solution of a differential equation, numerically solved with NDSolve, if there's a singularity on the curve ?
I know how to use NDSolve, but Mathematica gives a bad curve after some point (singularity jumping). I don't want to show this part, just the regular curve BEFORE the singularity.
More precisely, the curve should be strictcly positive : a[t] > 0. The NDSolve should stop the resolution if a <= 0. I added the command StoppingTest -> (a[t] < 0.001) or StoppingTest -> (a[t] <= 0) but it doesn't work. I'm still getting wrong curves with a[t] < 0.
Any idea ?