N1=1000;
N2=6;dcolor=1/N2;
tao0=6.77/10^5;
V0Eext=7.2 10^6;Epxi0=8.85/10^12;Epxir=6;a1=3.1/10^2;\[Rho]=10^3;R=a1/2;\[CurlyPhi]=(5 \[Pi])/12;omigaw=100 \[Pi];
tlist={0.1,1,10,100,1000};
\[Mu]=1/10^3;
Bc=(Epxi0 (Epxir-1) V0Eext)/(2 a1 Epxir);
\[CapitalDelta]1=Bc Cos[\[CurlyPhi]]-2 tao0;
For[i=1,i<N1+1,{Subscript[k, i]=listBJZs10000[[i]];Subscript[a, i]=(\[Mu] (Subscript[k, i]/R)^2)/\[Rho];Subscript[\[Gamma], i]=ArcTan[(2 omigaw)/Subscript[a, i]]};i++];
VelocityAV[r_]:=\!\(
\*UnderoverscriptBox[\(\[Sum]\), \(n = 1\), \(N1\)]\(
\*FractionBox[\(1\), \(\[Mu]\ \((
\*SubsuperscriptBox[\(k\), \(n\), \(3\)]\
\*SuperscriptBox[\(BesselJ[0,
\*SubscriptBox[\(k\), \(n\)]]\), \(2\)])\)\)]\(\((2\ R)\)\ \((1 - BesselJ[0,
\*SubscriptBox[\(k\), \(n\)]])\)\ BesselJ[1,
\*FractionBox[\(
\*SubscriptBox[\(k\), \(n\)]\ r\), \(R\)]]\ \((\[CapitalDelta]1\ \((1 -
\*SuperscriptBox[\(E\), \(t\ \((\(-
\*SubscriptBox[\(a\), \(n\)]\))\)\)])\) -
\*FractionBox[\(Bc\ \((Cos[\(-
\*SubscriptBox[\(\[Gamma]\), \(n\)]\) + 2\ omigaw\ t + \[CurlyPhi]] -
\*SuperscriptBox[\(E\), \(t\ \((\(-
\*SubscriptBox[\(a\), \(n\)]\))\)\)]\ Cos[\[CurlyPhi] -
\*SubscriptBox[\(\[Gamma]\), \(n\)]])\)\),
SqrtBox[\(
\*SuperscriptBox[\((
\*FractionBox[\(2\ omigaw\),
SubscriptBox[\(a\), \(n\)]])\), \(2\)] + 1\)]])\)\)\)\);
\[CapitalDelta]1
For[j=1,j<N2,{t=tlist[[j]];Subscript[graph, j]=Plot[VelocityAV[r],{r,0,R},Frame->True,FrameLabel->{"r ","\[Omega](r)"},PlotStyle->{RGBColor[dcolor j,1-dcolor j,dcolor j]}]};j++]
Show[Subscript[graph, 1],Subscript[graph, 2],Subscript[graph, 3],Subscript[graph, 4],Subscript[graph, 5]]
N2=6;dcolor=1/N2;
tao0=6.77/10^5;
V0Eext=7.2 10^6;Epxi0=8.85/10^12;Epxir=6;a1=3.1/10^2;\[Rho]=10^3;R=a1/2;\[CurlyPhi]=(5 \[Pi])/12;omigaw=100 \[Pi];
tlist={0.1,1,10,100,1000};
\[Mu]=1/10^3;
Bc=(Epxi0 (Epxir-1) V0Eext)/(2 a1 Epxir);
\[CapitalDelta]1=Bc Cos[\[CurlyPhi]]-2 tao0;
For[i=1,i<N1+1,{Subscript[k, i]=listBJZs10000[[i]];Subscript[a, i]=(\[Mu] (Subscript[k, i]/R)^2)/\[Rho];Subscript[\[Gamma], i]=ArcTan[(2 omigaw)/Subscript[a, i]]};i++];
VelocityAV[r_]:=\!\(
\*UnderoverscriptBox[\(\[Sum]\), \(n = 1\), \(N1\)]\(
\*FractionBox[\(1\), \(\[Mu]\ \((
\*SubsuperscriptBox[\(k\), \(n\), \(3\)]\
\*SuperscriptBox[\(BesselJ[0,
\*SubscriptBox[\(k\), \(n\)]]\), \(2\)])\)\)]\(\((2\ R)\)\ \((1 - BesselJ[0,
\*SubscriptBox[\(k\), \(n\)]])\)\ BesselJ[1,
\*FractionBox[\(
\*SubscriptBox[\(k\), \(n\)]\ r\), \(R\)]]\ \((\[CapitalDelta]1\ \((1 -
\*SuperscriptBox[\(E\), \(t\ \((\(-
\*SubscriptBox[\(a\), \(n\)]\))\)\)])\) -
\*FractionBox[\(Bc\ \((Cos[\(-
\*SubscriptBox[\(\[Gamma]\), \(n\)]\) + 2\ omigaw\ t + \[CurlyPhi]] -
\*SuperscriptBox[\(E\), \(t\ \((\(-
\*SubscriptBox[\(a\), \(n\)]\))\)\)]\ Cos[\[CurlyPhi] -
\*SubscriptBox[\(\[Gamma]\), \(n\)]])\)\),
SqrtBox[\(
\*SuperscriptBox[\((
\*FractionBox[\(2\ omigaw\),
SubscriptBox[\(a\), \(n\)]])\), \(2\)] + 1\)]])\)\)\)\);
\[CapitalDelta]1
For[j=1,j<N2,{t=tlist[[j]];Subscript[graph, j]=Plot[VelocityAV[r],{r,0,R},Frame->True,FrameLabel->{"r ","\[Omega](r)"},PlotStyle->{RGBColor[dcolor j,1-dcolor j,dcolor j]}]};j++]
Show[Subscript[graph, 1],Subscript[graph, 2],Subscript[graph, 3],Subscript[graph, 4],Subscript[graph, 5]]