{ diff(x1, t) = (1*k12*x4 + (-1)*k8*x2*x1 + 1*k9*x3 + (-1)*k22*x1*x17 + 1*k21*x28*x21 + 1*k18*x27*x18 + (-1)*k31*x1*x6 + 1*k32*x9)/k59,
  diff(x2, t) = ((-1)*k8*x2*x1 + 1*k9*x3 + (-1)*k10*x39*x2 + 1*k11*x40 + 1*k6*x5 + (-1)*k29*x2*x6 + 1*k30*x8)/k59,
  diff(x3, t) = (1*k8*x2*x1 + (-1)*k9*x3 + (-1)*k16*x17*x3 + 1*k19*x23*x20)/k59,
  diff(x4, t) = (1*k3*x2 + 1*k3*x8 + 1*k5*x40 + 1*k5*x41 + (-1)*k4*x4)/k59,
  diff(x5, t) = (1*k1*k63 + (-1)*k2*x5 + 1*k52*x42)/k59,
  diff(x6, t) = (1*k27*x13*x7 + (-1)*k28*x6)/k59,
  diff(x7, t) = ((-1)*k27*x13*x7 + 1*k28*x6)/k59,
  diff(x8, t) = ((-1)*k10*x39*x8 + 1*k11*x41 + 1*k29*x2*x6 + (-1)*k30*x8)/k59,
  diff(x9, t) = (1*k31*x1*x6 + (-1)*k32*x9 + (-1)*k23*x9*x17 + 1*k21*x32*x21)/k59,
  diff(x11, t) = (1*k55*x44 + 1*k55*x43)/k59,
  diff(x12, t) = 0/k59,
  diff(x13, t) = (1*k25*x10 + (-1)*k26*x13 + 1*k53*x11 + 1*k54*x12)/k59,
  diff(x14, t) = ((-1)*k14*x14*x19*k60/(5000 + k60) + 1*k15*x15*x16)/k59,
  diff(x15, t) = ((-1)*k15*x15*x16 + 1*k22*x1*x17 + 1*k24*x28*x17 + 1*k24*x29*x17 + 1*k24*x30*x17 + 1*k16*x17*x3 + 1*k17*x23*x17 + 1*k17*x24*x17 + 1*k17*x25*x17 + 1*k23*x9*x17 + 1*k24*x32*x17 + 1*k24*x33*x17 + 1*k24*x34*x17)/k59,
  diff(x16, t) = (1*k14*x14*x19*k60/(5000 + k60) + (-1)*k15*x15*x16)/k59,
  diff(x17, t) = (1*k15*x15*x16 + (-1)*k22*x1*x17 + (-1)*k24*x28*x17 + (-1)*k24*x29*x17 + (-1)*k24*x30*x17 + (-1)*k16*x17*x3 + (-1)*k17*x23*x17 + (-1)*k17*x24*x17 + (-1)*k17*x25*x17 + (-1)*k23*x9*x17 + (-1)*k24*x32*x17 + (-1)*k24*x33*x17 + (-1)*k24*x34*x17)/k59,
  diff(x18, t) = ((-1)*k18*x31*x18 + 1*k13*x37*k56 + (-1)*k18*x27*x18 + 1*k7*k56*x36*k60/(5000 + k60) + (-1)*k48*x50*x18 + (-1)*k48*x44*x18 + (-1)*k18*x35*x18 + 1*k13*x38*k56 + (-1)*k41*x46*x18 + 1*k42*x52)/k59,
  diff(x19, t) = ((-1)*k14*x14*x19*k60/(5000 + k60) + 1*k21*x31*x21 + 1*k21*x30*x21 + 1*k21*x29*x21 + 1*k21*x28*x21 + 4*k13*x37*k56 + 1*k19*x26*x20 + 1*k19*x25*x20 + 1*k19*x24*x20 + 1*k19*x23*x20 + 4*k7*k56*x36*k60/(5000 + k60) + 1*k21*x35*x21 + 1*k21*x34*x21 + 1*k21*x33*x21 + 1*k21*x32*x21 + 4*k13*x38*k56)/k59,
  diff(x20, t) = 0/k59,
  diff(x21, t) = 0/k59,
  diff(x22, t) = 0/k59,
  diff(x23, t) = (1*k16*x17*x3 + (-1)*k17*x23*x17 + 1*k19*x24*x20 + (-1)*k19*x23*x20)/k59,
  diff(x24, t) = (1*k17*x23*x17 + (-1)*k17*x24*x17 + 1*k19*x25*x20 + (-1)*k19*x24*x20)/k59,
  diff(x25, t) = (1*k17*x24*x17 + (-1)*k17*x25*x17 + 1*k19*x26*x20 + (-1)*k19*x25*x20)/k59,
  diff(x26, t) = (1*k17*x25*x17 + (-1)*k19*x26*x20 + (-1)*k33*x26*x39 + (-1)*k34*x26*x40 + (-1)*k34*x26*x41)/k59,
  diff(x27, t) = (1*k33*x26*x39 + 1*k34*x26*x40 + 1*k34*x26*x41 + (-1)*k18*x27*x18)/k59,
  diff(x28, t) = (1*k22*x1*x17 + (-1)*k24*x28*x17 + 1*k21*x29*x21 + (-1)*k21*x28*x21)/k59,
  diff(x29, t) = (1*k24*x28*x17 + (-1)*k24*x29*x17 + 1*k21*x30*x21 + (-1)*k21*x29*x21)/k59,
  diff(x30, t) = (1*k24*x29*x17 + (-1)*k24*x30*x17 + 1*k21*x31*x21 + (-1)*k21*x30*x21)/k59,
  diff(x31, t) = (1*k24*x30*x17 + (-1)*k21*x31*x21 + (-1)*k18*x31*x18)/k59,
  diff(x32, t) = (1*k23*x9*x17 + (-1)*k24*x32*x17 + 1*k21*x33*x21 + (-1)*k21*x32*x21)/k59,
  diff(x33, t) = (1*k24*x32*x17 + (-1)*k24*x33*x17 + 1*k21*x34*x21 + (-1)*k21*x33*x21)/k59,
  diff(x34, t) = (1*k24*x33*x17 + (-1)*k24*x34*x17 + 1*k21*x35*x21 + (-1)*k21*x34*x21)/k59,
  diff(x35, t) = (1*k24*x34*x17 + (-1)*k21*x35*x21 + (-1)*k18*x35*x18)/k59,
  diff(x36, t) = (1*k18*x27*x18 + (-1)*k7*k56*x36*k60/(5000 + k60))/k59,
  diff(x37, t) = (1*k18*x31*x18 + (-1)*k13*x37*k56)/k59,
  diff(x38, t) = (1*k18*x35*x18 + (-1)*k13*x38*k56)/k59,
  diff(x39, t) = ((-1)*k10*x39*x2 + 1*k11*x40 + (-1)*k10*x39*x8 + 1*k11*x41)/k59,
  diff(x40, t) = (1*k10*x39*x2 + (-1)*k11*x40)/k59,
  diff(x41, t) = (1*k10*x39*x8 + (-1)*k11*x41)/k59,
  diff(x42, t) = (1*k47*x40 + 1*k47*x41 + (-1)*k49*x42 + (-2)*k50*x42*(x42 - 1)*1/2 + (-1)*k50*x42*x44)/k59,
  diff(x43, t) = 1*k48*x44*x18/k59,
  diff(x44, t) = ((-1)*k48*x44*x18 + 1*k50*x42*(x42 - 1)*1/2 + 1*k50*x42*x44 + (-2)*k51*x44*(x44 - 1)*1/2 + (-1)*k51*x44*x45)/k59,
  diff(x45, t) = (2*k51*x44*(x44 - 1)*1/2 + 1*k51*x44*x45)/k59,
  diff(x46, t) = ((-1)*k38*x46 + 1*k39*x49 + (-1)*k35*x40*x46 + (-1)*k35*x41*x46 + (-1)*k36*x39*x46 + 1*k37*x47*x53 + (-2)*k43*x46*(x46 - 1)*1/2 + (-1)*k43*x46*x50 + 1*k40*k63 + (-1)*k41*x46*x18)/k59,
  diff(x47, t) = (1*k35*x40*x46 + (-1)*k35*x40*x47 + 1*k35*x41*x46 + (-1)*k35*x41*x47 + 1*k36*x39*x46 + (-1)*k36*x39*x47 + 1*k37*x48*x53 + (-1)*k37*x47*x53 + (-2)*k44*x47*(x47 - 1)*1/2 + (-1)*k44*x47*x50)/k59,
  diff(x48, t) = (1*k35*x40*x47 + 1*k35*x41*x47 + 1*k36*x39*x47 + (-1)*k37*x48*x53 + (-2)*k45*x48*(x48 - 1)*1/2 + (-1)*k45*x48*x50)/k59,
  diff(x49, t) = (1*k38*x46 + (-1)*k39*x49)/k59,
  diff(x50, t) = (2*k44*x47*(x47 - 1)*1/2 + 1*k44*x47*x50 + 2*k45*x48*(x48 - 1)*1/2 + 1*k45*x48*x50 + 2*k43*x46*(x46 - 1)*1/2 + 1*k43*x46*x50 + (-2)*k46*x50*(x50 - 1)*1/2 + (-1)*k46*x50*x54 + (-1)*k48*x50*x18)/k59,
  diff(x51, t) = 1*k48*x50*x18/k59,
  diff(x52, t) = (1*k41*x46*x18 + (-1)*k42*x52)/k59,
  diff(x53, t) = 0/k59,
  diff(x54, t) = (2*k46*x50*(x50 - 1)*1/2 + 1*k46*x50*x54)/k59,
  diff(x55, t) = 0,
  diff(x56, t) = 0,
  diff(x57, t) = 0,
  diff(x58, t) = 0,
  diff(x59, t) = 0 }