Hello,
Im working on a differential equation based model some of its equations has power term (x^a) which is major issue because without this model works fine. Now issue is that I have used stiff solver(RodasP5), different tolerence levels (1e-2 - 1e-15), maxiters, passing dt
but still its of no help.
I see NaN issues, negative values although system should be always positive. So, I also used callbacks which helps with positivity but still system is very unstable and just produces NaN. Solver fails even at t = 0.0.
Can you guys suggest something or point me to any related discussion if you know about.
Thank you