In[]:=
$Version
Out[]=
14.1.0 for Microsoft Windows (64-bit) (July 16, 2024)
In[]:=
ClearAll["`*"];Remove["`*"];
L=200;Ω=1/10;σ=4/10;λ=;Iplus[a_,m_]:=NIntegrate-+π-+,{τ,-Infinity,Infinity};Imins[a_,m_]:=NIntegrate,{τ,-Infinity,Infinity};
1
100
τ
10
25
2
τ
8
1
100
m(-1+Cosh[aτ])
a
Abs[m]Sinh[aτ]
a
2
π
Abs[m]
--+π-+
τ
10
25
2
τ
8
1
100
m(-1+Cosh[aτ])
a
Abs[m]Sinh[aτ]
a
2
π
Abs[m]
In[]:=
M=100;
In[]:=
η0[a_]:=λSum[Iplus[a,m],{m,-M,-1}]+λSum[Iplus[a,m],{m,1,M}];η1[a_]:=λSum[Imins[a,m],{m,-M,-1}]+λSum[Imins[a,m],{m,1,M}];F[a_]:=-(η0[a]^2/(1+η0[a]^2))+1/(1+η1[a]^2);
ListLinePlot[Table[{a,F[a]},{a,1/100,1,1/10}],RotateLabel->False,Frame->True,FrameStyle->Directive[Black,12,Bold],FrameLabel->{Style["a",13,Bold,Black],Style["C",13,Bold,Black]},PlotStyle->{Block,Thickness[0.015],Blue},DataRange->{0,1},PlotRangePadding->0,ImageSize->Medium]
Out[]=
(*-------------------------------------------------------------------------------------------------------------------------------------------------------------------*)
In[]:=
ClearAll["`*"];Remove["`*"];
In[]:=
L=200;Ω=1/10;σ=4/10;λ=;Iplus[a_,m_]:=NIntegrate-+π-+,{τ,-Infinity,Infinity};Imins[a_,m_]:=NIntegrate,{τ,-Infinity,Infinity};
1
100
τ
10
25
2
τ
8
1
100
m(-1+Cosh[aτ])
a
Abs[m]Sinh[aτ]
a
2
π
Abs[m]
--+π-+
τ
10
25
2
τ
8
1
100
m(-1+Cosh[aτ])
a
Abs[m]Sinh[aτ]
a
2
π
Abs[m]
In[]:=
M=200;(*isnotasinfinty*)η0[a_]:=λTotal@Table[Iplus[a,m],{m,-M,-1}]+λTotal@Table[Iplus[a,m],{m,1,M}];η1[a_]:=λTotal@Table[Imins[a,m],{m,-M,-1}]+λTotal@Table[Imins[a,m],{m,1,M}];F[a_]:=-(η0[a]^2/(1+η0[a]^2))+1/(1+η1[a]^2);
In[]:=
ListLinePlot[ParallelTable[{a,F[a]},{a,1/100,1,1/20}],RotateLabel->False,Frame->True,FrameStyle->Directive[Black,12,Bold],FrameLabel->{Style["a",13,Bold,Black],Style["C",13,Bold,Black]},PlotStyle->{Block,Thickness[0.015],Blue},DataRange->{0,1},PlotRangePadding->0,ImageSize->Medium]
Out[]=
(*-------------------------------------------------------------------------------------------------------------------------------------------------------------------*)