Double integral of a Log (natural)

Click For Summary

Homework Help Overview

The discussion revolves around evaluating a double integral of the logarithmic function over a rectangular region defined by variables u and v. The integral involves parameters a and b, which are positive constants, and the variables x and y, which are real numbers.

Discussion Character

  • Exploratory, Assumption checking, Mathematical reasoning

Approaches and Questions Raised

  • Participants discuss the differences in results when integrating with respect to u first versus v first, noting discrepancies in the outcomes. Some participants question the conditions under which the integral may diverge, particularly when the logarithm's argument approaches zero. Others explore the implications of the arctangent function in the context of the integral.

Discussion Status

The conversation is ongoing, with participants sharing their findings and questioning the behavior of the integral under various conditions. Some have provided insights into potential divergence and the behavior of the integral at the boundaries, while others are seeking clarification on how to achieve consistent results regardless of the order of integration.

Contextual Notes

There are discussions about the integral's behavior within the specified boundaries and the implications of approaching those boundaries. Participants are also considering the effects of specific values of x and y on the integral's convergence.

Mbert
Messages
63
Reaction score
0

Homework Statement


I wish to find the following integral over the rectangle [-a,a] in u and [-b,b] in v using Mathematica. The constants a and b are positive (and non-zero). The variables x and y are in Reals.

[tex]A(x,y)=\int_{-a}^a{\int_{-b}^b{\log{\left[(u-x)^2+(v-y)^2\right]}{\rm d}v}{\rm d}u}[/tex]

Homework Equations



see above

The Attempt at a Solution



I find using Mathematica a different integral depending if I integrate on u or v first. I use the following code:

To integrate on v first:
Code:
Integrate[Log[(u - x)^2 + (v - y)^2], v]
(% /. v -> b) - (% /. v -> -b)
Integrate[%, u]
(% /. u -> a) - (% /. u -> -a)

To integrate on u first:
Code:
Integrate[Log[(u - x)^2 + (v - y)^2], u]
(% /. u -> a) - (% /. u -> -a)
Integrate[%, v]
(% /. v -> b) - (% /. v -> -b)

The difference between these two (definite) double integrals is:
[tex](x-y) (x+y) \left(\text{ArcTan}\left[\frac{a-x}{b-y}\right]+\text{ArcTan}\left[\frac{a+x}{b-y}\right]+\text{ArcTan}\left[\frac{b-y}{a-x}\right]+\text{ArcTan}\left[\frac{b-y}{a+x}\right]+\text{ArcTan}\left[\frac{a-x}{b+y}\right]+\text{ArcTan}\left[\frac{a+x}{b+y}\right]+\text{ArcTan}\left[\frac{b+y}{a-x}\right]+\text{ArcTan}\left[\frac{b+y}{a+x}\right]\right)[/tex]

Could someone please tell me what I am doing wrong?
 
Physics news on Phys.org
You know, there is a relation
[tex]\arctan(\frac{1}{x}) = \frac{\pi}{2} sgn(x) - \arctan(x)[/tex]
By looking at your fractions, you will see that there is always one that is reciprocal to another. So I'm guessing that under certain conditions, you will get 0.

I'm saying under certain conditions, because I suspect the integral to diverge for -a<y<a and -b<x<b and then you can't expect to get 0 anymore, but it's just a guessing, I could be wrong, you should apply the relation and see if and when you get 0.
 
Thank you.

Wouldn't it diverge for y=+ or - b and x=+ or - a? Why would it diverge for values in between? I think this should be defined up to the boundary and the term sgn(x) depends on which side of the boundary we are approaching the limit.
 
Well, my idea was just that, since
[tex]\log 0 = -\infty[/tex]
the integral might diverge in case the argument of the logarithm becomes zero. And as long as [itex]-a \leq y \leq a[/itex] and [itex]-b \leq x \leq b[/itex], you have this zero in your logarithm. Anyhow, I'm not saying that it is definitely divergent. If you want to be sure, you should do the integration manually. Cheers
 
It seems to give 2*pi inside (−b<y<b and −a<x<a) and Indeterminate on the boundary. Outside the boundary, it gives 0.
 
Does that mean that if I am interested in the domain (−b<y<b and −a<x<a), I can just add the difference between the 2 definite double integrals (which is 2*pi*(x−y)*(x+y))? But then, if I integrate the other double integral first, the difference between them is 2*pi*(y−x)*(y+x)...

How can I make the definite double integral give the same result for (−b<y<b and −a<x<a) whether I integrate u or v first?
 
Well, it certainly seems like there is something fishy if x,y are inside the boundary.
Can you check if you get a finite result for x,y inside the boundary?
 
The result of both integrals inside the boundary is finite. The result inside the boundary of the difference between them is 2*pi*(x−y)*(x+y) or 2*pi*(y−x)*(y+x), depending on which is integrated first.
 
Can you post what you got for one of the integrals?
 
  • #10
[tex]\begin{align*}A(x,y,a,b)=-12 a b &+(a-x) (b-y) \text{Log}\left[(a-x)^2+(b-y)^2\right]\\ &+(a+x) (b-y) \text{Log}\left[(a+x)^2+(b-y)^2\right]\\ &+(a-x) (b+y) \text{Log}\left[(a-x)^2+(b+y)^2\right]\\ &+(a+x) (b+y) \text{Log}\left[(a+x)^2+(b+y)^2\right]\\ &+(a-x)^2 \left(-\text{ArcTan}\left[\frac{a-x}{b-y}\right]+\text{ArcTan}\left[\frac{-a+x}{b+y}\right]\right)\\ &+(a+x)^2 \left(-\text{ArcTan}\left[\frac{a+x}{b-y}\right]-\text{ArcTan}\left[\frac{a+x}{b+y}\right]\right)\\ &-(b-y)^2 \left(\text{ArcTan}\left[\frac{b-y}{a-x}\right]+\text{ArcTan}\left[\frac{b-y}{a+x}\right]\right)\\ &-(b+y)^2 \left(\text{ArcTan}\left[\frac{b+y}{a-x}\right]+\text{ArcTan}\left[\frac{b+y}{a+x}\right]\right)\\ &+x^2\left(-\text{ArcTan}\left[\frac{a-x}{b-y}\right]-\text{ArcTan}\left[\frac{a+x}{b-y}\right]-\text{ArcTan}\left[\frac{b-y}{a-x}\right]-\text{ArcTan}\left[\frac{b-y}{a+x}\right]+\text{ArcTan}\left[\frac{-a+x}{b+y}\right]-\text{ArcTan}\left[\frac{a+x}{b+y}\right]-\text{ArcTan}\left[\frac{b+y}{a-x}\right]-\text{ArcTan}\left[\frac{b+y}{a+x}\right]\right)\end{align*}[/tex]
 
  • #11
If I integrate with respect to the other variable first, the last term is multiplied by y^2 instead of x^2.
 

Similar threads

Replies
10
Views
2K
  • · Replies 11 ·
Replies
11
Views
3K
Replies
4
Views
3K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 7 ·
Replies
7
Views
3K
  • · Replies 9 ·
Replies
9
Views
3K
  • · Replies 2 ·
Replies
2
Views
2K
Replies
20
Views
2K
  • · Replies 5 ·
Replies
5
Views
3K
  • · Replies 5 ·
Replies
5
Views
3K