d_parent = - (alpha/beta) * 1/((time/beta) + 1) * parent

d_parent = - ((k1 * g * exp(-k1 * time) + k2 * (1 - g) * exp(-k2 * time)) / (g * exp(-k1 * time) + (1 - g) * exp(-k2 * time))) * parent

d_parent = - ifelse(time <= tb, k1, k2) * parent

d_A1 = + f_parent_to_A1 * k_parent * parent - k_A1 * A1

d_parent = - (alpha/beta) * 1/((time/beta) + 1) * parent
d_A1 = + f_parent_to_A1 * (alpha/beta) * 1/((time/beta) + 1) * parent - k_A1 * A1

d_parent = - ifelse(time <= tb, k1, k2) * parent
d_A1 = + f_parent_to_A1 * ifelse(time <= tb, k1, k2) * parent - k_A1 * A1

d_parent = - k_parent * parent
d_A1 = + f_parent_to_A1 * k_parent * parent - k_A1 * A1
d_B1 = + f_parent_to_B1 * k_parent * parent - k_B1 * B1

d_parent = - (alpha/beta) * 1/((time/beta) + 1) * parent
d_A1 = + f_parent_to_A1 * (alpha/beta) * 1/((time/beta) + 1) * parent - k_A1 * A1
d_B1 = + f_parent_to_B1 * (alpha/beta) * 1/((time/beta) + 1) * parent - k_B1 * B1

d_parent = - ifelse(time <= tb, k1, k2) * parent
d_A1 = + f_parent_to_A1 * ifelse(time <= tb, k1, k2) * parent - k_A1 * A1
d_B1 = + f_parent_to_B1 * ifelse(time <= tb, k1, k2) * parent - k_B1 * B1

d_parent = - (alpha/beta) * 1/((time/beta) + 1) * parent
d_A1 = + f_parent_to_A1 * (alpha/beta) * 1/((time/beta) + 1) * parent - k_A1 * A1
d_A2 = + f_A1_to_A2 * k_A1 * A1 - k_A2 * A2

d_parent = - ifelse(time <= tb, k1, k2) * parent
d_A1 = + f_parent_to_A1 * ifelse(time <= tb, k1, k2) * parent - k_A1 * A1
d_A2 = + f_A1_to_A2 * k_A1 * A1 - k_A2 * A2

d_parent = - k_parent * parent
d_A1 = + f_parent_to_A1 * k_parent * parent - k_A1 * A1
d_B1 = + f_parent_to_B1 * k_parent * parent + f_A1_to_B1 * k_A1 * A1 - k_B1 * B1

d_parent = - (alpha/beta) * 1/((time/beta) + 1) * parent
d_A1 = + f_parent_to_A1 * (alpha/beta) * 1/((time/beta) + 1) * parent - k_A1 * A1
d_B1 = + f_parent_to_B1 * (alpha/beta) * 1/((time/beta) + 1) * parent + f_A1_to_B1 * k_A1 * A1 - k_B1 * B1

d_parent = - ifelse(time <= tb, k1, k2) * parent
d_A1 = + f_parent_to_A1 * ifelse(time <= tb, k1, k2) * parent - k_A1 * A1
d_B1 = + f_parent_to_B1 * ifelse(time <= tb, k1, k2) * parent + f_A1_to_B1 * k_A1 * A1 - k_B1 * B1