"Parameters"
In[]:=
R=0.05;​​L=0.1;​​a=0.0005;
"Ecuation"
"
M0
2
∞
∫
0
​​
J0(qρ)
-q(z)
e
(RJ1(qR)-aJ1(qa))
q
​q"
"Consider
M0
2
=constant"
"J0 and J1 reffers to JBessel"
"Solution"
In[]:=
Integrate[BesselJ[0,qρ]Exp[L-qz](RBesselJ[1,qR]-aBesselJ[1,qa])/q,{q,0,∞}]
Out[]=
∞
∫
0
0.1-qz

BesselJ[0,qρ](-0.0005BesselJ[1,0.0005q]+0.05BesselJ[1,0.05q])
q
q