{ diff(x1, t) = (1*k29*(k1*x1*(1 - x1/k4) + k10*x1*x5/(k13 + x1)) + (-1)*k29*k7*x4*x1)/k29,
  diff(x2, t) = (1*k29*(k2*x1*x1/k4*(1 - x2/k5) + k11*x2*x5/(k14 + x2)) + (-1)*k29*(k16*x2 + k8*x4*x2))/k29,
  diff(x3, t) = (1*k29*k3*x3*(1 - x3/k6) + (-1)*k29*(k18*x3*x2 + k12*x3*x5/(k15 + x3)))/k29,
  diff(x4, t) = (1*k29*(k19 + k20*x4*x2/(k21 + x2)) + (-1)*k29*(k9*x4*x2 + k17*x4 + k22*x4*x5/(k23 + x5)))/k29,
  diff(x5, t) = (1*k29*k24 + (-1)*k29*(k28 + k25*x1/(k13 + x1) + k26*x2/(k14 + x2) + k27*x3/(k15 + x3))*x5)/k29 }