Mathematica 9 is now available
Student Support Forum
Student Support Forum: 'How to stop a curve when singularity in NDSolve ?' topicStudent Support Forum > General > Archives > "How to stop a curve when singularity in NDSolve ?"

Help | Reply To Topic
Author Comment/Response
02/14/12 09:25am

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 ?

URL: ,
Help | Reply To Topic