{ diff(x1, t) = ((-1)*k57*(k10*k1 + k11*x10 + k12*x11)*x1 + 1*k57*k13*x2 + 1*k57*k9*x4 + (-1)*k57*(k7*x1*x3 - k8*x4))/k57,
  diff(x2, t) = (1*k57*(k10*k1 + k11*x10 + k12*x11)*x1 + 1*k57*(k10*k1 + k11*x10 + k12*x11)*x4 + (-1)*k57*k13*x2)/k57,
  diff(x3, t) = (1*k57*(k10*k1 + k11*x10 + k12*x11)*x4 + 1*k57*(k4 + k5*x3/(k6 + x3)) + (-1)*k57*k9*x3 + (-1)*k57*(k7*x1*x3 - k8*x4))/k57,
  diff(x4, t) = ((-1)*k57*(k10*k1 + k11*x10 + k12*x11)*x4 + (-1)*k57*k9*x4 + 1*k57*(k7*x1*x3 - k8*x4))/k57,
  diff(x5, t) = (1*k57*k14*x3 + (-1)*k57*k16*x5 + (-1)*k57*(k21*x5*x6 - k22*x8) + 1*k57*(k17 + k18*x10 + k19*x11)*x8)/k57,
  diff(x6, t) = ((-1)*k57*(k21*x5*x6 - k22*x8) + 1*k57*k15*x8 + (-1)*k57*(k17 + k18*x10 + k19*x11)*x6 + 1*k57*k20*x7)/k57,
  diff(x7, t) = (1*k57*(k17 + k18*x10 + k19*x11)*x6 + 1*k57*(k17 + k18*x10 + k19*x11)*x8 + (-1)*k57*k20*x7)/k57,
  diff(x8, t) = (1*k57*(k21*x5*x6 - k22*x8) + (-1)*k57*k15*x8 + (-1)*k57*(k17 + k18*x10 + k19*x11)*x8)/k57,
  diff(x9, t) = (1*k57*(k23 + k24*x21) + (-1)*k57*(k30*x10*x9 - k31*x12) + (-1)*k57*(k30*x11*x9 - k31*x13) + 1*k57*(k32 + k34*x11)*x12 + 1*k57*(k33 + k35*x6)*x13 + (-1)*k57*(k25 + k26*k2*(x10 + x11) + k27*k3*x18)*x9 + 1*k57*k37*x15 + (-1)*k57*(k44*x14*x9 - k45*x15) + (-1)*k57*(k44*x18*x9 - k45*x19))/k57,
  diff(x10, t) = (1*k57*k28*x3 + (-1)*k57*(k30*x10*x9 - k31*x12) + 1*k57*(k25 + k26*k2*(x10 + x11) + k27*k3*x18)*x12 + (-1)*k57*(k32 + k34*x11)*x10)/k57,
  diff(x11, t) = (1*k57*k29*x3 + (-1)*k57*(k30*x11*x9 - k31*x13) + 1*k57*(k25 + k26*k2*(x10 + x11) + k27*k3*x18)*x13 + (-1)*k57*(k33 + k35*x6)*x11)/k57,
  diff(x12, t) = (1*k57*(k30*x10*x9 - k31*x12) + (-1)*k57*(k25 + k26*k2*(x10 + x11) + k27*k3*x18)*x12 + (-1)*k57*(k32 + k34*x11)*x12)/k57,
  diff(x13, t) = (1*k57*(k30*x11*x9 - k31*x13) + (-1)*k57*(k25 + k26*k2*(x10 + x11) + k27*k3*x18)*x13 + (-1)*k57*(k33 + k35*x6)*x13)/k57,
  diff(x14, t) = (1*k57*k36 + (-1)*k57*k37*x14 + (-1)*k57*(k44*x14*x9 - k45*x15) + 1*k57*(k25 + k26*k2*(x10 + x11) + k27*k3*x18)*x15 + (-1)*k57*(k42*x14*x17 - k43*x18) + 1*k57*piecewise(0, x20 < 1, piecewise(1*x18, x20 > 1, 1/2*x18)))/k57,
  diff(x15, t) = ((-1)*k57*k37*x15 + 1*k57*(k44*x14*x9 - k45*x15) + (-1)*k57*(k25 + k26*k2*(x10 + x11) + k27*k3*x18)*x15 + (-1)*k57*(k42*x15*x17 - k43*x19) + 1*k57*piecewise(0, x20 < 1, piecewise(1*x19, x20 > 1, 1/2*x19)))/k57,
  diff(x16, t) = ((-1)*k57*k38*(x10 + x11)^k41/(k40^k41 + (x10 + x11)^k41)*x16 + 1*k57*k39*x17 + (-1)*k57*piecewise(0, x20 < 1, piecewise(1*x16, x20 > 1, 1/2*x16)))/k57,
  diff(x17, t) = ((-1)*k57*(k42*x14*x17 - k43*x18) + (-1)*k57*(k42*x15*x17 - k43*x19) + 1*k57*k38*(x10 + x11)^k41/(k40^k41 + (x10 + x11)^k41)*x16 + (-1)*k57*k39*x17 + (-1)*k57*piecewise(0, x20 < 1, piecewise(1*x17, x20 > 1, 1/2*x17)))/k57,
  diff(x18, t) = (1*k57*(k42*x14*x17 - k43*x18) + (-1)*k57*(k44*x18*x9 - k45*x19) + 1*k57*(k25 + k26*k2*(x10 + x11) + k27*k3*x18)*x19 + (-1)*k57*piecewise(0, x20 < 1, piecewise(1*x18, x20 > 1, 1/2*x18)))/k57,
  diff(x19, t) = (1*k57*(k42*x15*x17 - k43*x19) + 1*k57*(k44*x18*x9 - k45*x19) + (-1)*k57*(k25 + k26*k2*(x10 + x11) + k27*k3*x18)*x19 + (-1)*k57*piecewise(0, x20 < 1, piecewise(1*x19, x20 > 1, 1/2*x19)))/k57,
  diff(x20, t) = 1*k57*k46*x18/k57,
  diff(x21, t) = (1*k57*k47 + (-1)*k57*k48/(k49 + x22)*x21)/k57,
  diff(x22, t) = (1*k57*k50 + 1*k57*k51*x18 + (-1)*k57*(k52 + k53*x21/(k54 + x22))*x22)/k57,
  diff(x23, t) = (1*k57*k55 + (-1)*k57*(k56 + k35*x6)*x23)/k57 }