Parameter optimization for system of ODEs using scipy.differential_evolution fails due to solve_ivp throwing runtime_warnings

Question

I\'m trying to solve a model in python and to fit unknown parameters of the model to experimental data. The model consists of 2 ODEs and I solve it using scipy.integrate.sol