{ diff(x1, t) = (-1)*(k2*x1*x2 - k3*x3)/k1,
  diff(x2, t) = ((-1)*(k2*x1*x2 - k3*x3) + 1*k94*x58)/k1,
  diff(x3, t) = (1*(k2*x1*x2 - k3*x3) + (-1)*(k4*x3*x3 - k5*x4))/k1,
  diff(x4, t) = (1*(k4*x3*x3 - k5*x4) + (-1)*k6*x4 + 1*k12*x10 + 1*k23*x16 + 1*k34*x22 + 1*k56*x33)/k1,
  diff(x5, t) = (1*k6*x4 + (-1)*(k7*x5*x7 - k8*x8) + (-1)*(k9*x5*x9 - k10*x10) + 1*(k13*x11 - k14*x5*x12) + (-1)*(k30*x5*x14 - k31*x21) + (-1)*(k37*x5*x20 - k38*x23) + (-1)*(k49*x5*x29 - k50*x31) + 1*k89*x55 + 1*k92*x57 + (-1)*(k105*x5*x6 - k106*x62) + (-1)*(k110*x5*x64 - k111*x65) + 1*(k114*x65 - k115*x5*x67) + (-1)*(k151*x5*x90 - k152*x91) + 1*k153*x91 + (-1)*(k154*x5*x92 - k155*x93))/k1,
  diff(x6, t) = ((-1)*(k95*x18*x6 - k96*x19) + 1*k99*x59 + (-1)*(k100*x23*x6 - k101*x24) + 1*k104*x61 + (-1)*(k105*x5*x6 - k106*x62) + 1*k109*x63 + (-1)*(k178*x91*x6 - k179*x106) + 1*k182*x107 + (-1)*(k183*x65*x6 - k184*x108))/k1,
  diff(x7, t) = ((-1)*(k7*x5*x7 - k8*x8) + 1*k17*x13 + 1*k18*x12)/k1,
  diff(x8, t) = (1*(k7*x5*x7 - k8*x8) + (-1)*k11*x8)/k1,
  diff(x9, t) = ((-1)*(k9*x5*x9 - k10*x10) + 1*k12*x10 + (-1)*(k15*x12*x9 - k16*x13) + 1*k17*x13 + (-1)*(k21*x15*x9 - k22*x16) + 1*k23*x16 + (-1)*(k32*x21*x9 - k33*x22) + 1*k34*x22 + (-1)*(k54*x31*x9 - k55*x33) + 1*k56*x33)/k1,
  diff(x10, t) = (1*(k9*x5*x9 - k10*x10) + (-1)*k12*x10)/k1,
  diff(x11, t) = (1*k11*x8 + (-1)*(k13*x11 - k14*x5*x12) + (-1)*(k19*x11*x14 - k20*x15) + (-1)*(k28*x11*x20 - k29*x18))/k1,
  diff(x12, t) = (1*(k13*x11 - k14*x5*x12) + (-1)*(k15*x12*x9 - k16*x13) + (-1)*k18*x12 + 1*k23*x16 + 1*k89*x55 + 1*k99*x59)/k1,
  diff(x13, t) = (1*(k15*x12*x9 - k16*x13) + (-1)*k17*x13)/k1,
  diff(x14, t) = ((-1)*(k19*x11*x14 - k20*x15) + 1*k23*x16 + (-1)*(k26*x14*x17 - k27*x20) + (-1)*(k30*x5*x14 - k31*x21) + 1*k34*x22 + 1*k89*x55 + 1*k92*x57)/k1,
  diff(x15, t) = (1*(k19*x11*x14 - k20*x15) + (-1)*(k21*x15*x9 - k22*x16) + (-1)*(k24*x15*x17 - k25*x18))/k1,
  diff(x16, t) = (1*(k21*x15*x9 - k22*x16) + (-1)*k23*x16)/k1,
  diff(x17, t) = ((-1)*(k24*x15*x17 - k25*x18) + (-1)*(k26*x14*x17 - k27*x20) + (-1)*(k35*x21*x17 - k36*x23) + 1*k93*x56)/k1,
  diff(x18, t) = (1*(k24*x15*x17 - k25*x18) + 1*(k28*x11*x20 - k29*x18) + (-1)*(k39*x18*x25 - k40*x26) + 1*k42*x26 + (-1)*(k87*x46*x18 - k88*x55) + (-1)*(k95*x18*x6 - k96*x19))/k1,
  diff(x19, t) = (1*(k95*x18*x6 - k96*x19) + (-1)*(k97*x19*x60 - k98*x59))/k1,
  diff(x20, t) = (1*(k26*x14*x17 - k27*x20) + (-1)*(k28*x11*x20 - k29*x18) + (-1)*(k37*x5*x20 - k38*x23) + 1*k99*x59 + 1*k104*x61)/k1,
  diff(x21, t) = (1*(k30*x5*x14 - k31*x21) + (-1)*(k32*x21*x9 - k33*x22) + (-1)*(k35*x21*x17 - k36*x23))/k1,
  diff(x22, t) = (1*(k32*x21*x9 - k33*x22) + (-1)*k34*x22)/k1,
  diff(x23, t) = (1*(k35*x21*x17 - k36*x23) + 1*(k37*x5*x20 - k38*x23) + 1*k43*x28 + (-1)*(k44*x23*x25 - k45*x28) + (-1)*(k90*x46*x23 - k91*x57) + (-1)*(k100*x23*x6 - k101*x24))/k1,
  diff(x24, t) = (1*(k100*x23*x6 - k101*x24) + (-1)*(k102*x24*x60 - k103*x61))/k1,
  diff(x25, t) = ((-1)*(k39*x18*x25 - k40*x26) + 1*k41*x27 + (-1)*(k44*x23*x25 - k45*x28) + 1*k48*x30 + 1*k53*x32)/k1,
  diff(x26, t) = (1*(k39*x18*x25 - k40*x26) + (-1)*k42*x26)/k1,
  diff(x27, t) = ((-1)*k41*x27 + 1*k42*x26 + 1*k43*x28 + (-1)*(k46*x27*x29 - k47*x30) + (-1)*(k51*x31*x27 - k52*x32) + (-1)*(k57*x34*x27 - k58*x35) + 1*k59*x35)/k1,
  diff(x28, t) = ((-1)*k43*x28 + 1*(k44*x23*x25 - k45*x28))/k1,
  diff(x29, t) = ((-1)*(k46*x27*x29 - k47*x30) + 1*k48*x30 + (-1)*(k49*x5*x29 - k50*x31) + 1*k56*x33)/k1,
  diff(x30, t) = (1*(k46*x27*x29 - k47*x30) + (-1)*k48*x30)/k1,
  diff(x31, t) = (1*(k49*x5*x29 - k50*x31) + (-1)*(k51*x31*x27 - k52*x32) + 1*k53*x32 + (-1)*(k54*x31*x9 - k55*x33))/k1,
  diff(x32, t) = (1*(k51*x31*x27 - k52*x32) + (-1)*k53*x32)/k1,
  diff(x33, t) = (1*(k54*x31*x9 - k55*x33) + (-1)*k56*x33)/k1,
  diff(x34, t) = ((-1)*(k57*x34*x27 - k58*x35) + 1*k74*x48)/k1,
  diff(x35, t) = (1*(k57*x34*x27 - k58*x35) + (-1)*k59*x35)/k1,
  diff(x36, t) = (1*k59*x35 + (-1)*(k60*x36*x37 - k61*x38) + 1*k62*x38 + (-1)*(k63*x36*x39 - k64*x40) + 1*k65*x40 + (-1)*(k72*x36*x47 - k73*x48) + (-1)*k148*x88*x36/(k149 + x36) + 1*k150*x89)/k1,
  diff(x37, t) = ((-1)*(k60*x36*x37 - k61*x38) + 1*k80*x51)/k1,
  diff(x38, t) = (1*(k60*x36*x37 - k61*x38) + (-1)*k62*x38)/k1,
  diff(x39, t) = (1*k62*x38 + (-1)*(k63*x36*x39 - k64*x40) + 1*k77*x50 + (-1)*(k78*x39*x49 - k79*x51))/k1,
  diff(x40, t) = (1*(k63*x36*x39 - k64*x40) + (-1)*k65*x40)/k1,
  diff(x41, t) = (1*k65*x40 + (-1)*(k66*x41*x42 - k67*x43) + 1*k68*x43 + (-1)*(k69*x41*x44 - k70*x45) + 1*k71*x45 + (-1)*(k75*x41*x49 - k76*x50))/k1,
  diff(x42, t) = ((-1)*(k66*x41*x42 - k67*x43) + 1*k84*x53)/k1,
  diff(x43, t) = (1*(k66*x41*x42 - k67*x43) + (-1)*k68*x43)/k1,
  diff(x44, t) = (1*k68*x43 + (-1)*(k69*x41*x44 - k70*x45) + 1*k83*x54 + (-1)*(k85*x44*x52 - k86*x53))/k1,
  diff(x45, t) = (1*(k69*x41*x44 - k70*x45) + (-1)*k71*x45)/k1,
  diff(x46, t) = (1*k71*x45 + (-1)*(k81*x46*x52 - k82*x54) + (-1)*(k87*x46*x18 - k88*x55) + 1*k89*x55 + (-1)*(k90*x46*x23 - k91*x57) + 1*k92*x57)/k1,
  diff(x47, t) = ((-1)*(k72*x36*x47 - k73*x48) + 1*k74*x48)/k1,
  diff(x48, t) = (1*(k72*x36*x47 - k73*x48) + (-1)*k74*x48)/k1,
  diff(x49, t) = ((-1)*(k75*x41*x49 - k76*x50) + 1*k77*x50 + (-1)*(k78*x39*x49 - k79*x51) + 1*k80*x51)/k1,
  diff(x50, t) = (1*(k75*x41*x49 - k76*x50) + (-1)*k77*x50)/k1,
  diff(x51, t) = (1*(k78*x39*x49 - k79*x51) + (-1)*k80*x51)/k1,
  diff(x52, t) = ((-1)*(k81*x46*x52 - k82*x54) + 1*k83*x54 + 1*k84*x53 + (-1)*(k85*x44*x52 - k86*x53))/k1,
  diff(x53, t) = ((-1)*k84*x53 + 1*(k85*x44*x52 - k86*x53))/k1,
  diff(x54, t) = (1*(k81*x46*x52 - k82*x54) + (-1)*k83*x54)/k1,
  diff(x55, t) = (1*(k87*x46*x18 - k88*x55) + (-1)*k89*x55)/k1,
  diff(x56, t) = (1*k89*x55 + 1*k92*x57 + (-1)*k93*x56)/k1,
  diff(x57, t) = (1*(k90*x46*x23 - k91*x57) + (-1)*k92*x57)/k1,
  diff(x58, t) = (-1)*k94*x58/k1,
  diff(x59, t) = (1*(k97*x19*x60 - k98*x59) + (-1)*k99*x59)/k1,
  diff(x60, t) = ((-1)*(k97*x19*x60 - k98*x59) + 1*k99*x59 + (-1)*(k102*x24*x60 - k103*x61) + 1*k104*x61 + (-1)*(k107*x62*x60 - k108*x63) + 1*k109*x63 + (-1)*(k180*x106*x60 - k181*x107) + 1*k182*x107 + (-1)*(k185*x108*x60 - k186*x109) + 1*k187*x109)/k1,
  diff(x61, t) = (1*(k102*x24*x60 - k103*x61) + (-1)*k104*x61)/k1,
  diff(x62, t) = (1*(k105*x5*x6 - k106*x62) + (-1)*(k107*x62*x60 - k108*x63))/k1,
  diff(x63, t) = (1*(k107*x62*x60 - k108*x63) + (-1)*k109*x63)/k1,
  diff(x64, t) = ((-1)*(k110*x5*x64 - k111*x65) + 1*(k119*x70 - k120*x68*x64) + 1*k187*x109)/k1,
  diff(x65, t) = (1*(k110*x5*x64 - k111*x65) + (-1)*(k112*x65 - k113*x66) + (-1)*(k114*x65 - k115*x5*x67) + (-1)*(k183*x65*x6 - k184*x108))/k1,
  diff(x66, t) = 1*(k112*x65 - k113*x66)/k1,
  diff(x67, t) = (1*(k114*x65 - k115*x5*x67) + (-1)*(k116*x68*x67 - k117*x69) + (-1)*(k121*x67*x71 - k122*x72) + 1*k123*x72)/k1,
  diff(x68, t) = ((-1)*(k116*x68*x67 - k117*x69) + 1*(k119*x70 - k120*x68*x64))/k1,
  diff(x69, t) = (1*(k116*x68*x67 - k117*x69) + (-1)*k118*x69)/k1,
  diff(x70, t) = (1*k118*x69 + (-1)*(k119*x70 - k120*x68*x64))/k1,
  diff(x71, t) = ((-1)*(k121*x67*x71 - k122*x72) + 1*k188*x73)/k1,
  diff(x72, t) = (1*(k121*x67*x71 - k122*x72) + (-1)*k123*x72)/k1,
  diff(x73, t) = (1*k123*x72 + (-1)*(k124*x74*x73 - k125*x75) + 1*(k131*x79 - k132*x80*x73) + (-1)*k188*x73)/k1,
  diff(x74, t) = ((-1)*(k124*x74*x73 - k125*x75) + 1*(k134*x82 - k135*x74*x83))/k1,
  diff(x75, t) = (1*(k124*x74*x73 - k125*x75) + (-1)*(k126*x75*x76 - k127*x77) + 1*(k141*x85 - k142*x75*x83))/k1,
  diff(x76, t) = ((-1)*(k126*x75*x76 - k127*x77) + 1*(k129*x78 - k130*x79*x76))/k1,
  diff(x77, t) = (1*(k126*x75*x76 - k127*x77) + (-1)*k128*x77 + 1*(k146*x87 - k147*x77*x83))/k1,
  diff(x78, t) = (1*k128*x77 + (-1)*(k129*x78 - k130*x79*x76) + (-1)*(k143*x78*x83 - k144*x86))/k1,
  diff(x79, t) = (1*(k129*x78 - k130*x79*x76) + (-1)*(k131*x79 - k132*x80*x73) + (-1)*(k138*x79*x83 - k139*x84))/k1,
  diff(x80, t) = (1*(k131*x79 - k132*x80*x73) + (-1)*(k136*x80*x83 - k137*x81))/k1,
  diff(x81, t) = ((-1)*k133*x81 + 1*(k136*x80*x83 - k137*x81))/k1,
  diff(x82, t) = (1*k133*x81 + (-1)*(k134*x82 - k135*x74*x83))/k1,
  diff(x83, t) = (1*(k134*x82 - k135*x74*x83) + (-1)*(k136*x80*x83 - k137*x81) + (-1)*(k138*x79*x83 - k139*x84) + 1*(k141*x85 - k142*x75*x83) + (-1)*(k143*x78*x83 - k144*x86) + 1*(k146*x87 - k147*x77*x83))/k1,
  diff(x84, t) = (1*(k138*x79*x83 - k139*x84) + (-1)*k140*x84)/k1,
  diff(x85, t) = (1*k140*x84 + (-1)*(k141*x85 - k142*x75*x83))/k1,
  diff(x86, t) = (1*(k143*x78*x83 - k144*x86) + (-1)*k145*x86)/k1,
  diff(x87, t) = (1*k145*x86 + (-1)*(k146*x87 - k147*x77*x83))/k1,
  diff(x89, t) = (1*k148*x88*x36/(k149 + x36) + (-1)*k150*x89)/k1,
  diff(x90, t) = ((-1)*(k151*x5*x90 - k152*x91) + 1*k158*x95 + (-1)*(k164*x90*x92 - k165*x98) + 1*k177*x104 + 1*k182*x107)/k1,
  diff(x91, t) = (1*(k151*x5*x90 - k152*x91) + (-1)*k153*x91 + (-1)*(k178*x91*x6 - k179*x106))/k1,
  diff(x92, t) = (1*k153*x91 + (-1)*(k154*x5*x92 - k155*x93) + (-1)*(k156*x92*x94 - k157*x95) + (-1)*(k159*x92*x92 - k160*x97) + (-1)*(k164*x90*x92 - k165*x98))/k1,
  diff(x93, t) = 1*(k154*x5*x92 - k155*x93)/k1,
  diff(x94, t) = ((-1)*(k156*x92*x94 - k157*x95) + 1*k158*x95 + (-1)*(k161*x97*x94 - k162*x96) + 1*k163*x96)/k1,
  diff(x95, t) = (1*(k156*x92*x94 - k157*x95) + (-1)*k158*x95)/k1,
  diff(x96, t) = (1*(k161*x97*x94 - k162*x96) + (-1)*k163*x96)/k1,
  diff(x97, t) = (1*(k159*x92*x92 - k160*x97) + (-1)*(k161*x97*x94 - k162*x96) + (-1)*k166*x97)/k1,
  diff(x98, t) = (1*k163*x96 + 1*(k164*x90*x92 - k165*x98))/k1,
  diff(x99, t) = (1*k166*x97 + 1*(k167*x100*x100 - k168*x99) + (-1)*(k169*x99*x101 - k170*x102))/k1,
  diff(x100, t) = ((-1)*(k167*x100*x100 - k168*x99) + (-1)*(k172*x104*x100 - k173*x103) + (-1)*(k174*x100*x101 - k175*x105))/k1,
  diff(x101, t) = ((-1)*(k169*x99*x101 - k170*x102) + 1*k171*x102 + (-1)*(k174*x100*x101 - k175*x105) + 1*k176*x105)/k1,
  diff(x102, t) = (1*(k169*x99*x101 - k170*x102) + (-1)*k171*x102)/k1,
  diff(x103, t) = (1*k171*x102 + 1*(k172*x104*x100 - k173*x103))/k1,
  diff(x104, t) = ((-1)*(k172*x104*x100 - k173*x103) + 1*k176*x105 + (-1)*k177*x104)/k1,
  diff(x105, t) = (1*(k174*x100*x101 - k175*x105) + (-1)*k176*x105)/k1,
  diff(x106, t) = (1*(k178*x91*x6 - k179*x106) + (-1)*(k180*x106*x60 - k181*x107))/k1,
  diff(x107, t) = (1*(k180*x106*x60 - k181*x107) + (-1)*k182*x107)/k1,
  diff(x108, t) = (1*(k183*x65*x6 - k184*x108) + (-1)*(k185*x108*x60 - k186*x109))/k1,
  diff(x109, t) = (1*(k185*x108*x60 - k186*x109) + (-1)*k187*x109)/k1 }