{ diff(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,
  diff(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,
  diff(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,
  diff(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,
  diff(x5, t) = (1*k41*k1*x14*x15 + (-1)*k41*k2*x5)/k41,
  diff(x6, t) = (1*k41*k9*x7 + (-1)*k41*k24*x6 + 1*k41*k25 + 1*k41*k12*x8 + (-1)*k41*k8*x6*x4)/k41,
  diff(x7, t) = ((-1)*k41*k9*x7 + 1*k41*k8*x6*x4 + (-1)*k41*k39*x7)/k41,
  diff(x8, t) = (1*k41*k22 + (-1)*k41*k23*x8)/k41,
  diff(x9, t) = (1*k41*k34*x8 + 1*k41*k35 + (-1)*k41*k36*x9 + (-1)*k41*k37*x9*x10 + 1*k41*k38*x11)/k41,
  diff(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,
  diff(x11, t) = (1*k41*k37*x9*x10 + (-1)*k41*k38*x11)/k41,
  diff(x12, t) = ((-1)*k41*(k19*(x1 + x2)*x12 + k20*x12) + 1*k41*k21*x10*x13 + 1*k41*k26/(1 + k27*x8))/k41,
  diff(x13, t) = (1*k41*k29 + (-1)*k41*(k16*x10*x13 + k17*x13) + 1*k41*k18*x16)/k41,
  diff(x14, t) = ((-1)*k41*k1*x14*x15 + 1*k41*k2*x5)/k41,
  diff(x15, t) = (1*k41*k10 + (-1)*k41*k11*x15)/k41,
  diff(x16, t) = (1*k41*(k16*x10*x13 + k17*x13) + (-1)*k41*k18*x16 + (-1)*k41*k30*x16)/k41 }