Python: SciPy Curve_Fit for Fitting System of ODEs to Data for Chemical Reaction Kinetics

Question

I am trying to build a way to fit a system of ODEs to my mock chemical reaction data. I am trying to find a way to build a tool to find certain parameters (rate constants e.g. k