{ x1 = (1 - 1/(k1*(exp(k2*t) - 1) + 1))*k3*k4*k5*(k8 - k5)/k5 }