Estimation parameters using GEKKO (a second-order differential equation)

Question

I am using a second-order parallel model:

It has four parameters (initial values Cf0 and Cs0, kf and ks) to be estimated using GEKKO and six experiments. However