{ diff(x1, t) = (-1)*k14*k1*x2*x1/k14,
  diff(x2, t) = ((-1)*k14*k1*x2*x1 + 1*k14*k2*x3 + (-1)*k14*k7*k9*(k10 + (1 - k10)*(1 - exp(-t/k11)^3))*x2 + 1*k14*k4*x4)/k14,
  diff(x3, t) = (1*k14*k1*x2*x1 + (-1)*k14*k2*x3 + (-1)*k14*k5*k7*x3*x6 + 1*k14*k6*x5 + (-1)*k14*k7*k9*(k10 + (1 - k10)*(1 - exp(-t/k11)^3))*x3)/k14,
  diff(x4, t) = (1*k14*k7*k9*(k10 + (1 - k10)*(1 - exp(-t/k11)^3))*x2 + 1*k14*k7*k9*(k10 + (1 - k10)*(1 - exp(-t/k11)^3))*x3 + 1*k14*k9*(k10 + (1 - k10)*(1 - exp(-t/k11)^3))*x5 + (-1)*k14*k4*x4)/k14,
  diff(x5, t) = (1*k14*k5*k7*x3*x6 + (-1)*k14*k6*x5 + (-1)*k14*k9*(k10 + (1 - k10)*(1 - exp(-t/k11)^3))*x5)/k14,
  diff(x6, t) = (-1)*k14*k5*k7*x3*x6/k14,
  diff(x7, t) = (1*k14*k7*k9*(k10 + (1 - k10)*(1 - exp(-t/k11)^3))*x3 + 1*k14*k9*(k10 + (1 - k10)*(1 - exp(-t/k11)^3))*x5 + (-1)*k14*k8*x7)/k14 }