XX model at finite temperature

Nave[h_,T_]:=(1/(2π))NIntegrate[(1/(Exp[(-h-Cos[k])/T]+1)),{k,-π,π},Method->{Automatic,"SymbolicProcessing"->False}];
Have[h_,T_]:=(1/(2π))NIntegrate[(((-h - Cos[k]))/(Exp[(-h-Cos[k])/T]+1)),{k,-π,π},Method->{Automatic,"SymbolicProcessing"->False}];
mag[h_,T_]:=(1/(2π))NIntegrate[Tanh[((h+Cos[k])/(2T))],{k,-π,π},Method->{Automatic,"SymbolicProcessing"->False}];
χ[h_,T_]:=(1/(2π))(1/(2T))NIntegrate[((Sech[((h+Cos[k])/(2T))])^2),{k,-π,π},Method->{Automatic,"SymbolicProcessing"->False}];
c[h_,T_]:=(1/(2π))(1/(2T))NIntegrate[(((((-h - Cos[k])/T))^2)/(((Exp[(-h-Cos[k])/T]+1))^2))Exp[(-h-Cos[k])/T],{k,-π,π},Method->{Automatic,"SymbolicProcessing"->False}];

Show[
Quiet@ListLinePlot[Table[{h,mag[h,T]},{T,{0.05,0.1,0.2}},{h,-2,2,0.05}],
PlotLegends->LineLegend[Automatic,{"T = 0.05", "T = 0.1", "T = 0.2"}]],
Plot[(2/π)ArcCos[-h]-1,{h,-1,1},PlotStyle->Directive[Black,Dashed]],
FrameLabel->{"h","m"}
]

Show[
Quiet@ListLinePlot[Table[{h,χ[h,T]},{T,{0.05,0.1,0.2}},{h,-2,2,0.03}],
PlotLegends->LineLegend[Automatic,{"T = 0.05", "T = 0.1", "T = 0.2"}]],
Plot[(2/π)(1/(Sqrt[1-(h^2)])),{h,-1,1},PlotStyle->Directive[Black,Dashed]],
FrameLabel->{"h","χ"}
]

Show[
Quiet@ListLinePlot[Table[{h,Have[h,T]},{T,{0.02,0.1,0.2}},{h,-2,2,0.03}],
PlotLegends->LineLegend[Automatic,{"T = 0.05", "T = 0.1", "T = 0.2"}]],
FrameLabel->{"h","⟨H⟩"}
]

Show[
Quiet@ListLinePlot[Table[{T,Have[h,T]},{h,{0.5,0.95,1.0,1.05,1.5}},{T,0.02,2.,0.2}],
PlotLegends->LineLegend[Automatic,{"h = 0.5", "h= 1", "h = 1.5"}]],
FrameLabel->{"T","⟨H⟩"}
]

Show[
Quiet@ListLinePlot[Table[{T,c[h,T]},{h,{0.5,1.0,1.5}},{T,0.02,1.,0.01}],
PlotLegends->LineLegend[Automatic,{"h = 0.5", "h = 1" ,"h = 1.5"}]],
FrameLabel->{"T","c"}
]

Show[
Quiet@ListLinePlot[Table[{h,c[h,T]},{T,{0.05,0.1,0.2}},{h,-2,2,0.03}],
PlotLegends->LineLegend[Automatic,{"T = 0.05", "T = 0.1", "T = 0.2"}]],
Plot[(2/π)(1/(Sqrt[1-(h^2)])),{h,-1,1},PlotStyle->Directive[Black,Dashed]],
FrameLabel->{"h","χ"}
]