We plug the values for the rate constants into the solved a, b, and c functions.

Input := 


cfit[t_] = csol[t]/.solk[[2]];

bfit[t_] = bsol[t]/.solk[[2]];

afit[t_] = asol[t]/.solk[[2]];