it’s the exactly same with RobustMultiNewton : Diverges at step 0.05 (exactly the same picture as here MethodOfLines - numerical instability for example in docs), and hangs at 0.025 (Julia continues to run in the background even is I “trash” it in the VSCode)