Tensor part of the stress

kp[t_]=1(((t))^(-1/2));
magspec[k_,t_]=(1/150)((((k/kp[t]))^4)/(1+(((k/kp[t]))^(4+5/3))))(t^(-1/2));
stress[k_,t_]:=NIntegrate[magspec[k1,t]magspec[√((k^2)+(k1^2)-2k k1 μ),t]((k^2)/((k^2)+(k1^2)-2k k1 μ))(1+(μ^2))(1+(((k-k1 μ))^2)/((k^2)+(k1^2)-2k k1 μ))(2π/4π),{k1,0.001,(10^5)},{μ,-1,1}];

kph[t_]=1(((t))^(-2/3));
magspech[k_,t_]=(1/150)((((k/kph[t]))^4)/(1+(((k/kph[t]))^(4+5/3))));
stressh[k_,t_]:=NIntegrate[magspech[k1,t]magspech[√((k^2)+(k1^2)-2k k1 μ),t]((k^2)/((k^2)+(k1^2)-2k k1 μ))((1+(μ^2))(1+(((k-k1 μ))^2)/((k^2)+(k1^2)-2k k1 μ))+4μ((k-k1 μ)/((((k^2)+(k1^2)-2k k1 μ))^(1/2))))(2π/4π),{k1,0.001,(10^5)},{μ,-1,1}];

findLogDivisions[{xmin_,xmax_},n_Integer]:=10^FindDivisions[Log10@{xmin,xmax},n]
interval=findLogDivisions[{0.01,100},100];
dd=Dimensions[interval];

plots1=Table[{interval[[i]],stress[interval[[i]],1]},{i,1,dd[[1]]}];
plots2=Table[{interval[[i]],stress[interval[[i]],11.1111]},{i,1,dd[[1]]}];
plots3=Table[{interval[[i]],stress[interval[[i]],100]},{i,1,dd[[1]]}];
plotm1=Table[{interval[[i]],magspec[interval[[i]],1]},{i,1,dd[[1]]}];
plotm2=Table[{interval[[i]],magspec[interval[[i]],11.1111]},{i,1,dd[[1]]}];
plotm3=Table[{interval[[i]],magspec[interval[[i]],100]},{i,1,dd[[1]]}];

plotsh1=Table[{interval[[i]],stressh[interval[[i]],1]},{i,1,dd[[1]]}];
plotsh2=Table[{interval[[i]],stressh[interval[[i]],6.08581]},{i,1,dd[[1]]}];
plotsh3=Table[{interval[[i]],stressh[interval[[i]],31.6228]},{i,1,dd[[1]]}];
plotmh1=Table[{interval[[i]],magspech[interval[[i]],1]},{i,1,dd[[1]]}];
plotmh2=Table[{interval[[i]],magspech[interval[[i]],6.08581]},{i,1,dd[[1]]}];
plotmh3=Table[{interval[[i]],magspech[interval[[i]],31.6228]},{i,1,dd[[1]]}];

Figure 7 and 8

SetDirectory["/home/ramkishor/Desktop/GWspectra_from_mhd_turbulence (Version 1229)/mathematica/figures"];
plot1=ListLogLogPlot[{plotm1,plotm2,plotm3,plots1,plots2,plots3},Joined->True,PlotRange->{{0.005,100},{(10^(-10)),(10^(-1))}},PlotStyle->{Blue,{Blue,Dotted},{Blue ,Dashed},Red,{Red,Dotted},{Red ,Dashed},Black,{Dashed,Black},{Dashed,Blue},{Dashed,Black}},Frame->True,FrameStyle->Directive[Black,Bold,12],FrameLabel->{Style["k",Black,Bold,14],Style["(Subscript[E, M])(k) and Sp(T)",Black,12]},PlotLegends->Placed[LineLegend[{Style["(Subscript[k, p])=1",10],Style["(Subscript[k, p])=0.3",10],Style["(Subscript[k, p])=0.1",10]},LegendFunction->"Frame",LegendMarkerSize->12],{0.4,0.17}]];
plot2=ListLogLogPlot[{plots1,{plots2[[All,1]],((((1/(11.1111^(1/2)))))^(-1))plots2[[All,2]]}//Transpose,{plots3[[All,1]],((((1/(100^(1/2)))))^(-1))plots3[[All,2]]}//Transpose,plots1,{plots2[[All,1]],((((1/(11.1111^(1/2)))))^(-14/3))plots2[[All,2]]}//Transpose,{plots3[[All,1]],((((1/(100^(1/2)))))^(-14/3))plots3[[All,2]]}//Transpose},Joined->True,PlotRange->{{0.005,100},{(10^(-10)),(10^(-1))}},PlotStyle->{Blue,{Blue,Dotted},{Blue ,Dashed},{Red,Opacity[0.5]},{Red,Dotted},{Red ,Dashed},Black,{Dashed,Black},{Dashed,Blue},{Dashed,Black}},Frame->True,FrameStyle->Directive[Black,Bold,12],FrameLabel->{Style["k",Black,Bold,14],Style["Sp(T)/(Subscript[k, p])and Sp(T)/((Subscript[k, p])^(14/3))",Black,12]},PlotLegends->Placed[LineLegend[{Style["(Subscript[k, p])=1",10],Style["(Subscript[k, p])=0.3",10],Style["(Subscript[k, p])=0.1",10]},LegendFunction->"Frame",LegendMarkerSize->12],{0.42,0.25}]];
plot3=ListLogLogPlot[{plotmh1,plotmh2,plotmh3,plotsh1,plotsh2,plotsh3},Joined->True,PlotRange->{{0.005,100},{(10^(-10)),(10^(-1))}},PlotStyle->{Blue,{Blue,Dotted},{Blue ,Dashed},Red,{Red,Dotted},{Red ,Dashed},Black,{Dashed,Black},{Dashed,Blue},{Dashed,Black}},Frame->True,FrameStyle->Directive[Black,Bold,12],FrameLabel->{Style["k",Black,Bold,14],Style["(Subscript[E, M])(k) and Sp(T)",Black,12]},PlotLegends->Placed[LineLegend[{Style["(Subscript[k, p])=1",10],Style["(Subscript[k, p])=0.3",10],Style["(Subscript[k, p])=0.1",10]},LegendFunction->"Frame",LegendMarkerSize->12],{0.48,0.25}]];
plot4=ListLogLogPlot[{plotsh1,{plotsh2[[All,1]],((1/(6.08581^(2/3))))plotsh2[[All,2]]}//Transpose,{plotsh3[[All,1]],((1/(31.6228^(2/3))))plotsh3[[All,2]]}//Transpose,plots1,{plotsh2[[All,1]],((((1/(6.08581^(2/3)))))^(-8/3))plotsh2[[All,2]]}//Transpose,{plotsh3[[All,1]],((((1/(31.6228^(2/3)))))^(-8/3))plotsh3[[All,2]]}//Transpose},Joined->True,PlotRange->{{0.005,100},{(10^(-10)),(10^(-1))}},PlotStyle->{Blue,{Blue,Dotted},{Blue ,Dashed},Red,{Red,Dotted},{Red ,Dashed},Black,{Dashed,Black},{Dashed,Blue},{Dashed,Black}},Frame->True,FrameStyle->Directive[Black,Bold,12],FrameLabel->{Style["k",Black,Bold,14],Style["(Subscript[k, p])Sp(T) and Sp(T)/((Subscript[k, p])^(8/3))" ,Black,12]},PlotLegends->Placed[LineLegend[{Style["(Subscript[k, p])=1",10],Style["(Subscript[k, p])=0.3",10],Style["(Subscript[k, p])=0.1",10]},LegendFunction->"Frame",LegendMarkerSize->12],{0.42,0.25}]];
Export["nonhelical1.eps",plot1];
Export["nonhelical2.eps",plot2];
Export["helical1.eps",plot3];
Export["helical2.eps",plot4];

Figure 6

plotnn=ListLogLogPlot[{plotm1,plotmh1,plots1,plotsh1},Joined->True,PlotStyle->{{Blue,Dashed},{Red,Dotted},Blue,Red,Black,{Dashed,Black},{Dashed,Blue},{Dashed,Black}},Frame->True,FrameStyle->Directive[Black,Bold,12],FrameLabel->{Style["k",Black,Bold,14],Style["(Subscript[E, M])(k) and Sp(T)",Black,Bold,14]}];
Export["helical_vs_nonhelical.eps",plotnn];