ode23tb stops solving

Hello,

I'm using ode solvers for my nonlinear problem in the form of:

Mx''+Cx'+Kx=F

in which M, C, and K are matrices and F is a vector, with nonlinear terms. I tried all ODE solvers, but they ran very slowly, and the best and fastest one was ODE23tb.

This solver was 100 times faster than others, but in some conditions of the problem, if the time interval is [0 10] for example, ode23tb quickly runs up to 5th second (for example in 1 minute) but stops at this second and do not go ahead, and spends one minute for every 0.001 seconds. It seems suddenly the time steps reduces to 10e-6.

I changed different options of ODE but were not effective.

May you help me?

Thanks.

[Jan]

If the function to be integrated contains a pole at this point or a discontinuity, this is the expected behavior.

[Sadegh A]

Dear jan

Thank you for your comment.

No, there is not any pole or discontinuity.

May it has any other reason?

[Jan]

Without seeing any details it is impossible to guess the cause of the problem. Can you post your code?

Answer 1

[Nikhil Sonavane]

You may refer to the following documentation which will help you choose a solver for your problem more accurately-

• Choose a Solver