a2 a1 b k1 1.0 -1.0 b Tau1 a3 a2 -1.0 c b Tau2 Cai Ca_dia iteration_time tstart a1 t k1 1.0 -1.0 t Tau1 Ca_dia iteration_time tstart iteration_time b tstart a2 -1.0 t b Tau2 Ca_dia iteration_time b tstart iteration_time c tstart a3 -1.0 t c Tau3 Ca_dia iteration_time c tstart