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