NDSolve to solve a system of coupled ODEs for a range of different parameters. For most points in the parameter space, the integration breaks down prematurely. However, I would still like to plot the resulting
InterpolatingFunctions up to the point where
NDSolve broke down.
Of course, that point
NDSolve changes based on the chosen parameters so I find myself constantly having to adjust the plot bounds manually.
It seems to me that there should be an easier way. Can I tell Mathematica to simply plot an
InterpolatingFunction for all values inside the domain (i.e. without using extrapolation)?