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