AlexNascou

Excellent suggestions. I'll give it a try when I get a chance and report back.

Sorry for the delay in getting back. I looked at the location of the discontinuity again and you're right, it doesn't correspond to initTauU. I'm not sure why that is...

I played around more with the tolerances and dtmax and I now have it working very well. I still find the behaviour a bit strange however. Using Rodas5, I can...