Software

To do the modelling, you need to use Python. To solve the system of ODEs you can use scipy.integrate.solve_ivp. To fit the model with the data, you can use lmfit.