Timeline for Error with NIntegrate of a two variables integration
Current License: CC BY-SA 4.0
9 events
| when toggle format | what | by | license | comment | |
|---|---|---|---|---|---|
| Mar 19, 2021 at 12:08 | history | edited | Dominic | CC BY-SA 4.0 |
added def of function to updated code
|
| Mar 19, 2021 at 12:06 | comment | added | Dominic | @Daniele: Thanks. That reduces the time on my 3.5 GHz machine to about 24 seconds. | |
| Mar 19, 2021 at 9:20 | comment | added | Daniele Binosi |
I get an enormous speed up of your code when: removing everywhere MaxRecursion -> 25, WorkingPrecision -> 30; changing the integration method to Method -> "LocalAdaptive". Then the result I get are (*-2.6879624115154894*) and (*-3.134223337637785*)
|
|
| Mar 19, 2021 at 6:23 | comment | added | user64494 | +1. No comments. | |
| Mar 18, 2021 at 21:00 | history | edited | Dominic | CC BY-SA 4.0 |
Added second attempt.
|
| Mar 18, 2021 at 18:10 | comment | added | user64494 |
Look at Plot3D[{Im[((8 - 8 I) x y (1 + x^2 + y^2) ((3 + 4 I) + (2 + 4 I) x^2 - x^4 - 2 ((1 + 2 I) + x^2) y^2 + 3 y^4))/(\[Pi] (-((-2 - I) + x)^2 + y^2) (-((2 + I) + x)^2 + y^2) ((1 + x^2)^2 - 2 (-1 + x^2) y^2 + y^4)^2)]}, {y, 0, 100}, {x, -100, 0}, PlotRange -> All, PlotPoints -> 100]. Maybe, the ravine along y==-x is badly handled.
|
|
| Mar 18, 2021 at 17:52 | comment | added | Dominic | Sounds so for sure. Would be nice to understand why my code is not exhibiting that. | |
| Mar 18, 2021 at 17:37 | comment | added | user64494 | If I correctly count, the total degree of the numerator equals 4 and the total degree of the denominator equals 12. In view of it I believe that the double integral under consideration converges. | |
| Mar 18, 2021 at 17:20 | history | answered | Dominic | CC BY-SA 4.0 |