I am trying to evaluate
I=∫∞1(erf(a−blog(x))2x2−erf(a+blog(x))2x)dx
Let log(x)=v gives x=ev and dx=evdv
I=2√π∫∞0e−v∫a−bv0e−u2dudv−2√π∫∞0∫a+bv0e−u2dudv
I am not sure how to proceed from here though. Any help would be greatly appreciated.
Answer
Let me consider the two integrals I=∫e−v∫a−bv0e−u2dudv J=∫∫a+bv0e−u2dudv First ∫a−bv0e−u2du=√π2erf(a−bv) which makes I=√π2∫e−verf(a−bv)dv This one can be integrated by parts and, after completing the square, we get I=√π2(e1−4ab4b2erf(a−bv−12b)−e−verf(a−bv)) In a similar manner ∫a+bv0e−u2du=√π2erf(a+bv) which makes J=√π2∫erf(a+bv)dv Again, an obvious change of variable and integration by parts lead to J=√π2((ab+v)erf(a+bv)+e−(a+bv)2√πb) There is no problem for the lower bound (v=0); for the upper bound, consider v=t and I suppose that you need to look at the asymptotics for the erf(z) function when z→±∞ (have a look here). For infinitely large values of z erf(z)=1−1√πe−z2(1z−12z3+O(1z4))
No comments:
Post a Comment