[ Diff(x[1], t) = (1*k[40]*k[31]*x[1]*(1 - x[1]/k[32]) + (-1)*k[40]*k[36]*k[34]*x[1]/(k[35] + k[33]*x[1]) + (-1)*k[40]*k[1]*x[1]*x[4]*(1 + k[25]*k[38]) + (-1)*k[40]*k[2]*x[1]*x[3]/(1 + (x[4]/k[3])^2 + (k[37]/k[29])^2)*(1 + k[27]*k[38]))/k[40],
  Diff(x[2], t) = ((-1)*k[40]*k[4]*x[2]*x[4]*(1 + k[25]*x[4]) + (-1)*k[40]*k[5]*x[2]*x[3]/(1 + (x[4]/k[3])^2 + (k[37]/k[29])^2)*(1 + k[27]*k[38]) + 1*k[40]*k[11]*x[4] + (-1)*k[40]*k[6]*x[2])/k[40],
  Diff(x[3], t) = (1*k[40]*k[17]*(k[16]*x[2] + k[13]*x[1] + k[14]/(1 + (k[38]/k[26])^2)*x[3] + k[15]*x[4])/(1 + (k[37]/k[30])^2)/(k[18] + (k[16]*x[2] + k[13]*x[1] + k[14]/(1 + (k[38]/k[26])^2)*x[3] + k[15]*x[4])/(1 + (k[37]/k[30])^2)) + (-1)*k[40]*k[19]*x[3])/k[40],
  Diff(x[4], t) = ((-1)*k[40]*k[11]*x[4] + 1*k[40]*k[9]*(k[8]*x[2] + k[7]*x[1])/(1 + (k[37]/k[28])^2)/(k[10] + (k[8]*x[2] + k[7]*x[1])/(1 + (k[37]/k[28])^2))/(1 + (k[38]/k[24])^2) + (-1)*k[40]*k[12]*x[4]*x[3]/(1 + (x[4]/k[3])^2 + (k[37]/k[29])^2))/k[40] ]