Solving ODE that uses a max() function

[ijustlovemath]

After much, much debugging on my end, and further examination of how the equations should work, I think I've got it. Thanks so much for the help!

[johnsold]

Can you post a brief summary of the things you learned or found so that someone who comes across this thread in the future can learn from your experience?

Thanks,

Lynn

[ijustlovemath]

Sure thing.

I found that the biggest bugs were squashed by rectifying unit discrepancies between A pump, B pump, and the units of the differential equation. This also reduced the solutions sensitivity to h. As always, unit tests and modularity helped a lot with debugging.

Finally, I was using the wrong values for X0, so I shot an email to the author of the paper in question to get reasonable estimates.