First, cut off the tail towards infinity:
f(x)=int1xfracPhi(t)t5dt+intinfty1fracPhi(t)t5dt.
The second term is a constant, so you can compute it numerically once and for all.
Write
eipiu2/2=1+fracipi2u2+R(u)
and
intt0eipiu2/2du=t+fracipi6t3+intt0R(u)du.
So
fracPhi(t)t5=left(t−4+fracipi6t−2+t−5intt0R(u)duright)left(1+fracipi2t2+R(t)right)=
t−4+frac2pii3t−2+left(t−4R(t)−fracpi212+intt0R(u)duright).
So
int1xfracPhi(t)t5dt=frac13left(x−3−1right)+frac2pii3left(x−1−1right)+int1xleft(t−4R(t)−fracpi212+intt0R(u)duright)du.
The integrands in the last term are bounded functions, and they are being integrated over bounded domains, so there is no problem approximating them numerically.
If you want an asymptotic formula, instead of a numerical approximation, you should be able to keep taking more terms out to get a formula like
f(x)=frac13x−3+frac2pii3x−1+C+a1x+a2x2+cdots+anxn+O(xn+1)quadmathrmasxto0.
You probably won't be able to get the constant C in closed form, because it involves all those convergent integrals. The other ai will be gettable in closed form, although they will get worse and worse as you compute more of them.
No comments:
Post a Comment