Skip to main content
4 of 4
added 144 characters in body
Fisher
  • 179
  • 6

MinRecursion changes result of integration by 100 orders of magnitude

I want to compute a triple integral. I did it using two codes: the first is NIntegrate[2 1/2 ((50^2 - (((u - v)/2)^2 + t^2 + m^2) Log[ (50^2 + (u - v)^2/4 + t^2 + m^2)/((u - v)^2/4 + t^2 + m^2)])/Log [ ((u + v)^2/4 + t^2 + m^2 + 50^2)/(((u + v)/2)^2 + t^2 + m^2)] ) Exp[-(u)^2/(2 10^-3)], {v, 0, 50 Sqrt[2]},{u, v - 2 50 , -v + 2 50 }, {t, -50, 50}, WorkingPrecision -> 15, MinRecursion -> 8, PrecisionGoal -> 5, AccuracyGoal -> 5]

Where m is a constant equal to 0.1. Notice that the integrand is bounded and smooth in the domain of integration. The only thing is that it has a smooth peak at the u=0 plane.

This gives me no error message, and the result 5675.53769275483. Now i try with the same code, but without specifying MinRecursion, and the result is(again with no error messages) 2.44759996103605*10^(-132) !

How can I rely in such a result?? Anyone knows what is going on?

Also, setting MinRecursion to 4 and varying PrecisionGoal and AccuracyGoal i get different results. E.g. for Precision and Accuracy equal to 8 I get 75853.6468454678, and for Precision=Accu=10 I get 134823.121278000.

Fisher
  • 179
  • 6