In[]:=
$Version
Out[]=
14.0.0 for Microsoft Windows (64-bit) (December 13, 2023)
In[]:=
ClearAll["`*"];Remove["`*"];
In[]:=
func=-//FullSimplify
1+
-
2
(-m+r)
2
2
s
π
2
2
r
2BesselJ1,r
1/4
(-1)
v
A
3/4
(-1)
v
A
rBesselJ0,r
3/4
(-1)
v
A
Asρ
Out[]=
-
2
(m-r)
2
2
s
π
2
4
r
A
2
r
4v
4svρHypergeometric0F1Regularized1,
A
2
r
4v
In[]:=
IntegrateChangeVariables[Inactive[Integrate][func,{r,0,Infinity}],t,t10^5*r]//ExpandAll(*Rescaling*)
Out[]=
∞
∫
0
-+-
2
m
2
2
s
mt
100000
2
s
2
t
20000000000
2
s
π
2
4
t
A
2
t
40000000000v
40000000000000000000000000svρHypergeometric0F1Regularized1,
A
2
t
40000000000v
In[]:=
ρ=1000;v=10^(-6);m=10^(-4);s=2*10^(-5);newfunc=Hypergeometric0F1Regularized3,;
-+-
2
m
2
2
s
mt
100000
2
s
2
t
20000000000
2
s
π
2
4
t
A
2
t
40000000000v
40000000000000000000000000svρHypergeometric0F1Regularized1,
A
2
t
40000000000v
In[]:=
LogLogPlot[Abs[NIntegrate[newfunc,{t,0,Infinity}]],{A,10^-10,10^10},PlotRangeAll]
Out[]=
