Whenever I get problems with the solver hitting dtmin, I try force_dtmin=true
res = step(P, 5, force_dtmin=true)
display(plot(res));
I don’t see why this particular problem would be hard on the solver, but the force_dtmin hack appears to have worked.