Finally the RP equation was solved in its dimensional form with DPRKN6(). For low values of acoustic pressure (< 1.28e5 Pa), there is no problem.
For higher values, the collapse is very sharp and (of course) I get:
Warning: dt <= dtmin. Aborting. There is either an error in your model specification or the true solution is unstable.
What solver would you suggest. Is there an example of use?
Thank you in advance,
Thierry
PS: the code was changed on 2 points: (1) + 2*σ/R₀ (instead of -), (2) (R₀/R)^(4.0).