{ diff(x1, t) = ((-1)*k1*x1*x33*k40 + 1*k7*x10*x13*k40)/k40,
  diff(x2, t) = ((-1)*k2*x2*x10*k40 + 1*k6*x11*x13*k40)/k40,
  diff(x3, t) = ((-1)*k3*x3*x11*k40 + 1*k5*x13*k40)/k40,
  diff(x4, t) = ((-1)*k8*x4*x11*k40 + 1*k9*x14*k40 + 1*k25*x20*k40)/k40,
  diff(x5, t) = ((-1)*k10*x5*x14*k40 + 1*k11*x15*k40)/k40,
  diff(x6, t) = ((-1)*k12*x6*x15*k40 + 1*k17*x21*k40)/k40,
  diff(x7, t) = ((-1)*k14*x7*x15*k40 + 1*k17*x22*k40)/k40,
  diff(x8, t) = ((-1)*k18*x8*x16*k40 + (-1)*k18*x8*x17*k40 + 1*k23*x23*k40)/k40,
  diff(x9, t) = ((-1)*k20*x9*x16*k40 + (-1)*k20*x9*x17*k40 + 1*k23*x24*k40)/k40,
  diff(x10, t) = (1*k1*x1*x33*k40 + (-1)*k7*x10*x13*k40)/k40,
  diff(x11, t) = (1*k2*x2*x10*k40 + (-1)*k6*x11*x13*k40)/k40,
  diff(x12, t) = (1*k3*x3*x11*k40 + (-1)*k4*x12*k40)/k40,
  diff(x13, t) = (1*k4*x32*k40 + (-1)*k5*x13*k40)/k40,
  diff(x14, t) = (1*k8*x4*x11*k40 + (-1)*k9*x14*k40 + (-1)*k24*x14*x18*k40 + (-1)*k24*x14*x19*k40)/k40,
  diff(x15, t) = (1*k10*x5*x14*k40 + (-1)*k11*x15*k40)/k40,
  diff(x16, t) = (1*k13*x21*x15*k40 + (-1)*k16*x16*k40)/k40,
  diff(x17, t) = (1*k15*x22*x15*k40 + (-1)*k16*x17*k40)/k40,
  diff(x18, t) = (1*k19*x23*x16*k40 + 1*k19*x23*x17*k40 + (-1)*k22*x18*k40)/k40,
  diff(x19, t) = (1*k21*x24*x16*k40 + 1*k21*x24*x17*k40 + (-1)*k22*x19*k40)/k40,
  diff(x20, t) = (1*k24*x14*x18*k40 + 1*k24*x14*x19*k40 + (-1)*k25*x20*k40)/k40,
  diff(x21, t) = (1*k12*x6*x15*k40 + (-1)*k13*x21*x15*k40 + 1*k16*x16*k40 + (-1)*k17*x21*k40)/k40,
  diff(x22, t) = (1*k14*x7*x15*k40 + (-1)*k15*x22*x15*k40 + 1*k16*x17*k40 + (-1)*k17*x22*k40)/k40,
  diff(x23, t) = (1*k18*x8*x16*k40 + 1*k18*x8*x17*k40 + (-1)*k19*x23*x16*k40 + (-1)*k19*x23*x17*k40 + 1*k22*x18*k40 + (-1)*k23*x23*k40)/k40,
  diff(x24, t) = (1*k20*x9*x16*k40 + 1*k20*x9*x17*k40 + (-1)*k21*x24*x16*k40 + (-1)*k21*x24*x17*k40 + 1*k22*x19*k40 + (-1)*k23*x24*k40)/k40,
  diff(x25, t) = (1*k4*x12*k40 + (-1)*k4*x25*k40)/k40,
  diff(x26, t) = (1*k4*x25*k40 + (-1)*k4*x26*k40)/k40,
  diff(x27, t) = (1*k4*x26*k40 + (-1)*k4*x27*k40)/k40,
  diff(x28, t) = (1*k4*x27*k40 + (-1)*k4*x28*k40)/k40,
  diff(x29, t) = (1*k4*x28*k40 + (-1)*k4*x29*k40)/k40,
  diff(x30, t) = (1*k4*x29*k40 + (-1)*k4*x30*k40)/k40,
  diff(x31, t) = (1*k4*x30*k40 + (-1)*k4*x31*k40)/k40,
  diff(x32, t) = (1*k4*x31*k40 + (-1)*k4*x32*k40)/k40,
  diff(x33, t) = 0/k40 }