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