My code solves a 1st order differential equation using Scipy’s solve_ivp function and the standard RK45 solver. The code itself is reasonably large, and I will certainly attempt to produce a minimal working example if all else fails. However, there might already be a problem at a conceptual that stage, so I am hoping this is not a bug but just my lack of understanding.