{ diff(x1, t) = 1*k43*x25/k2,
  diff(x2, t) = (-1)*k44*x25/(x25 + k45)/k2,
  diff(x3, t) = (-1)*k46/k2,
  diff(x4, t) = 1*k46/k2,
  diff(x5, t) = ((-1)*x7*x5*k56/(k55 + x5) + 1*k58*x26/(x26 + k57))/k3,
  diff(x6, t) = ((-1)*k47*x6/(x6 + k48) + 1*k49*x19*x20/(k50 + x19))/k2,
  diff(x7, t) = ((-1)*k51*x7/(k52 + x7) + 1*k53*x23*x6/(k54 + x23))/k2,
  diff(x8, t) = 1*exp(k59*(k62 - x4/k60))/(1 + 2*exp(k59*(k62 - x4/k60)) + exp(2*k59*(k62 - x4/k60)))/k61*3657/125/k2,
  diff(x9, t) = 1*k44*x25/(x25 + k45)/k2,
  diff(x10, t) = ((-1)*(k13*x10*x12 - k14*x28) + 1*k21*x17/(k22 + x17))/k2,
  diff(x11, t) = 1*exp(k11 - (x4/k10)^(13/10))*k9*(x4/k12)^(3/10)*(k11 - (x4/k10)^(13/10))/k2,
  diff(x12, t) = (1*exp(k6 - (x4/k5)^(7/20))*k4*((x4 + k8)/k7)^(-13/20)*(k6 - (x4/k5)^(7/20)) + (-1)*(k13*x10*x12 - k14*x28) + 1*(k19*x31 - k20*x12*x32))/k2,
  diff(x13, t) = (-1)*exp(k11 - (x4/k10)^(13/10))*k9*(x4/k12)^(3/10)*(k11 - (x4/k10)^(13/10))/k2,
  diff(x14, t) = (-1)*exp(k6 - (x4/k5)^(7/20))*k4*((x4 + k8)/k7)^(-13/20)*(k6 - (x4/k5)^(7/20))/k2,
  diff(x15, t) = ((-1)*k35*x20*x15/(k36 + x15) + 1*k38*x19/(x19 + k37))/k2,
  diff(x16, t) = ((-1)*k31*x21*x16/(k32 + x16) + 1*k34*x20/(k33 + x20))/k2,
  diff(x17, t) = ((-1)*k21*x17/(k22 + x17) + 1*k23*x32)/k2,
  diff(x18, t) = ((-1)*k39*x18*x6/(x18 + k40) + 1*k42*x23/(x23 + k41))/k2,
  diff(x19, t) = (1*k35*x20*x15/(k36 + x15) + (-1)*k38*x19/(x19 + k37) + 1*k47*x6/(x6 + k48) + (-1)*k49*x19*x20/(k50 + x19))/k2,
  diff(x20, t) = (1*k31*x21*x16/(k32 + x16) + (-1)*k34*x20/(k33 + x20))/k2,
  diff(x21, t) = (1*k27*x32*x22/(x22 + k28) + (-1)*k29*x21/(k30 + x21))/k2,
  diff(x22, t) = ((-1)*k27*x32*x22/(x22 + k28) + 1*k29*x21/(k30 + x21))/k2,
  diff(x23, t) = (1*k39*x18*x6/(x18 + k40) + (-1)*k42*x23/(x23 + k41) + 1*k51*x7/(k52 + x7) + (-1)*k53*x23*x6/(k54 + x23))/k2,
  diff(x24, t) = (-1)*exp(k59*(k62 - x4/k60))/(1 + 2*exp(k59*(k62 - x4/k60)) + exp(2*k59*(k62 - x4/k60)))/k61*3657/125/k2,
  diff(x25, t) = (1*k24*x33 + (-1)*k43*x25)/k2,
  diff(x26, t) = (1*x7*x5*k56/(k55 + x5) + (-1)*k58*x26/(x26 + k57))/k3,
  diff(x27, t) = (-1)*(k25*x26 + k26*x8)/k3,
  diff(x28, t) = (1*(k13*x10*x12 - k14*x28) + (-1)*(k15*x28*x11 - k16*x30))/k2,
  diff(x29, t) = ((-1)*(k17*x29*x30 - k18*x31) + 1*k23*x32)/k2,
  diff(x30, t) = (1*(k15*x28*x11 - k16*x30) + (-1)*(k17*x29*x30 - k18*x31))/k2,
  diff(x31, t) = (1*(k17*x29*x30 - k18*x31) + (-1)*(k19*x31 - k20*x12*x32))/k2,
  diff(x32, t) = (1*(k19*x31 - k20*x12*x32) + (-1)*k23*x32)/k2,
  diff(x33, t) = ((-1)*k24*x33 + 1*(k25*x26 + k26*x8))/k3,
  diff(x34, t) = 0/k1 }