Discussion about math, puzzles, games and fun. Useful symbols: ÷ × ½ √ ∞ ≠ ≤ ≥ ≈ ⇒ ± ∈ Δ θ ∴ ∑ ∫ π -¹ ² ³ °

You are not logged in.

- Topics: Active | Unanswered

Can anyone determine whether or not the following integral converges?

where is large (but independent of and ), and denotes the standard Euclidean norm on i.e. and denotes the Bessel function of the first kind.I would like for this integral to be at least Equivalently, this integral can be written (less compactly) as:Offline

**bobbym****bumpkin**- From: Bumpkinland
- Registered: 2009-04-12
- Posts: 109,606

Is this symmetric around either axis?

**In mathematics, you don't understand things. You just get used to them.****If it ain't broke, fix it until it is.**** Always satisfy the Prime Directive of getting the right answer above all else.**

Offline

**bobbym****bumpkin**- From: Bumpkinland
- Registered: 2009-04-12
- Posts: 109,606

You say p>>0, any idea about the size of it?

**In mathematics, you don't understand things. You just get used to them.****If it ain't broke, fix it until it is.**** Always satisfy the Prime Directive of getting the right answer above all else.**

Offline

**bobbym****bumpkin**- From: Bumpkinland
- Registered: 2009-04-12
- Posts: 109,606

When we encountered this integral before we were concerned about the many singularities. It is possible that this double integral does not exist because of that. What if that is the case?

**In mathematics, you don't understand things. You just get used to them.****If it ain't broke, fix it until it is.**** Always satisfy the Prime Directive of getting the right answer above all else.**

Offline

There is a singularity at (0,0) and another at (b,c), yes. But in dimensions 3 and higher, it can be shown that the integral converges (the terms have slightly different powers, which depend on d, but the locations of the singularities are the same). To prove convergence for d > 2, we can cut R^d into overlapping pieces -- but we can't do that for d = 2 because we encounter a 1-dimensional integral of r^[(d-3)/2] through r = 0. (And that integral converges iff d > 3.) I can write up the details, but the case d = 2 may be more delicate.

If the double integral diverges, then we will be forced not to estimate the sum with an integral and try to bound the sum in a different way.

*Last edited by zetafunc (2016-11-02 09:42:22)*

Offline

**bobbym****bumpkin**- From: Bumpkinland
- Registered: 2009-04-12
- Posts: 109,606

Since I am not getting anywhere fast with that integral perhaps the sum might be better.

**In mathematics, you don't understand things. You just get used to them.****If it ain't broke, fix it until it is.**** Always satisfy the Prime Directive of getting the right answer above all else.**

Offline

The sum was:

where (\substack does not seem to work on here.)We could also use the asymptotic bound (That was the bound I used for d > 2.) Doing that gives us the integral:or, as a double integral,

*Last edited by zetafunc (2016-11-02 23:02:16)*

Offline

**bobbym****bumpkin**- From: Bumpkinland
- Registered: 2009-04-12
- Posts: 109,606

Hi;

b and c will be a problem for M as before.

**In mathematics, you don't understand things. You just get used to them.****If it ain't broke, fix it until it is.**** Always satisfy the Prime Directive of getting the right answer above all else.**

Offline

**bobbym****bumpkin**- From: Bumpkinland
- Registered: 2009-04-12
- Posts: 109,606

What are b and c again? Reals? Integers?

**In mathematics, you don't understand things. You just get used to them.****If it ain't broke, fix it until it is.**** Always satisfy the Prime Directive of getting the right answer above all else.**

Offline

**bobbym****bumpkin**- From: Bumpkinland
- Registered: 2009-04-12
- Posts: 109,606

So that means an infinite number of singularities.

I need to eat something and that means I have to cook it up first. Sorry, to break off but I am starving see you later.

**In mathematics, you don't understand things. You just get used to them.****If it ain't broke, fix it until it is.**** Always satisfy the Prime Directive of getting the right answer above all else.**

Offline

**bobbym****bumpkin**- From: Bumpkinland
- Registered: 2009-04-12
- Posts: 109,606

The double integral you posted in post #9,

that one does not converge.

**In mathematics, you don't understand things. You just get used to them.****If it ain't broke, fix it until it is.**** Always satisfy the Prime Directive of getting the right answer above all else.**

Offline

**bobbym****bumpkin**- From: Bumpkinland
- Registered: 2009-04-12
- Posts: 109,606

I do not know but M does not think so.

```
ClearAll[b, c];
Integrate[Integrate[1/((x^2 + y^2)^(-3/4) ((b - x)^2 + (c - y)^2)^(-3/4)), {x, -\[Infinity], \[Infinity]}], {y, -\[Infinity], \[Infinity]}]
```

If that is the correct integral then it is possible that M might be wrong but it is at least 10 to 1 against.

**In mathematics, you don't understand things. You just get used to them.****If it ain't broke, fix it until it is.**** Always satisfy the Prime Directive of getting the right answer above all else.**

Offline

**anonimnystefy****Real Member**- From: Harlan's World
- Registered: 2011-05-23
- Posts: 16,039

bobbym wrote:

The double integral you posted in post #9,

that one does not converge.

There shouldn't minuses in those exponents.

*Last edited by anonimnystefy (2016-11-02 20:32:40)*

Here lies the reader who will never open this book. He is forever dead.

Taking a new step, uttering a new word, is what people fear most. ― Fyodor Dostoyevsky, Crime and Punishment

The knowledge of some things as a function of age is a delta function.

Offline

anonimnystefy wrote:

There shouldn't minuses in those exponents.

Thanks for pointing that out -- I've removed them from my post.

bobbym wrote:

I do not know but M does not think so.

`ClearAll[b, c]; Integrate[Integrate[1/((x^2 + y^2)^(-3/4) ((b - x)^2 + (c - y)^2)^(-3/4)), {x, -\[Infinity], \[Infinity]}], {y, -\[Infinity], \[Infinity]}]`

If that is the correct integral then it is possible that M might be wrong but it is at least 10 to 1 against.

I am running it now in M with the negative signs removed.

*Last edited by zetafunc (2016-11-02 20:50:09)*

Offline

**bobbym****bumpkin**- From: Bumpkinland
- Registered: 2009-04-12
- Posts: 109,606

That will make a big difference. Howdy anonimnystefy!

**In mathematics, you don't understand things. You just get used to them.****If it ain't broke, fix it until it is.**** Always satisfy the Prime Directive of getting the right answer above all else.**

Offline

If we get something which looks to be divergent, all is not lost: it just means that we lose too much information by bounding the Bessel functions in modulus. In which case, it may be necessary to use the asymptotics of the Bessel functions, which can get very messy. (But it might be necessary in order to take advantage of the positive-negative cancellation that occurs.)

Offline

**bobbym****bumpkin**- From: Bumpkinland
- Registered: 2009-04-12
- Posts: 109,606

Hi;

The problem is that just like Mathematica, we find that if an expression is too complicated it is hard to do anything analytical to it.

**In mathematics, you don't understand things. You just get used to them.****If it ain't broke, fix it until it is.**** Always satisfy the Prime Directive of getting the right answer above all else.**

Offline

We may need to take the integral at the bottom of post #1, set rho = 1000, say, and sum over all b,c such that b,c are non-zero and b+c=0. If we can find roughly what that should be for some ranges of b,c then we might be able to guess what the answer should be by varying rho.

*Last edited by zetafunc (2016-11-20 22:02:23)*

Offline