[ x[4] = x[1] - x[2],
  x[6] = k[18] - x[3],
  x[9] = 1/2*(x[2] - k[19] - k[2]/k[3]*(k[3] + x[3] + k[6]*(k[18] - x[3])) + sqrt((k[2]/k[3]*(k[3] + x[3] + k[6]*(k[18] - x[3])) - x[2] + k[19])^2 + 4*x[2]*k[2]/k[3]*(k[3] + x[3] + k[6]*(k[18] - x[3])))),
  x[10] = 1/2*(x[1] - k[20] - k[4] + sqrt((k[4] - x[1] + k[20])^2 + 4*x[1]*k[4])),
  x[12] = (k[2] + x[9])/(k[2]/k[3]*(k[3] + x[3] + k[6]*(k[18] - x[3])) + x[9]),
  x[13] = x[2]/(k[5] + x[1]),
  x[14] = x[1]/(k[5] + x[1]),
  x[15] = (x[1] - x[2])/(k[5] + x[1]),
  x[16] = x[3]/k[18] ]