[ Eq(Derivative(x1, t), ((-1)*k3*x4*x1 + 1*k4*x5 + 1*k5*x5*x14 + (-1)*k8*x1*x3 + 1*k9*x2 + 1*k16*x11 + (-1)*k17*x1*x14)/k22),
  Eq(Derivative(x2, t), (1*k8*x1*x3 + (-1)*k9*x2)/k22),
  Eq(Derivative(x3, t), ((-1)*k8*x1*x3 + 1*k9*x2 + (-2)*(x3 - 1)*k10*x3/2 + (-1)*k11*x3*x7 + 1*k12*x6 + 2*k13*x7)/k22),
  Eq(Derivative(x4, t), (1*k2*x8*x13 + (-1)*k3*x4*x1 + 1*k4*x5 + (-1)*k6*x4*x14 + (-2)*(x4 - 1)*k7*x4/2 + (-1)*k7*x4*x9)/k22),
  Eq(Derivative(x5, t), (1*k3*x4*x1 + (-1)*k4*x5 + (-1)*k5*x5*x14)/k22),
  Eq(Derivative(x6, t), (1*k11*x3*x7 + (-1)*k12*x6 + (-1)*k14*x10*x6 + 1*k15*x11)/k22),
  Eq(Derivative(x7, t), (1*(x3 - 1)*k10*x3/2 + (-1)*k11*x3*x7 + 1*k12*x6 + (-1)*k13*x7)/k22),
  Eq(Derivative(x8, t), (1*k1 + (-1)*k2*x8*x13 + 1*k5*x5*x14)/k22),
  Eq(Derivative(x9, t), (1*(x4 - 1)*k7*x4/2 + 1*k7*x4*x9)/k22),
  Eq(Derivative(x10, t), ((-1)*k14*x10*x6 + 1*k15*x11)/k22),
  Eq(Derivative(x11, t), (1*k14*x10*x6 + (-1)*k15*x11)/k22),
  Eq(Derivative(x12, t), 1*1/k22),
  Eq(Derivative(x13, t), (1*k20 + (-1)*k21*x13)/k22),
  Eq(Derivative(x14, t), ((-1)*k5*x5*x14 + (-1)*k6*x4*x14 + (-1)*k17*x1*x14 + 1*k18*x15 + (-1)*k19*x14)/k22),
  Eq(Derivative(x15, t), (1*k5*x5*x14 + 1*k6*x4*x14 + 1*k17*x1*x14 + (-1)*k18*x15 + 1*k19*x14)/k22),
  Eq(Derivative(x16, t), 0) ]