{ diff(x1, t) = (1*k15 + (-1)*(k13*(k5 - x9) + k14*x9)*x1 + (-1)*k16*x1*x10)/k37,
  diff(x2, t) = (1*k16*x1*x10 + 1*k17*x3 + (-1)*(k10*x7 + k11*(k3 - x7))*x2 + (-1)*k18*x2 + (-1)*(k13*(k5 - x9) + k14*x9)*x2 + 1*(k7*(k2 - x6) + k8*x6)*x4)/k37,
  diff(x3, t) = ((-1)*k17*x3 + 1*k18*x2 + (-1)*(k10*x7 + k11*(k3 - x7))*x3 + 1*(k7*(k2 - x6) + k8*x6)*x5 + (-1)*(k13*(k5 - x9) + k14*x9)*x3)/k37,
  diff(x4, t) = (1*(k10*x7 + k11*(k3 - x7))*x2 + (-1)*(k7*(k2 - x6) + k8*x6)*x4 + (-1)*k18*x4 + (-1)*(k13*(k5 - x9) + k14*x9)*x4 + 1*k17*x5)/k37,
  diff(x5, t) = (1*k18*x4 + (-1)*k17*x5 + 1*(k10*x7 + k11*(k3 - x7))*x3 + (-1)*(k7*(k2 - x6) + k8*x6)*x5 + (-1)*(k13*(k5 - x9) + k14*x9)*x5)/k37,
  diff(x6, t) = (1*k19*x3*(k2 - x6)/(k20 + k2 - x6) + (-1)*k21*x6/(k22 + x6))/k37,
  diff(x7, t) = (1*k23*x3*(k3 - x7)/(k24 + k3 - x7) + (-1)*k25*x7/(k26 + x7))/k37,
  diff(x8, t) = (1*k27*x3*(k4 - x8)/(k28 + k4 - x8) + (-1)*k29*x8/(k30 + x8))/k37,
  diff(x9, t) = (1*k31*x8*(k5 - x9)/(k32 + k5 - x9) + (-1)*k33*x9/(k34 + x9))/k37 }