{ diff(x1, t) = (-1)*k1*x1*x29*k3/k1,
  diff(x2, t) = ((-1)*k2*x2*x1*k4 + 1*k2*x3*k5 + (-1)*k2*x2*k7)/k2,
  diff(x3, t) = (1*k2*x2*x1*k4 + (-1)*k2*x3*k5 + (-1)*k2*x3*k6)/k2,
  diff(x4, t) = ((-1)*k2*x3*x4*k8 + 1*k2*x7*x6*k11)/k2,
  diff(x5, t) = (1*k2*x3*x4*k8 + (-1)*k2*x5*k9 + (-1)*k2*x5*x36*k10)/k2,
  diff(x6, t) = (1*k2*x3*x4*k8 + (-1)*k2*x7*x6*k11 + (-1)*k2*x6*x8*k12 + 1*k2*x9*k13 + 1*k2*x17*k23 + 1*k2*x18*k25 + 1*k2*x19*k27 + 1*k2*x20*k29 + 1*k2*x22*k34 + (-1)*k2*x6*x31*k44 + 1*k2*x33*k45)/k2,
  diff(x7, t) = (1*k2*x5*k9 + 1*k2*x5*x36*k10 + (-1)*k2*x7*x6*k11)/k2,
  diff(x8, t) = ((-1)*k2*x6*x8*k12 + 1*k2*x9*k13 + 1*k2*x12*x15*k18 + (-1)*k2*x8*k19)/k2,
  diff(x9, t) = (1*k2*x6*x8*k12 + (-1)*k2*x9*k13 + (-1)*k2*x9*x16*k20 + 1*k2*x17*k21)/k2,
  diff(x10, t) = ((-1)*k2*x10*x11*k14 + 1*k2*x12*k15 + 1*k2*x8*k19 + 1*k2*x17*k23 + 1*k2*x18*k25 + 1*k2*x19*k27 + 1*k2*x20*k29 + 1*k2*x22*k34)/k2,
  diff(x11, t) = ((-1)*k2*x10*x11*k14 + 1*k2*x12*k15 + 1*k2*x8*k19 + 1*k2*x17*k23 + 1*k2*x18*k25 + 1*k2*x19*k27 + 1*k2*x20*k29 + 1*k2*x22*k34)/k2,
  diff(x12, t) = (1*k2*x10*x11*k14 + (-1)*k2*x12*k15 + (-1)*k2*x12*x15*k18)/k2,
  diff(x13, t) = ((-1)*k2*x13*x14*k16 + 1*k2*x15*k17 + 1*k2*x8*k19 + 1*k2*x17*k23 + 1*k2*x18*k25 + 1*k2*x19*k27 + 1*k2*x20*k29 + 1*k2*x22*k34)/k2,
  diff(x14, t) = ((-1)*k2*x13*x14*k16 + 1*k2*x15*k17 + 1*k2*x8*k19 + 1*k2*x17*k23 + 1*k2*x18*k25 + 1*k2*x19*k27 + 1*k2*x20*k29 + (-1)*k2*x22*x14*k31 + 1*k2*x24*k32 + 1*k2*x23*k35)/k2,
  diff(x15, t) = (1*k2*x13*x14*k16 + (-1)*k2*x15*k17 + (-1)*k2*x12*x15*k18)/k2,
  diff(x16, t) = ((-1)*k2*x9*x16*k20 + 1*k2*x17*k21 + 1*k2*x17*k23 + 1*k2*x18*k25 + 1*k2*x19*k27 + 1*k2*x20*k29 + 1*k2*x22*k34)/k2,
  diff(x17, t) = (1*k2*x9*x16*k20 + (-1)*k2*x17*k21 + (-1)*k2*x17*k22 + (-1)*k2*x17*k23)/k2,
  diff(x18, t) = (1*k2*x17*k22 + (-1)*k2*x18*k24 + (-1)*k2*x18*k25)/k2,
  diff(x19, t) = (1*k2*x18*k24 + (-1)*k2*x19*k26 + (-1)*k2*x19*k27)/k2,
  diff(x20, t) = (1*k2*x19*k26 + (-1)*k2*x20*k28 + (-1)*k2*x20*k29)/k2,
  diff(x21, t) = (1*k2*x20*k28 + (-1)*k2*x21*k30 + 1*k2*x24*k33)/k2,
  diff(x22, t) = (1*k2*x21*k30 + (-1)*k2*x22*x14*k31 + 1*k2*x24*k32 + (-1)*k2*x22*k34)/k2,
  diff(x23, t) = (1*k2*x21*k30 + (-1)*k2*x23*k35 + (-1)*k2*x25*x23*k36 + 1*k2*x26*k37)/k2,
  diff(x24, t) = (1*k2*x22*x14*k31 + (-1)*k2*x24*k32 + (-1)*k2*x24*k33)/k2,
  diff(x25, t) = ((-1)*k2*x25*x23*k36 + 1*k2*x26*k37)/k2,
  diff(x26, t) = (1*k2*x25*x23*k36 + (-1)*k2*x26*k37)/k2,
  diff(x27, t) = ((-1)*k2*x26*x27*k38 + 1*k2*x28*k39)/k2,
  diff(x28, t) = (1*k2*x26*x27*k38 + (-1)*k2*x28*k39 + (-1)*k2*x28*k40)/k2,
  diff(x29, t) = 1*k2*x28*k40/k1,
  diff(x30, t) = ((-1)*k2*x30*x23*x23/(100*100 + x23*x23)*k41 + 1*k2*x31*k42 + (-1)*k2*x30*x35*k43)/k2,
  diff(x31, t) = (1*k2*x30*x23*x23/(100*100 + x23*x23)*k41 + (-1)*k2*x31*k42 + (-1)*k2*x6*x31*k44 + 1*k2*x33*k45 + 1*k2*x34*k46 + (-1)*k2*x31*x35*k47)/k2,
  diff(x32, t) = 1*k2*x30*x35*k43/k2,
  diff(x33, t) = (1*k2*x6*x31*k44 + (-1)*k2*x33*k45)/k2,
  diff(x34, t) = ((-1)*k2*x34*k46 + 1*k2*x31*x35*k47)/k2,
  diff(x35, t) = (1*k2*x34*k46 + (-1)*k2*x31*x35*k47)/k2,
  diff(x36, t) = (1*k2*x23^2/(4^2 + x23^2)*k48 + (-1)*k2*x36*k49)/k2,
  diff(x37, t) = 0 }