{ diff(x1, t) = 0,
  diff(x2, t) = 0,
  diff(x3, t) = 0,
  diff(x4, t) = 0,
  diff(x5, t) = 1*k1*(k2*k92*k91*k89*k95*k93*k90*k94 - k3*x5)/k1,
  diff(x6, t) = 0,
  diff(x7, t) = 0,
  diff(x8, t) = 0,
  diff(x9, t) = (-1)*k1*(k4*x5*x9*k96 - k5*x11)/k1,
  diff(x10, t) = 0,
  diff(x11, t) = (1*k1*(k4*x5*x9*k96 - k5*x11) + (-1)*k1*(k6*x11*k97 - k7*x13))/k1,
  diff(x12, t) = 0,
  diff(x13, t) = 1*k1*(k6*x11*k97 - k7*x13)/k1,
  diff(x14, t) = ((-1)*k1*k8*x13*x14/(k9 + x14) + 1*k1*k10*x15 + (-1)*k1*k23*x28*x14/(k24 + x14))/k1,
  diff(x15, t) = (1*k1*k8*x13*x14/(k9 + x14) + (-1)*k1*k10*x15 + 1*k1*k23*x28*x14/(k24 + x14) + (-1)*k1*(k27*x15*x30 - k28*x32) + (-1)*k1*(k31*x15*x34 - k32*x35) + (-1)*k1*(k35*x15*x37 - k36*x38) + 1*k1*k37*x32 + 1*k1*k38*x35 + 1*k1*k39*x38 + (-1)*k1*(k58*x15*x29 - k59*x49) + (-1)*k1*(k64*x15*x33 - k65*x50) + (-1)*k1*(k70*x15*x36 - k71*x51) + 1*k1*k85*x49 + 1*k1*k86*x50 + 1*k1*k87*x51)/k1,
  diff(x16, t) = 0,
  diff(x17, t) = 0,
  diff(x18, t) = 0,
  diff(x19, t) = ((-1)*k1*k13*x52*x19/(k14 + x19) + 1*k1*k15*x20)/k1,
  diff(x20, t) = (1*k1*k13*x52*x19/(k14 + x19) + (-1)*k1*k15*x20 + (-1)*k1*(k16*x20 - k17*x21))/k1,
  diff(x21, t) = 1*k1*(k16*x20 - k17*x21)/k1,
  diff(x22, t) = 0,
  diff(x23, t) = 0,
  diff(x24, t) = (1*k1*k18*x21^2 + 1*k1*k19*k101 + (-1)*k1*k20*x24 + (-1)*k1*(k21*x24*k103*k105*k104*k98 - k22*x28) + 1*k1*k88*x39^2)/k1,
  diff(x25, t) = 0,
  diff(x26, t) = 0,
  diff(x27, t) = 0,
  diff(x28, t) = 1*k1*(k21*x24*k103*k105*k104*k98 - k22*x28)/k1,
  diff(x29, t) = ((-1)*k1*(k25*x29*x31 - k26*x30) + (-1)*k1*(k58*x15*x29 - k59*x49) + 1*k1*k60*x46 + (-1)*k1*k61*x29 + (-1)*k1*(k62*x29 - k63*x40))/k1,
  diff(x30, t) = (1*k1*(k25*x29*x31 - k26*x30) + (-1)*k1*(k27*x15*x30 - k28*x32) + (-1)*k1*k40*x30 + 1*k1*k78*x41)/k1,
  diff(x31, t) = ((-1)*k1*(k25*x29*x31 - k26*x30) + (-1)*k1*(k29*x33*x31 - k30*x34) + (-1)*k1*(k33*x36*x31 - k34*x37) + 1*k1*k37*x32 + 1*k1*k38*x35 + 1*k1*k39*x38 + 1*k1*k40*x30 + 1*k1*k41*x34 + 1*k1*k42*x37 + (-1)*k1*(k43*x31 - k44*x39) + (-1)*k1*(k76*x49*x31 - k77*x32) + (-1)*k1*(k79*x50*x31 - k80*x35) + (-1)*k1*(k82*x51*x31 - k83*x38))/k1,
  diff(x32, t) = (1*k1*(k27*x15*x30 - k28*x32) + (-1)*k1*k37*x32 + 1*k1*(k76*x49*x31 - k77*x32))/k1,
  diff(x33, t) = ((-1)*k1*(k29*x33*x31 - k30*x34) + (-1)*k1*(k64*x15*x33 - k65*x50) + 1*k1*k66*x47 + (-1)*k1*k67*x33 + (-1)*k1*(k68*x33 - k69*x42))/k1,
  diff(x34, t) = (1*k1*(k29*x33*x31 - k30*x34) + (-1)*k1*(k31*x15*x34 - k32*x35) + (-1)*k1*k41*x34 + 1*k1*k81*x43)/k1,
  diff(x35, t) = (1*k1*(k31*x15*x34 - k32*x35) + (-1)*k1*k38*x35 + 1*k1*(k79*x50*x31 - k80*x35))/k1,
  diff(x36, t) = ((-1)*k1*(k33*x36*x31 - k34*x37) + (-1)*k1*(k70*x15*x36 - k71*x51) + 1*k1*k72*x48 + (-1)*k1*k73*x36 + (-1)*k1*(k74*x36 - k75*x44))/k1,
  diff(x37, t) = (1*k1*(k33*x36*x31 - k34*x37) + (-1)*k1*(k35*x15*x37 - k36*x38) + (-1)*k1*k42*x37 + 1*k1*k84*x45)/k1,
  diff(x38, t) = (1*k1*(k35*x15*x37 - k36*x38) + (-1)*k1*k39*x38 + 1*k1*(k82*x51*x31 - k83*x38))/k1,
  diff(x39, t) = (1*k1*(k43*x31 - k44*x39) + (-1)*k1*(k45*x40*x39 - k46*x41) + (-1)*k1*(k47*x42*x39 - k48*x43) + (-1)*k1*(k49*x44*x39 - k50*x45))/k1,
  diff(x40, t) = ((-1)*k1*(k45*x40*x39 - k46*x41) + 1*k1*(k62*x29 - k63*x40))/k1,
  diff(x41, t) = (1*k1*(k45*x40*x39 - k46*x41) + (-1)*k1*k78*x41)/k1,
  diff(x42, t) = ((-1)*k1*(k47*x42*x39 - k48*x43) + 1*k1*(k68*x33 - k69*x42))/k1,
  diff(x43, t) = (1*k1*(k47*x42*x39 - k48*x43) + (-1)*k1*k81*x43)/k1,
  diff(x44, t) = ((-1)*k1*(k49*x44*x39 - k50*x45) + 1*k1*(k74*x36 - k75*x44))/k1,
  diff(x45, t) = (1*k1*(k49*x44*x39 - k50*x45) + (-1)*k1*k84*x45)/k1,
  diff(x46, t) = (1*k1*k51*k101 + 1*k1*k52*x39^2 + (-1)*k1*k53*x46)/k1,
  diff(x47, t) = (1*k1*k54*k101 + (-1)*k1*k55*x47)/k1,
  diff(x48, t) = (1*k1*k56*k101 + (-1)*k1*k57*x48)/k1,
  diff(x49, t) = (1*k1*(k58*x15*x29 - k59*x49) + (-1)*k1*(k76*x49*x31 - k77*x32) + (-1)*k1*k85*x49)/k1,
  diff(x50, t) = (1*k1*(k64*x15*x33 - k65*x50) + (-1)*k1*(k79*x50*x31 - k80*x35) + (-1)*k1*k86*x50)/k1,
  diff(x51, t) = (1*k1*(k70*x15*x36 - k71*x51) + (-1)*k1*(k82*x51*x31 - k83*x38) + (-1)*k1*k87*x51)/k1,
  diff(x52, t) = 1*k1*(k11*k92*k91*k89*k95*k100*k99*k98*k106 - k12*x52)/k1,
  diff(x53, t) = 0 }