I am using solve_ivp to solve a series of integrated differential equations. The below code is a simplified representation of a much larger code. I would like to be able to
