{ diff(x1, t) = (1*k58*k1 + (-1)*k58*k25*x1)/k58,
  diff(x2, t) = (1*k58*k2 + 1*k58*k11*x1 + 1*k58*k18*x8^k19/(k20^k19 + x8^k19) + (-1)*k58*k26*x2)/k58,
  diff(x3, t) = (1*k58*(k3 + k12*x2)*k47*(1 - k51)/(k48*x6) + (-1)*k58*k27*x3)/k58,
  diff(x4, t) = (1*k58*k44*x9 + (-1)*k58*k28*x4 + (-1)*k58*k37*x4*x5 + 1*k58*k13*x6 + (-1)*k58*k42*k10*x4*(k21*x6 + k7)*(1 - k9*k50)/(k23 + x4) + (-1)*k58*k14*k10*k46*x4*(1 - k9*k50)/(k22 + x4) + 1*k58*k15*x10)/k58,
  diff(x5, t) = ((-1)*k58*k37*x4*x5 + 1*k58*k13*x6 + 1*k58*k24 + (-1)*k58*k29*x5)/k58,
  diff(x6, t) = (1*k58*k37*x4*x5 + (-1)*k58*k13*x6 + (-1)*k58*k30*x6)/k58,
  diff(x7, t) = (1*k58*k4 + (-1)*k58*k31*x7 + 1*k58*k39*x1 + 1*k58*k17*x3 + 1*k58*k43*x6)/k58,
  diff(x8, t) = ((-1)*k58*k40*x8 + 1*k58*(k5 + k40*x3 + k38*(1 - k35*k50)*x1)*(1 - k34*k52))/k58,
  diff(x9, t) = (1*k58*k6 + (-1)*k58*k32*x9 + 1*k58*k41*x7 + 1*k58*k36*x8)/k58,
  diff(x10, t) = (1*k58*k14*k10*k46*x4*(1 - k9*k50)/(k22 + x4) + (-1)*k58*k15*x10 + (-1)*k58*k16*k10*x10*(k21*x6 + k7)/(k45 + x10) + (-1)*k58*k33*x10)/k58 }