Sorry for the lack of a minimal working example but the code I am concerned about is rather large and I have no idea what the problem is. All I know is that when I run the code in Mathematica versions 9, 10, 11, version 9 gives the correct result while versions 10 and 11 do not. There are no error messages. My best guess is that the problem is related to changes to NDSolve between the different versions since I have seen a similar problem before. However, with previous problems, I was able to stumble across a fix by playing with different choices for Method. This time I am having no such luck.
So my question is, is there a way to force NDSolve in version 11 to work precisely like the default settings of version 9?