{ diff(x1, t) = ((-1)*k1*k2*k15*x1 + 1*k1*k8*x17 + 1*k1*k10)/k1,
  diff(x2, t) = 0,
  diff(x3, t) = (1*k1*k2*k15*x1 + (-1)*k1*k3*x5*x3)/k1,
  diff(x4, t) = 0,
  diff(x5, t) = ((-1)*k1*k3*x5*x3 + 1*k1*k8*x17 + 1*k1*k9)/k1,
  diff(x6, t) = 0,
  diff(x7, t) = (1*k1*k3*x5*x3 + (-1)*k1*k4*k18*k15*x7)/k1,
  diff(x8, t) = 0,
  diff(x9, t) = 0,
  diff(x10, t) = 0,
  diff(x11, t) = (1*k1*k4*k18*k15*x7 + (-1)*k1*k7*x11*k24)/k1,
  diff(x12, t) = 0,
  diff(x13, t) = 0,
  diff(x14, t) = 0,
  diff(x15, t) = 0,
  diff(x16, t) = ((-1)*k1*(k5*x16*k21 - k6*x17*k22) + 1*k1*k7*x11*k24 + (-1)*k1*k12*x16/(k11 + x16))/k1,
  diff(x17, t) = (1*k1*(k5*x16*k21 - k6*x17*k22) + (-1)*k1*k8*x17 + (-1)*k1*k14*x17/(k13 + x17))/k1 }