{ diff(x1, t) = ((-1)*k52*k6*x1 + (-1)*k52*k8*x1 + 1*k52*3/500000*k51*x14*x3 + 1*k52*3/125000*k51*x14*x2)/k52,
  diff(x2, t) = (1*k52*k6*x1 + (-1)*k52*k9*x2 + 1*k52*3/250000*k51*x14*x4 + (-1)*k52*3/125000*k51*x14*x2 + 1*k52*3/25000000*k51*x14*x5)/k52,
  diff(x3, t) = (1*k52*k8*x1 + (-1)*k52*k7*x3 + (-1)*k52*3/500000*k51*x14*x3 + 1*k52*3/250000*k51*x14*x4)/k52,
  diff(x4, t) = (1*k52*k7*x3 + (-1)*k52*3/250000*k51*x14*x4 + (-1)*k52*3/250000*k51*x14*x4)/k52,
  diff(x5, t) = (1*k52*k9*x2 + (-1)*k52*3/25000000*k51*x14*x5)/k52,
  diff(x8, t) = ((-1)*k52*k11*x7*x8/(k12 + x8)*(1 - (x9 + x10)/(x6*k10)) + 1*k52*k15*x11)/k52,
  diff(x9, t) = (1*k52*k11*x7*x8/(k12 + x8)*(1 - (x9 + x10)/(x6*k10)) + (-1)*k52*k11*x7*x9/(k12 + x9) + (-1)*k52*(k13*x9*x12 - k14*x11) + 1*k52*k15*x13)/k52,
  diff(x10, t) = (1*k52*k11*x7*x9/(k12 + x9) + (-1)*k52*(k13*x10*x12 - k14*x13))/k52,
  diff(x11, t) = (1*k52*(k13*x9*x12 - k14*x11) + 1*k52*k16*x13 + (-1)*k52*k15*x11)/k52,
  diff(x12, t) = ((-1)*k52*(k13*x9*x12 - k14*x11) + (-1)*k52*(k13*x10*x12 - k14*x13) + 1*k52*k15*x13 + 1*k52*k15*x11)/k52,
  diff(x13, t) = (1*k52*(k13*x10*x12 - k14*x13) + (-1)*k52*k15*x13 + (-1)*k52*k16*x13)/k52,
  diff(x15, t) = (1*k52*k17*x7*x16/(k18 + x16) + (-1)*k52*k19*x14*x15/(k20 + x15) + (-1)*k52*(k21*x15*x17 - k22*x18))/k52,
  diff(x16, t) = ((-1)*k52*k17*x7*x16/(k18 + x16) + 1*k52*k19*x14*x15/(k20 + x15))/k52,
  diff(x17, t) = (-1)*k52*(k21*x15*x17 - k22*x18)/k52,
  diff(x18, t) = 1*k52*(k21*x15*x17 - k22*x18)/k52,
  diff(x20, t) = (1*k52*k23*x7*x21/(k24 + x21) + (-1)*k52*k25*x14*x20/(k26 + x20) + (-1)*k52*(k27*x20*x24 - k28*x23))/k52,
  diff(x21, t) = ((-1)*k52*k23*x7*x21/(k24 + x21) + 1*k52*k25*x14*x20/(k26 + x20))/k52,
  diff(x22, t) = (1*k52*k29*x7*x23 + (-1)*k52*k30*x7*x22)/k52,
  diff(x23, t) = (1*k52*(k27*x20*x24 - k28*x23) + (-1)*k52*k29*x7*x23)/k52,
  diff(x24, t) = (-1)*k52*(k27*x20*x24 - k28*x23)/k52,
  diff(x25, t) = 1*k52*k30*x7*x22/k52,
  diff(x28, t) = (1*k52*k31*x25*x29 + (-1)*k52*(k32*x28*x31 - k33*x38))/k52,
  diff(x29, t) = (-1)*k52*k31*x25*x29/k52,
  diff(x31, t) = (-1)*k52*(k32*x28*x31 - k33*x38)/k52,
  diff(x32, t) = (-1)*k52*(k34*x38*x32 - k35*x33)/k52,
  diff(x33, t) = (1*k52*(k34*x38*x32 - k35*x33) + (-1)*k52*k36*x25*x33)/k52,
  diff(x37, t) = 1*k52*k36*x25*x33/k52,
  diff(x38, t) = (1*k52*(k32*x28*x31 - k33*x38) + (-1)*k52*(k34*x38*x32 - k35*x33))/k52 }