Friday, 8 September 2017

integration - Having trouble evaluting error function integrals



I am trying to evaluate
I=1(erf(ablog(x))2x2erf(a+blog(x))2x)dx



Let log(x)=v gives x=ev and dx=evdv



I=2π0evabv0eu2dudv2π0a+bv0eu2dudv

I am not sure how to proceed from here though. Any help would be greatly appreciated.


Answer



Let me consider the two integrals I=evabv0eu2dudv J=a+bv0eu2dudv First abv0eu2du=π2erf(abv) which makes I=π2everf(abv)dv This one can be integrated by parts and, after completing the square, we get I=π2(e14ab4b2erf(abv12b)everf(abv)) In a similar manner a+bv0eu2du=π2erf(a+bv) which makes J=π2erf(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)=11πez2(1z12z3+O(1z4))


No comments:

Post a Comment

real analysis - How to find limhrightarrow0fracsin(ha)h

How to find lim without lhopital rule? I know when I use lhopital I easy get $$ \lim_{h\rightarrow 0}...