Advertisement
Guest User

Untitled

a guest
Mar 30th, 2017
61
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 1.52 KB | None | 0 0
  1. f1 = 1 - z[2, 3]^2 + z[1, 4]^2 (-1 + z[2, 3]^2) - z[2, 4]^2 +
  2. z[1, 3]^2 (-1 + z[2, 4]^2) - 2 z[2, 3] z[2, 4] z[3, 4] - z[3, 4]^2
  3. - 2 z[1, 3] z[1, 4] (z[2, 3] z[2, 4] + z[3, 4]) +
  4. z[1, 2]^2 (-1 + z[3, 4]^2) +
  5. 2 z[1, 2] (z[1, 4] (z[2, 4] + z[2, 3] z[3, 4]) +
  6. z[1, 3] (z[2, 3] + z[2, 4] z[3, 4]))
  7.  
  8. f2 = 1 - z[1, 2]^2 - z[1, 3]^2 + 2 z[1, 2] z[1, 3] z[2, 3] - z[2, 3]^2
  9.  
  10. f3 = 1 - z[1, 2]^2
  11.  
  12. Integrate[Boole[f1 > 0 && f2 > 0 && f3 > 0], {z[1, 2], -1, 1}, {z[1, 3], -1, 1},
  13. {z[1, 4], -1, 1}, {z[2, 3], -1, 1}, {z[2, 4], -1, 1}, {z[3, 4], -1, 1}]
  14.  
  15. J6 = GenericCylindricalDecomposition[
  16. f1 > 0 && f2 > 0 && f3 > 0, {z[1, 2], z[1, 2], z[1, 3], z[1, 4],
  17. z[2, 3], z[2, 4], z[3, 4]}];
  18.  
  19. Integrate[
  20. Boole[J6], {y[1, 2], -1, 1}, {z[1, 2], -1, 1}, {z[1, 3], -1, 1},
  21. {z[1, 4], -1, 1}, {z[2, 3], -1, 1}, {z[2, 4], -1, 1}, {z[3, 4], -1, 1}]
  22.  
  23. f4 = -u^4 z[1, 4]^2 - z[2, 3]^2 +
  24. 2 u z[2, 3] (z[1, 2] z[2, 4] - z[1, 3] z[3, 4])
  25. + 2 u^3 z[1, 4] (z[1, 2] z[1, 3] - z[2, 4] z[3, 4])
  26. - u^2 (-1 - z[1, 4]^2 z[2, 3]^2 + 2 z[1, 3] z[1, 4] z[2, 3] z[2, 4] +
  27. z[2, 4]^2 - z[1, 3]^2 (-1 + z[2, 4]^2) -
  28. 2 z[1, 2] (z[1, 4] z[2, 3] + z[1, 3] z[2, 4]) z[3, 4] +
  29. z[3, 4]^2 - z[1, 2]^2 (-1 + z[3, 4]^2))
  30.  
  31. f5 = 1 - z[1, 2]^2 - z[1, 3]^2 + 2 u z[1, 2] z[1, 3] z[1, 4] - u^2 z[1, 4]^2,
  32.  
  33. f6 = u^2 - u^2 z[1, 2]^2 - z[2, 3]^2 + 2 u z[1, 2] z[2, 3] z[2, 4] - u^2 z[2, 4]^2
  34.  
  35. f7 = u^2 - u^2 z[1, 3]^2 - z[2, 3]^2 + 2 u z[1, 3] z[2, 3] z[3, 4] - u^2 z[3, 4]^2
  36.  
  37. f8 = 1 - u^2 z[1, 4]^2 - z[2, 4]^2 + 2 u z[1, 4] z[2, 4] z[3, 4] - z[3, 4]^2
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement