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