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