Guest User

Untitled

a guest
Oct 19th, 2017
73
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 2.62 KB | None | 0 0
  1. expr2 = (31 (52275 arb[1] + 17425 arb[2] -
  2. 8659 (3 arb[3] + arb[4])))/(
  3. 181178952 (-arb[2] arb[3] + arb[1] arb[4])) >= 0 && (
  4. 31 (6531 arb[1] + 2177 arb[2] - 4283 (3 arb[3] + arb[4])))/(
  5. 181178952 (-arb[2] arb[3] + arb[1] arb[4])) >= 0 && (
  6. 31 (129 arb[1] + 43 arb[2] - 13 (3 arb[3] + arb[4])))/(
  7. 4213464 (-arb[2] arb[3] + arb[1] arb[4])) >= 0 && (
  8. 505325 arb[1] + 5035825 arb[2] - 8659 (29 arb[3] + 289 arb[4]))/(
  9. 60392984 (arb[2] arb[3] - arb[1] arb[4])) >= 0 && (
  10. 63133 arb[1] + 629153 arb[2] - 4283 (29 arb[3] + 289 arb[4]))/(
  11. 60392984 (arb[2] arb[3] - arb[1] arb[4])) >=
  12. 0 && (-1247 arb[1] - 12427 arb[2] + 377 arb[3] + 3757 arb[4])/(
  13. 1404488 (-arb[2] arb[3] + arb[1] arb[4])) >=
  14. 0 && (-557600 arb[1] + 2247825 arb[2] +
  15. 8659 (32 arb[3] - 129 arb[4]))/(
  16. 90589476 (arb[2] arb[3] - arb[1] arb[4])) >=
  17. 0 && (-69664 arb[1] + 280833 arb[2] +
  18. 4283 (32 arb[3] - 129 arb[4]))/(
  19. 90589476 (arb[2] arb[3] - arb[1] arb[4])) >= 0 && (
  20. 1376 arb[1] - 5547 arb[2] - 416 arb[3] + 1677 arb[4])/(
  21. 2106732 (-arb[2] arb[3] + arb[1] arb[4])) >=
  22. 0 && ((93 - 745 arb[2]) arb[3] + (31 + 745 arb[1]) arb[4])/(
  23. 838 (-arb[2] arb[3] + arb[1] arb[4])) >= 0 && (
  24. 31 (-arb[2] (-1 + arb[3]) + arb[1] (3 + arb[4])))/(
  25. 838 (arb[2] arb[3] - arb[1] arb[4])) >= 0 && (
  26. 31 (arb[2] (29 - 774 arb[3]) + 31 (3 arb[3] + arb[4]) +
  27. arb[1] (87 + 774 arb[4])))/(
  28. 216204 (-arb[2] arb[3] + arb[1] arb[4])) >= 0 && (
  29. 87 (1 + arb[2]) arb[3] + 3 (289 - 29 arb[1]) arb[4])/(
  30. 838 (arb[2] arb[3] - arb[1] arb[4])) >=
  31. 0 && (-867 arb[2] (-1 + arb[3]) + arb[1] (87 + 867 arb[4]))/(
  32. 838 (-arb[2] arb[3] + arb[1] arb[4])) >= 0 && (
  33. arb[2] (8381 - 7482 arb[3]) + 899 arb[3] + 8959 arb[4] +
  34. 29 arb[1] (29 + 258 arb[4]))/(
  35. 72068 (arb[2] arb[3] - arb[1] arb[4])) >=
  36. 0 && ((32 - 387 arb[2]) arb[3] + 129 (-1 + 3 arb[1]) arb[4])/(
  37. 419 (-arb[2] arb[3] + arb[1] arb[4])) >= 0 && (
  38. 129 arb[2] (-1 + arb[3]) + arb[1] (32 - 129 arb[4]))/(
  39. 419 (arb[2] arb[3] - arb[1] arb[4])) >= 0 && (
  40. 992 arb[3] - 129 arb[2] (29 + 64 arb[3]) - 3999 arb[4] +
  41. 32 arb[1] (29 + 258 arb[4]))/(
  42. 108102 (-arb[2] arb[3] + arb[1] arb[4])) >= 0;
  43.  
  44. vars2 = arb /@ Range[4];
  45.  
  46. cons2 = And @@ Thread[-10 <= vars2 <= 10];
  47.  
  48. (soln3 = NMinimize[{Total[List @@ expr2[[All, 1]]], expr2 && cons2},
  49. vars2, Reals]) // Quiet // AbsoluteTiming
  50.  
  51. expr2 && cons2 /. soln3[[2]]
  52.  
  53. (soln4 = NMaximize[{Total[List @@ expr2[[All, 1]]], expr2 && cons2},
  54. vars2, Reals]) // Quiet // AbsoluteTiming
  55.  
  56. expr2 && cons2 /. soln4[[2]]
  57.  
  58. soleg1=Thread[vars2 -> {49/160, 46837/64500, 7071/4000, 25619/129000}];
  59. soleg2=Thread[vars2 -> {151/2000, 88147/258000, 4019/4000, -(67417/516000)}];
  60. expr2 && cons2 /. soleg1 (* Returns True *)
  61. expr2 && cons2 /. soleg2 (* Returns True *)
Add Comment
Please, Sign In to add comment