Needs["ErrorBarLogPlots`"]
PantheonAll=Import["Documents/Hubbletension/pantheon.txt","Table"];
PantheonAll[[2;;Length@PantheonAll]];
listErrorPTdm=Table[{%[[i,2]],%[[i,5]],%[[i,6]]},{i,1,Length@%,1}];
PantheonAll[[2;;Length@PantheonAll]];
listErrorPTdL=Table[{%[[i,3]],(10^((1/5)(-25+(%[[i,5]]+19.3))))(1/1000),(10^((1/5)(-25+(%[[i,6]]+19.3))))(1/1000)},{i,1,Length@%,1}];
obserrordm=ErrorListLogLinearPlot[listErrorPTdm,PlotStyle->RGBColor[0,0,0],Frame->True,FrameLabel->TraditionalForm/@{"z","(Subscript[m, B])"},BaseStyle->{FontFamily->"Times",FontSize->16},PlotRange->All]
obserrordL=ErrorListPlot[listErrorPTdL,PlotStyle->RGBColor[0,0,0],Frame->True,FrameLabel->TraditionalForm/@{"z","(Subscript[d, L]) (Gpc)"},BaseStyle->{FontFamily->"Times",FontSize->16},PlotRange->All]

paramsub={c->3  (10^8),H0->73,Omegar0->0.00009204227740000001,OmegaDE0->0.692,Omegak0H02->0};
paramsubn1={c->3  (10^8),H0->73,Omegar0->0.00009204227740000001,OmegaDE0->0.692,Omegak0H02->1};
paramsubp1={c->3  (10^8),H0->73,Omegar0->0.00009204227740000001,OmegaDE0->0.692,Omegak0H02->-1};
paramsubt={c->3  (10^8),H0->73,Omegar0->0.00009204227740000001,OmegaDE0->0,Omegak0H02->0};

dLΛCDM=ListPlot[Table[{i,(1+i)((c/.paramsub)/(1000*1000H0/.paramsub)) NIntegrate[1/(√((1-OmegaDE0-Omegar0-(Omegak0H02/(H0^2)))(((1+z))^3)+(OmegaDE0)+(Omegar0)(((1+z))^4)+((Omegak0H02/(H0^2))) (((1+z))^2))/.paramsub),{z,0,i}]},{i,0.,2.3,0.01}],Frame->True(*,PlotLabel->"Evolution of density parameters"*),FrameLabel->TraditionalForm/@{"z","(Subscript[d, L]) (Gpc)"},BaseStyle->{FontFamily->"Times",FontSize->16},PlotLegends->{"ΛCDM"},PlotRange->All(*,PlotRange->{{0,3.1},{60,100}}*),Joined->True]

dmΛCDM=ListLogLinearPlot[Table[{i,5Log10[(1+i)((c/.paramsub)/(1000H0/.paramsub))NIntegrate[1/(√((1-OmegaDE0-Omegar0-(Omegak0H02/(H0^2)))(((1+z))^3)+(OmegaDE0)+(Omegar0)(((1+z))^4)+((Omegak0H02/(H0^2))) (((1+z))^2)/.paramsub)),{z,0,i}]]+25-19.328},{i,0.,2.3,0.01}],Frame->True(*,PlotLabel->"Evolution of density parameters"*),FrameLabel->TraditionalForm/@{"z","(Subscript[m, B])"},BaseStyle->{FontFamily->"Times",FontSize->16},PlotLegends->{"ΛCDM"},PlotRange->All(*,PlotRange->{{0,3.1},{60,100}}*),Joined->True,PlotStyle->{Red}]

dmΛCDMtrial=ListLogLinearPlot[Table[{i,5Log10[(1+i)((c/.paramsubt)/(1000H0/.paramsubt))NIntegrate[1/(√((1-OmegaDE0-Omegar0-(Omegak0H02/(((1000 H0))^2)))(((1+z))^3)+(OmegaDE0)+(Omegar0)(((1+z))^4)+((Omegak0H02/(((1000 H0))^2)))(((1+z))^2))/.paramsubt),{z,0,i}]]+25-19.328},{i,0.,2.3,0.01}],Frame->True(*,PlotLabel->"Evolution of density parameters"*),FrameLabel->TraditionalForm/@{"z","(Subscript[m, B])"},BaseStyle->{FontFamily->"Times",FontSize->16},PlotLegends->{"ΛCDM0DE"},PlotRange->All(*,PlotRange->{{0,3.1},{60,100}}*),Joined->True,PlotStyle->{Blue}]

dmΛCDMn1=ListLogLinearPlot[Table[{i,5Log10[(1+i)((c/.paramsubn1)/(1000H0/.paramsubn1))Sinh[NIntegrate[1/(√((1-OmegaDE0-Omegar0-(Omegak0H02/(((1000H0))^2)))(((1+z))^3)+(OmegaDE0)+(Omegar0)(((1+z))^4)+((Omegak0H02/(((1000 H0))^2)))(((1+z))^2))/.paramsubn1),{z,0,i}]]]+25-19.328},{i,0.,2.3,0.01}],Frame->True(*,PlotLabel->"Evolution of density parameters"*),FrameLabel->TraditionalForm/@{"z","(Subscript[m, B])"},BaseStyle->{FontFamily->"Times",FontSize->16},PlotLegends->{"ΛCDMk-1"},PlotRange->All(*,PlotRange->{{0,3.1},{60,100}}*),Joined->True,PlotStyle->{Dashed,Green}]

dmΛCDMp1=ListLogLinearPlot[Table[{i,5Log10[(1+i)((c/.paramsubp1)/(1000H0/.paramsubp1))Sin[NIntegrate[1/(√((1-OmegaDE0-Omegar0-(Omegak0H02/(((1000 H0))^2)))(((1+z))^3)+(OmegaDE0)+(Omegar0)(((1+z))^4)+((Omegak0H02/(((1000 H0))^2)))(((1+z))^2))/.paramsubp1),{z,0,i}]]]+25-19.328},{i,0.,2.3,0.01}],Frame->True(*,PlotLabel->"Evolution of density parameters"*),FrameLabel->TraditionalForm/@{"z","(Subscript[m, B])"},BaseStyle->{FontFamily->"Times",FontSize->16},PlotLegends->{"ΛCDMk1"},PlotRange->All(*,PlotRange->{{0,3.1},{60,100}}*),Joined->True,PlotStyle->{Dashed,Orange}]

Show[obserrordm,dmΛCDM,dmΛCDMn1,dmΛCDMp1,dmΛCDMtrial]