[ Eq(Derivative(x1, t), (1*k139*k1 + (-1)*k139*k2*x1 + (-1)*k139*(k3*x1*x2 - k4*x3))/k139),
  Eq(Derivative(x2, t), ((-1)*k139*(k3*x1*x2 - k4*x3) + 1*k139*k13*x3)/k139),
  Eq(Derivative(x3, t), (1*k139*(k3*x1*x2 - k4*x3) + (-1)*k139*k13*x3 + (-1)*k139*k18*x15*x3 + 1*k139*k19*x16 + (-1)*k139*k20*x3*x17 + 1*k139*k21*x18 + (-1)*k139*k44*x3*x25)/k139),
  Eq(Derivative(x4, t), (1*k139*k5*x5 + (-1)*k139*k6*x4 + (-1)*k139*k7*x4*x6 + 1*k139*k8*x7)/k139),
  Eq(Derivative(x5, t), ((-1)*k139*k45*x5*x26 + 1*k139*k49*x14*x28 + 1*k139*k50*x28*x11 + 1*k139*k51*x5 + 1*k139*k52 + (-1)*k139*k53*x5 + (-1)*k139*k54*x5*x11)/k139),
  Eq(Derivative(x6, t), ((-1)*k139*k7*x4*x6 + 1*k139*k8*x7 + (-1)*k139*k11*x8*x6 + 1*k139*k12*x10 + 1*k139*k14*x11*(x12 + x13) + 1*k139*k15*x10*(x12 + x13) + 1*k139*k16*x7 + 1*k139*k17*x14*x14)/k139),
  Eq(Derivative(x7, t), (1*k139*k7*x4*x6 + (-1)*k139*k8*x7 + (-1)*k139*k16*x7 + (-1)*k139*k22*x19*x7 + 1*k139*k23*x14)/k139),
  Eq(Derivative(x8, t), (1*k139*k9*x9 + (-1)*k139*k10*x8 + (-1)*k139*k11*x8*x6 + 1*k139*k12*x10 + 1*k139*k75*x35 + (-1)*k139*k126*x8*x12 + (-1)*k139*k127*x8*x13 + 1*k139*k130*x5)/k139),
  Eq(Derivative(x9, t), (1*k139*k106*x50*x11 + (-1)*k139*k107*x9)/k139),
  Eq(Derivative(x10, t), (1*k139*k11*x8*x6 + (-1)*k139*k12*x10 + (-1)*k139*k15*x10*(x12 + x13) + (-1)*k139*k28*x19*x10 + 1*k139*k29*x11)/k139),
  Eq(Derivative(x11, t), ((-1)*k139*k14*x11*(x12 + x13) + 1*k139*k28*x19*x10 + (-1)*k139*k29*x11 + (-1)*k139*k30*x11*x17 + 1*k139*k31*x22 + (-1)*k139*k32*x11*x15 + 1*k139*k33*x23)/k139),
  Eq(Derivative(x12, t), ((-1)*k139*k122*x13*x12 + 1*k139*k123*x48*x54)/k139),
  Eq(Derivative(x13, t), ((-1)*k139*k124*x13*(x48 + x11) + 1*k139*k125*x55)/k139),
  Eq(Derivative(x14, t), ((-1)*k139*k17*x14*x14 + 1*k139*k22*x19*x7 + (-1)*k139*k23*x14 + (-1)*k139*k24*x14*x17 + 1*k139*k25*x20 + (-1)*k139*k26*x14*x15 + 1*k139*k27*x21)/k139),
  Eq(Derivative(x15, t), ((-1)*k139*k18*x15*x3 + 1*k139*k19*x16 + (-1)*k139*k26*x14*x15 + 1*k139*k27*x21 + (-1)*k139*k32*x11*x15 + 1*k139*k33*x23 + 1*k139*k37 + 1*k139*k38*x24 + (-1)*k139*k39*x15 + (-1)*k139*k103*x48*x15 + 1*k139*k104*x49)/k139),
  Eq(Derivative(x16, t), (1*k139*k18*x15*x3 + (-1)*k139*k19*x16)/k139),
  Eq(Derivative(x17, t), ((-1)*k139*k20*x3*x17 + 1*k139*k21*x18 + (-1)*k139*k24*x14*x17 + 1*k139*k25*x20 + (-1)*k139*k30*x11*x17 + 1*k139*k31*x22 + 1*k139*k34 + (-1)*k139*k35*x17*x14 + (-1)*k139*k36*x17*x11)/k139),
  Eq(Derivative(x18, t), (1*k139*k20*x3*x17 + (-1)*k139*k21*x18)/k139),
  Eq(Derivative(x19, t), (1*k139*k82*x36*(x14 + x11) + (-1)*k139*k84*x19*x37 + (-1)*k139*k85*x19 + (-1)*k139*k86*x19)/k139),
  Eq(Derivative(x20, t), (1*k139*k24*x14*x17 + (-1)*k139*k25*x20)/k139),
  Eq(Derivative(x21, t), (1*k139*k26*x14*x15 + (-1)*k139*k27*x21)/k139),
  Eq(Derivative(x22, t), (1*k139*k30*x11*x17 + (-1)*k139*k31*x22)/k139),
  Eq(Derivative(x23, t), (1*k139*k32*x11*x15 + (-1)*k139*k33*x23)/k139),
  Eq(Derivative(x24, t), (1*k139*k60 + 1*k139*k61*x30 + (-1)*k139*k62*x24 + (-1)*k139*x34*x31*x24)/k139),
  Eq(Derivative(x25, t), (1*k139*k40 + 1*k139*k41/(1 + k42*x26) + (-1)*k139*k43*x25 + (-1)*k139*k44*x3*x25)/k139),
  Eq(Derivative(x26, t), ((-1)*k139*k45*x5*x26 + 1*k139*k55*x29 + 1*k139*k56 + (-1)*k139*k57*x26 + 1*k139*k58/(1 + k59*x25))/k139),
  Eq(Derivative(x27, t), (1*k139*k45*x5*x26 + (-1)*k139*k46*x3*x27 + (-1)*k139*k47*x18*x27 + (-1)*k139*k48*x16*x27)/k139),
  Eq(Derivative(x28, t), (1*k139*k46*x3*x27 + 1*k139*k47*x18*x27 + 1*k139*k48*x16*x27 + (-1)*k139*k49*x14*x28 + (-1)*k139*k50*x28*x11)/k139),
  Eq(Derivative(x29, t), (1*k139*k49*x14*x28 + 1*k139*k50*x28*x11 + (-1)*k139*k55*x29)/k139),
  Eq(Derivative(x30, t), (1*k139*k78*x33 + (-1)*k139*k79*x30)/k139),
  Eq(Derivative(x31, t), (1*k139*k63 + (-1)*k139*k64*x31 + 1*k139*k66*x32**50/(k65**50 + x32**50))/k139),
  Eq(Derivative(x32, t), ((-1)*k139*k67*x32 + 1*k139*k70*x24*x33/(1 + k71*x24*x31))/k139),
  Eq(Derivative(x35, t), (1*k139*k89*x11 + (-1)*k139*k90*x35)/k139),
  Eq(Derivative(x36, t), (1*k139*k80*x5 + (-1)*k139*k81*x36*x37 + (-1)*k139*k82*x36*(x14 + x11) + (-1)*k139*k83*x36 + 1*k139*k85*x19)/k139),
  Eq(Derivative(x37, t), ((-1)*k139*k87*x37 + 1*k139*k88*x30*x38)/k139),
  Eq(Derivative(x38, t), (1*k139*k87*x37 + (-1)*k139*k88*x30*x38)/k139),
  Eq(Derivative(x39, t), (1*k139*k91*x35 + (-1)*k139*k92*x39 + (-1)*k93*x39*x40 + 1*k94*x41 + (-1)*k139*k128*x39*x12 + (-1)*k139*k129*x39*x13)/k139),
  Eq(Derivative(x40, t), ((-1)*k93*x39*x40 + 1*k94*x41 + 1*k97*x41*(x12 + x13) + 1*k98*x42*(x12 + x13))/k139),
  Eq(Derivative(x41, t), (1*k93*x39*x40 + (-1)*k94*x41 + (-1)*k95*x41*(x43 + x45) + 1*k140*k96*x42 + (-1)*k97*x41*(x12 + x13))/k140),
  Eq(Derivative(x42, t), (1*k95*x41*(x43 + x45) + (-1)*k140*k96*x42 + (-1)*k98*x42*(x12 + x13) + (-1)*k131*x42*x42 + 1*k132*x48)/k140),
  Eq(Derivative(x43, t), ((-1)*k139*k109*x43 + 1*k110*x51*(x42 + x48) + 1*k139*k112*x45 + (-1)*k139*k113*x37*x43 + (-1)*k139*k114*x43)/k139),
  Eq(Derivative(x44, t), (1*k139*k117 + 1*k139*k118*x24 + (-1)*k139*k119*x44 + (-1)*k139*k120*x52*x44)/k139),
  Eq(Derivative(x45, t), ((-1)*k139*k112*x45 + 1*k139*k113*x37*x43 + (-1)*k139*k115*x45 + 1*k116*x52*(x42 + x48))/k139),
  Eq(Derivative(x46, t), (1*k139*k99 + (-1)*k139*k100*x48*x46 + 1*k139*k101*x47)/k139),
  Eq(Derivative(x47, t), (1*k139*k100*x48*x46 + (-1)*k139*k101*x47 + (-1)*k139*k102*x47)/k139),
  Eq(Derivative(x48, t), ((-1)*k139*k103*x48*x15 + 1*k139*k104*x49 + 1*k131*x42*x42 + (-1)*k132*x48 + (-1)*k139*k133*x46*x48 + 1*k139*k134*x56*(x43 + x45) + (-1)*k139*k135*x48*x12 + (-1)*k139*k137*x48*x13)/k139),
  Eq(Derivative(x49, t), (1*k139*k103*x48*x15 + (-1)*k139*k104*x49)/k139),
  Eq(Derivative(x50, t), (1*k139*k105*x5 + (-1)*k139*k106*x50*x11)/k139),
  Eq(Derivative(x51, t), (1*k139*k108 + 1*k139*k109*x43 + (-1)*k110*x51*(x42 + x48) + (-1)*k139*k111*x37*x51)/k139),
  Eq(Derivative(x52, t), (1*k139*k111*x37*x51 + 1*k139*k115*x45 + (-1)*k116*x52*(x42 + x48) + (-1)*k139*k120*x52*x44)/k139),
  Eq(Derivative(x53, t), (1*k139*k120*x52*x44 + (-1)*k139*k121*x53)/k139),
  Eq(Derivative(x54, t), (1*k139*k122*x13*x12 + (-1)*k139*k123*x48*x54)/k139),
  Eq(Derivative(x55, t), (1*k139*k124*x13*(x48 + x11) + (-1)*k139*k125*x55)/k139),
  Eq(Derivative(x56, t), (1*k139*k133*x46*x48 + (-1)*k139*k134*x56*(x43 + x45) + (-1)*k139*k136*x56*(x12 + x13))/k139) ]