Guest User

Untitled

a guest
Mar 22nd, 2018
82
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 4.00 KB | None | 0 0
  1. I*y*(1/(x - I*y)^(5/2) - 1/(x + I*y)^(5/2))
  2.  
  3. I*Integrate[
  4. y*(1/(x - I*y)^(5/2) - 1/(x + I*y)^(5/2)), {x, -10, 10}, {y, 0, 10}]
  5. % // N
  6.  
  7. (* -(8/3) Sqrt[5] (2 Sqrt[2] + Sqrt[-1 + 5 Sqrt[2]] - [Pi] - ArcCosh[3]) *)
  8.  
  9. (* -2.31 *)
  10.  
  11. Quiet[Chop[
  12. I*NIntegrate[
  13. y *(1/(x - I y)^(5/2) - 1/(x + I y)^(5/2)), {x, -10, 10}, {y, 0,
  14. 10}, Method -> #]]] & /@ {"LocalAdaptive",
  15. {"EvenOddSubdivision", Method -> "LocalAdaptive"},
  16. "AdaptiveMonteCarlo", "QuasiMonteCarlo",
  17. "MonteCarlo", {"EvenOddSubdivision",
  18. Method -> "AdaptiveMonteCarlo"}, {"EvenOddSubdivision",
  19. Method -> "DuffyCoordinates"}}
  20.  
  21. expr = Simplify[
  22. TrigToExp[
  23. I*y*(1/(x - I*y)^(5/2) - 1/(x + I*y)^(5/2)) /. {x ->
  24. r*Cos[[CurlyPhi]], y -> r*Sin[[CurlyPhi]]}], {r > 0,
  25. 0 < [CurlyPhi] < [Pi]}]
  26.  
  27. (* (E^(-((7 I [CurlyPhi])/
  28. 2)) (-1 + E^(2 I [CurlyPhi])) (-1 + E^(5 I [CurlyPhi])))/(2 r^(
  29. 3/2)) *)
  30.  
  31. Integrate[expr2*r, {r, 0, 10}, {[CurlyPhi], 0, [Pi]}] // N
  32.  
  33. (* 2.41 *)
  34.  
  35. Quiet[NIntegrate[expr2*r, {r, 0, 10}, {[CurlyPhi], 0, [Pi]},
  36. Method -> #]] & /@ {"LocalAdaptive", {"EvenOddSubdivision",
  37. Method -> "LocalAdaptive"}, "AdaptiveMonteCarlo",
  38. "QuasiMonteCarlo",
  39. "MonteCarlo", {"EvenOddSubdivision",
  40. Method -> "AdaptiveMonteCarlo"}, {"EvenOddSubdivision",
  41. Method -> "DuffyCoordinates"}}
  42.  
  43.  
  44.  
  45. (* {2.41, 2.41, 2.43, 2.42, 2.28, 2.38, 2.41} *)
  46.  
  47. ceRe = FullSimplify[
  48. ComplexExpand[Re[I*y*(1/(x - I y)^(5/2) - 1/(x + I y)^(5/2))],
  49. TargetFunctions -> {Re, Im}], y >= 0 && x [Element] Reals]
  50.  
  51. NIntegrate[ceRe, {x, -10, 10}, {y, 0, 10}, WorkingPrecision -> 100,
  52. MaxRecursion -> 50, AccuracyGoal -> 6, PrecisionGoal -> 6,
  53. Exclusions -> {{0, 0}}]
  54.  
  55. (* 2.17329834703426526174107911927816608629584584526685923085449633834343
  56. 9718373530841484083310764572423 *)
  57.  
  58. cetr = TrigToExp[ceRe]
  59.  
  60. (* (I y)/(2 ((x - I y)/Sqrt[x^2 + y^2])^(5/2) (x^2 + y^2)^(5/4)) - (
  61. I y ((x - I y)/Sqrt[x^2 + y^2])^(5/2))/(2 (x^2 + y^2)^(5/4)) - (
  62. I y)/(2 ((x + I y)/Sqrt[x^2 + y^2])^(5/2) (x^2 + y^2)^(5/4)) + (
  63. I y ((x + I y)/Sqrt[x^2 + y^2])^(5/2))/(2 (x^2 + y^2)^(5/4)) *)
  64.  
  65.  
  66. Integrate[cetr, {x, -10, 10}, {y, 0, 10}]
  67. (N[#1, 20] &)[%]
  68.  
  69. (* (11/606 + (3 I)/
  70. 202) ((112 + 92 I) Sqrt[-10 + 10 I] + (22 - 18 I) Sqrt[-5 -
  71. 5 I] - (112 + 128 I) Sqrt[1 - 10 I] - (66 - 54 I) Sqrt[
  72. 5 - 5 I] + (176 - 144 I) Sqrt[10] + (38 + 24 I) Sqrt[
  73. 10 - 10 I] + (10 - 100 I) Sqrt[10 + 10 I] - (40 + 4 I) Sqrt[
  74. 170 - 70 I] +
  75. 16 Sqrt[1105 + 262 I] - (44 - 36 I) Sqrt[
  76. 5] [Pi] - (40 + 4 I) Sqrt[10]
  77. ArcCoth[2 Sqrt[50/101 - (5 I)/101]] + (6 - 60 I) Sqrt[10]
  78. ArcCoth[Sqrt[1 + I]] - (88 - 72 I) Sqrt[5]
  79. ArcCoth[Sqrt[2]] - (6 - 60 I) Sqrt[10]
  80. ArcCoth[2/Sqrt[2 - I/5]] + (40 + 4 I) Sqrt[10]
  81. ArcCoth[(11/101 + (9 I)/101) Sqrt[10 - 100 I]] + (16 -
  82. 160 I) Sqrt[5] ArcTan[(1 + I)/Sqrt[2]] + (16 - 160 I) Sqrt[5]
  83. ArcTanh[(1 + I)/Sqrt[2]] + (2 - 20 I) Sqrt[10]
  84. Log[1 - 1/2 (-1 + I)^(3/2)] + (18 + 22 I) (-1)^(1/4) Sqrt[5]
  85. Log[1 + 1/2 (-1 + I)^(3/2)] + (2 - 20 I) Sqrt[10]
  86. Log[1 - 1/Sqrt[1 + I]] - (2 - 20 I) Sqrt[10]
  87. Log[1 + 1/Sqrt[1 + I]] - (3 - 30 I) Sqrt[10]
  88. Log[((1 + 2 I) + 2 Sqrt[-1 + I]) (3 - 2 Sqrt[2])] + (3 -
  89. 30 I) Sqrt[10] Log[2 - Sqrt[2]] - (3 - 30 I) Sqrt[10]
  90. Log[2 + Sqrt[2]] + (1 - 10 I) Sqrt[10]
  91. Log[-Sqrt[1 - I] + Sqrt[2]] - (1 - 10 I) Sqrt[10]
  92. Log[Sqrt[1 - I] + Sqrt[2]] - (10 + I) Sqrt[-5 - 5 I] Sqrt[1 - I]
  93. Log[-Sqrt[1 + I] + Sqrt[2]] + (10 + I) Sqrt[-5 - 5 I] Sqrt[1 - I]
  94. Log[Sqrt[1 + I] + Sqrt[2]] - (3 - 30 I) Sqrt[10]
  95. Log[(-(1/2) + I/2) Sqrt[1 + 10 I] + Sqrt[10]] + (3 - 30 I) Sqrt[
  96. 10] Log[(1/2 - I/2) Sqrt[1 + 10 I] + Sqrt[10]] - (1 - 10 I) Sqrt[
  97. 10] Log[2 Sqrt[5] - Sqrt[10 - I]] + (1 - I) Sqrt[100 - 495 I]
  98. Log[2 Sqrt[5] - Sqrt[10 - I]] + (1 - 10 I) Sqrt[10]
  99. Log[2 Sqrt[5] + Sqrt[10 - I]] - (1 - I) Sqrt[100 - 495 I]
  100. Log[2 Sqrt[5] + Sqrt[10 - I]]) *)
  101.  
  102. (* 8/3 (2 Sqrt[10] - Sqrt[-5 + 25 Sqrt[2]]) *)
  103.  
  104. (* 2.1732996388253877021 + 0.*10^-20 I *)
  105.  
  106. I*Integrate[y*(1/(x - I*y)^(5/2) - 1/(x + I*y)^(5/2)), {y, 0, 10}, {x, -10, 10}]
  107. N[%]
  108. (*
  109. 8/3 Sqrt[5 (7 + 5 Sqrt[2] - 4 Sqrt[-2 + 10 Sqrt[2]])]
  110. 2.1733
  111. *)
Add Comment
Please, Sign In to add comment