sp6Flux0to000005=Import["sp6Flux0to000005.dat"];

sp6Flux0to000005v3=Delete[sp6Flux0to000005,1];

Ω[T_,i_,μ_]:=-(T/(2*169*169)) Sum[Log[1+Exp[-(1/T)(sp6Flux0to000005v3[[i,j]]-μ)]],{j,2,47322}]
χ[T_,μ_]:=Module[{tmp,fit},
tmp=Table[{(i-1)/100000,Ω[T,i,μ]},{i,1,6}];
fit=Fit[tmp,{1,x,(x^2)},x];
-2Coefficient[fit,(x^2)]]

δμ=0.01;

temp=0.1;

res=Monitor[Table[{μ,χ[temp,μ]},{μ,-5,+5,δμ}],μ];

ListLinePlot[res,PlotRange->All]

ListLinePlot[{res,chiv17},PlotRange->All]

chiv17=Import["chiv17.dat"];

ListLinePlot[chiv16,PlotRange->All]