- 29,404
- 21,100
PRASHANT KUMAR said:so where does the solution lie?
Okay. Let's look at
##I = \int f(x) \delta((x-a)(x-b))dx## where ##a < b##
Your problem is a special case of this, but it may be easier to look at the general case to see what is going on.
If we draw a graph of this quadratic, we get two roots at ##x = a## and ##x = b##. So, we can see that:
##I = \int_{a-\epsilon}^{a+\epsilon} f(x) \delta((x-a)(x-b))dx + \int_{b-\epsilon}^{b+\epsilon} f(x) \delta((x-a)(x-b))dx##
Where we have focused on the two zeroes and taken a small interval around each, using a small ##\epsilon##
Now, if we do the substitution ##y = (x-a)(x-b)## and ##dy = (2x - a -b) dx## and
##I = \int_{y_1}^{y_2} \frac{f(x)}{2x-a-b} \delta(y)dy + \int_{y_3}^{y_4} \frac{f(x)}{2x-a-b} \delta(y)dy##
Where ##y_1, y_2## etc. are the y-values for ##a-\epsilon, a+\epsilon## etc.
Now, from the graph, we have to note that ##y_1 > 0## and ##y_2 < 0##. So, we need to get these limits in the correct order. Also, as we have this in the Dirac function for ##y##, we are simply evaluating at y=0, which equates to ##x =a## and ##x=b## respectively. So:
##I = -\frac{f(x)}{2x-a-b}|_{x=a} + \frac{f(x)}{2x-a-b}|_{x=b} = \frac{f(a) +f(b)}{b-a}##