i=I;​​ρ=1000;​​v=10^(-6);​​s=20*10^-6;​​m=100^10^-6;​​a=1.02*10^-4;​​​​g=Table(r^2/(i*A))*1-2*BesselJ1,
3
2
i
*r*A
3
2
i
*r*A*BesselJ0,
3
2
i
*r*A,{A,0.4444444444,4,0.03591470258},{r,0.5,1.5,0.0101010101};​​Det[Abs[g]]​​g2=Inverse[g];
Out[]=
0.
---------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
In[]:=
func=
r^2
i*A
*1-
2*BesselJ1,
3
2
i
*r*
A

3
2
i
*r*
A
*BesselJ0,
3
2
i
*r*
A

*
a
s*
2*Pi
*
-
(a^2)*(r-0.981)^2
2*s^2

;
---------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
In[]:=
​​G=Table[NIntegrate[func,{r,0,+Infinity}],{A,1,100}];​​ListPlot[g2*G]
Out[]=
-1.0
-0.5
0.5
1.0
-1.0
-0.5
0.5
1.0