The idea is to plot in function of “r” “L(r)” by turning “g” into a matrix where “r” varies column wise and “A” line wise and multiplying “G” by the inverse matrix of “g”
G = g*L;
G = g*L;
i=I;ρ=1000;v=10^(-6);s=20*10^-6;m=100^10^-6;a=1.02*10^-4;g=*1-*r**r*;
r^2
i*A
2*BesselJ1,*r*
3
2
i
A
3
2
i
A
*BesselJ0,3
2
i
A
In[]:=
func=*1-*r**r***
r^2
i*A
2*BesselJ1,*r*
3
2
i
A
3
2
i
A
*BesselJ0,3
2
i
A
a
s*
2*Pi
-
(a^2)*(r-0.981)^2
2*s^2
Out[]=
-
(0.+2.03461)1+
-13.005
2
(-0.981+r)
2
r
2BesselJ1,
1/4
(-1)
3/4
(-1)
A
rA
rBesselJ0,3/4
(-1)
A
rA
G In matrix form
G In matrix form
In[]:=
Table[NIntegrate[func,{r,0,+Infinity}],{A,1,100}]//MatrixForm
Out[]//MatrixForm=
0.136994-0.0291638 |
0.120478-0.0498408 |
0.101869-0.0610854 |
0.0850454-0.0657029 |
0.071073-0.0664304 |
0.0598331-0.0650797 |
0.0508695-0.062711 |
0.0437064-0.0599135 |
0.0379414-0.0570049 |
0.0332583-0.0541509 |
0.0294151-0.0514332 |
0.026229-0.0488872 |
0.0235612-0.0465232 |
0.0213063-0.0443383 |
0.0193834-0.0423232 |
0.0177301-0.0404655 |
0.0162976-0.0387519 |
0.0150476-0.0371693 |
0.0139497-0.0357054 |
0.0129795-0.0343487 |
0.0121174-0.0330891 |
0.0113472-0.0319171 |
0.010656-0.0308245 |
0.0100328-0.0298039 |
0.00946868-0.0288487 |
0.00895602-0.0279529 |
0.00848848-0.0271114 |
0.00806068-0.0263195 |
0.00766802-0.0255729 |
0.00730657-0.024868 |
0.00697294-0.0242014 |
0.00666421-0.0235701 |
0.00637785-0.0229714 |
0.00611162-0.0224028 |
0.00586359-0.0218621 |
0.00563206-0.0213474 |
0.00541551-0.0208567 |
0.00521262-0.0203884 |
0.00502221-0.0199411 |
0.00484321-0.0195134 |
0.00467468-0.0191039 |
0.00451579-0.0187116 |
0.00436577-0.0183354 |
0.00422395-0.0179743 |
0.0040897-0.0176274 |
0.00396247-0.0172939 |
0.00384174-0.016973 |
0.00372708-0.016664 |
0.00361804-0.0163662 |
0.00351426-0.0160792 |
0.00341538-0.0158022 |
0.00332108-0.0155347 |
0.00323107-0.0152763 |
0.00314508-0.0150266 |
0.00306287-0.014785 |
0.0029842-0.0145512 |
0.00290886-0.0143249 |
0.00283665-0.0141056 |
0.00276741-0.013893 |
0.00270095-0.0136868 |
0.00263713-0.0134868 |
0.00257579-0.0132927 |
0.00251682-0.0131042 |
0.00246007-0.012921 |
0.00240543-0.012743 |
0.0023528-0.0125699 |
0.00230207-0.0124015 |
0.00225315-0.0122377 |
0.00220595-0.0120782 |
0.00216038-0.0119229 |
0.00211637-0.0117716 |
0.00207384-0.0116241 |
0.00203272-0.0114804 |
0.00199295-0.0113402 |
0.00195447-0.0112035 |
0.00191722-0.0110701 |
0.00188114-0.0109399 |
0.00184619-0.0108128 |
0.00181231-0.0106886 |
0.00177946-0.0105673 |
0.00174759-0.0104488 |
0.00171667-0.0103329 |
0.00168666-0.0102196 |
0.00165751-0.0101089 |
0.0016292-0.0100005 |
0.00160169-0.00989449 |
0.00157495-0.00979073 |
0.00154895-0.00968915 |
0.00152366-0.0095897 |
0.00149905-0.0094923 |
0.0014751-0.0093969 |
0.00145179-0.00930342 |
0.00142908-0.00921181 |
0.00140697-0.00912202 |
0.00138542-0.00903399 |
0.00136442-0.00894767 |
0.00134395-0.00886301 |
0.00132398-0.00877996 |
0.00130451-0.00869848 |
0.00128552-0.00861852 |