[ Eq(Derivative(x1, t), ((-1)*k41*(k3*x5*x1 + k6*x1) + 1*k41*k4*x2 + (-1)*k41*k15*x1 + 1*k41*(k13*x10*x3 + k14*x3))/k41),
  Eq(Derivative(x2, t), (1*k41*(k3*x5*x1 + k6*x1) + (-1)*k41*k4*x2 + 1*k41*(k31*x10*x4 + k32*x4) + (-1)*k41*k33*x2)/k41),
  Eq(Derivative(x3, t), (1*k41*k15*x1 + (-1)*k41*(k13*x10*x3 + k14*x3) + (-1)*k41*(k5*x5*x3 + k7*x3) + 1*k41*k9*x7 + 1*k41*k40)/k41),
  Eq(Derivative(x4, t), ((-1)*k41*(k31*x10*x4 + k32*x4) + 1*k41*k33*x2 + 1*k41*(k5*x5*x3 + k7*x3) + (-1)*k41*k8*x6*x4)/k41),
  Eq(Derivative(x5, t), (1*k41*k1*x14*x15 + (-1)*k41*k2*x5)/k41),
  Eq(Derivative(x6, t), (1*k41*k9*x7 + (-1)*k41*k24*x6 + 1*k41*k25 + 1*k41*k12*x8 + (-1)*k41*k8*x6*x4)/k41),
  Eq(Derivative(x7, t), ((-1)*k41*k9*x7 + 1*k41*k8*x6*x4 + (-1)*k41*k39*x7)/k41),
  Eq(Derivative(x8, t), (1*k41*k22 + (-1)*k41*k23*x8)/k41),
  Eq(Derivative(x9, t), (1*k41*k34*x8 + 1*k41*k35 + (-1)*k41*k36*x9 + (-1)*k41*k37*x9*x10 + 1*k41*k38*x11)/k41),
  Eq(Derivative(x10, t), ((-1)*k41*k37*x9*x10 + 1*k41*k38*x11 + 1*k41*(k19*(x1 + x2)*x12 + k20*x12) + (-1)*k41*k21*x10*x13 + (-1)*k41*k28*x10*x10)/k41),
  Eq(Derivative(x11, t), (1*k41*k37*x9*x10 + (-1)*k41*k38*x11)/k41),
  Eq(Derivative(x12, t), ((-1)*k41*(k19*(x1 + x2)*x12 + k20*x12) + 1*k41*k21*x10*x13 + 1*k41*k26/(1 + k27*x8))/k41),
  Eq(Derivative(x13, t), (1*k41*k29 + (-1)*k41*(k16*x10*x13 + k17*x13) + 1*k41*k18*x16)/k41),
  Eq(Derivative(x14, t), ((-1)*k41*k1*x14*x15 + 1*k41*k2*x5)/k41),
  Eq(Derivative(x15, t), (1*k41*k10 + (-1)*k41*k11*x15)/k41),
  Eq(Derivative(x16, t), (1*k41*(k16*x10*x13 + k17*x13) + (-1)*k41*k18*x16 + (-1)*k41*k30*x16)/k41) ]