Ode45, evolution giving NaNs after some time, with no evidence of singularities in the equations
1 view (last 30 days)
Show older comments
Mohammad Farhat
on 7 Apr 2020
Commented: Walter Roberson
on 8 Apr 2020
Hello,
Attached is a simple integrator for four coupled first order ode s.
after a certain time, the 4-vector solution turns into NaNs with no evidence of singularities in the equations.
This happens after a certain range where the solutions (starting as real) are written in complex notation although the imaginary part is zero. (why?)
can anyone please have a look and let me know if there's something im missing?
2 Comments
Accepted Answer
Walter Roberson
on 7 Apr 2020
Your first order, e1, gets to be more than 1.0, and then the sqrt(1-x^2) gets to be complex. Problems accumulate from there.
What are the original equations that you are trying to deal with? Do you have the Symbolic Toolbox ?
6 Comments
Walter Roberson
on 8 Apr 2020
I do not think an event function is appropriate here. Projection and failed steps is normal for ode45() .
Setting a MaxStep in the options might potentially help.
More Answers (0)
See Also
Categories
Find more on Loops and Conditional Statements in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!