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