We formulate a least squares function in the constants (now variables) k1 and k2 using the data from the chemical c.

Input := 


SSc[k1_,k2_] = Sum[(csol[cdata[[n,1]]] -  cdata[[n,2]])^2,

                 {n,1, Length[cdata]}];

Input := 


solck = FindMinimum[SSc[k1,k2],{k1,.1},{k2,.2}]

Output =


{0.000135185, {k1 -> 0.0210623, k2 -> 0.0640486}}