Numerical method to solve an integral that contains a singularity?

Click For Summary
The discussion centers on a numerical method for evaluating an integral with a singularity in a condensed matter physics context. The integral in question involves a function that needs to be visualized in terms of the variable ν, ranging from 0 to 1. Initial attempts using Simpson's method from 0 to 1 yielded incorrect results compared to the expected graph. It is noted that there is no universal method for integrating singular integrands, but one effective strategy is to evaluate the integral near the singularity without directly hitting it. Additionally, rewriting the integrand can help isolate the singular part, allowing for more accurate numerical integration.
amjad-sh
Messages
240
Reaction score
13
Homework Statement
//
Relevant Equations
//
In fact I'm working on a condensed matter physics paper, where I stumbled with an integral that I need to visualize.
The function, Ls I need to visualize is equal to:
$$Ls=4\nu^4 \dfrac{\int_{-1}^{1} \dfrac{( 1-u^2)}{(u+\sqrt{u^2-\nu^2})^3} \, du}{\int_{-1}^{1}-u \Big (\dfrac{u-\sqrt{u^2-\nu^2}}{u+\sqrt{u^2+\nu^2}}\Big)^2 \, du}$$

I must draw this function interms of ##\nu##, where ##0\leq\nu \leq1##.
I tried to use Simpson's method by just taking the integral from zero to 1,but it gave me a totally different graph to what is depicted in the paper.
Is there a numerical method that can be used to approximate such integrals?
 
Physics news on Phys.org
Why do you think integrals from 0 to 1 will be the same as integrals from -1 to 1?

Also, your title mentions a singularity. Could you tell us where you think the singularity is?

There is no general approach to numerically integrating a singular integrand. One approach is to simply evaluate the integral close to (but not at) the singularity; the challenge is that as you get closer and closer to the singularity you often need to use more points in your numerical integration rule to maintain accuracy. But sometimes a good first step is to see if cleverness can allow you to rewrite the integrand to allow you to extract the singular part out separately. In your case there is definitely a way to rewrite both of the integrals so that you don't actually need to numerically evaluate something that is singular.
 
  • Like
Likes amjad-sh
Question: A clock's minute hand has length 4 and its hour hand has length 3. What is the distance between the tips at the moment when it is increasing most rapidly?(Putnam Exam Question) Answer: Making assumption that both the hands moves at constant angular velocities, the answer is ## \sqrt{7} .## But don't you think this assumption is somewhat doubtful and wrong?

Similar threads

  • · Replies 7 ·
Replies
7
Views
2K
Replies
3
Views
2K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 7 ·
Replies
7
Views
1K
  • · Replies 105 ·
4
Replies
105
Views
6K
  • · Replies 14 ·
Replies
14
Views
1K
Replies
5
Views
1K
  • · Replies 10 ·
Replies
10
Views
2K
  • · Replies 14 ·
Replies
14
Views
2K
  • · Replies 5 ·
Replies
5
Views
1K