ClearAll["Global`*"]
tmax1[γf_,γc_]:=2/(γc*(2-γf)^2)
tmax2[γf_,γc_]:=(γc*(γf-1)^2+2*(γf-γc)-2*(γf-γc)*Sqrt[γf*(2-γf)])/(γf^2 *(2-γf)^2 /2)
q[γf_]:=(2-γf)/2
F1[γf_,γc_,t_]:=q[γf]*(1+γf*q[γf]*t/γc)^2+1/γc*(1-γf/γc)^2*q[γf]*t+(1+γf*q[γf]*t/γc)*(1/γc - γf/(γc*γc))
F2[γf_,γc_,t_]:=1+γf*q[γf]*t/γc
F3[γf_,γc_,t_]:=1/γc*(1-γf/γc)^2*(2*q[γf]*t+2*γf*q[γf]^2*t^2/γc)
F4[γf_,γc_,t_]:=1/(γc*γc)*(1-γf/γc)^2*(1+q[γf]*t)^2
F5[γf_,γc_,t_]:=FullSimplify[(F3[γf,γc,t]-2*F1[γf,γc,t]*F2[γf,γc,t] +Sign[1-γf/γc]*Sqrt[(2*F1[γf,γc,t]*F2[γf,γc,t]-F3[γf,γc,t])^2+4*F2[γf,γc,t]^2*(F4[γf,γc,t]-F1[γf,γc,t]^2)])/(2*F2[γf,γc,t]^2),Assumptions->{0<γf<1&&γc>0 && {γf,γc,t} ∈ℝ && -(-2+γc (((-2+γf))^2) t)>=0}]
F5[γf,γc,t]
rho[γf_,γc_,t_]:=FullSimplify[1/γf*(1-Sqrt[1+2*γf*F5[γf,γc,t]]), Assumptions->{0<γf<1&&γc>0 && {γf,γc,t} ∈ℝ && 1+2*γf*F5[γf,γc,t]>=0}]
rho[γf,γc,t]
TeXForm[rho[γf,γc,t]]

((2 ((γc)^2) (-2+γf)+4 γf+2 (γc-γf) (Sqrt[4-2 γc (((-2+γf))^2) t])+γc (-4-(((-2+γf))^2) γf t))/(4 ((γc)^2)))

((1-(Sqrt[1+((γf (2 ((γc)^2) (-2+γf)+4 γf+2 (γc-γf) (Sqrt[4-2 γc (((-2+γf))^2) t])+γc (-4-(((-2+γf))^2) γf t)))/(2 ((γc)^2)))]))/(γf))

IndefIntPress[γf_,γc_,t_]=FullSimplify[Integrate[rho[γf,γc,t],t,Assumptions->{{γf,γc,t} ∈ℝ &&0<γf<1&& γc>γf/(1-Sqrt[γf*(2-γf)])},PrincipalValue->True]]
IntPress1[γf_,γc_]=Evaluate[FullSimplify[IndefIntPress[γf,γc,tmax1[γf,γc]]-IndefIntPress[γf,γc,0], Assumptions->{{γf,γc} ∈ℝ &&0<γf<1&& γc>γf/(1-Sqrt[γf*(2-γf)])}]]

(1/(γf))(t+((γc (((4-((2 γf (-2 ((γc)^2) (-2+γf)+γc (4+(((-2+γf))^2) γf t-2 (Sqrt[4-2 γc (((-2+γf))^2) t]))+2 γf (-2+(Sqrt[4-2 γc (((-2+γf))^2) t]))))/((γc)^2))))^(3/2)))/(6 (((-2+γf))^2) ((γf)^2)))-(1/(γc (((-2+γf))^2) ((γf)^2)))(γc-γf) (((Sqrt[2]) γc+γf (-(Sqrt[2])+(Sqrt[2-γc (((-2+γf))^2) t]))) (Sqrt[2-((γf (-2 ((γc)^2) (-2+γf)+γc (4+(((-2+γf))^2) γf t-2 (Sqrt[4-2 γc (((-2+γf))^2) t]))+2 γf (-2+(Sqrt[4-2 γc (((-2+γf))^2) t]))))/((γc)^2))])+2 γc (-2+γf) γf ArcTanh[((2 γc+γf (-2+(Sqrt[4-2 γc (((-2+γf))^2) t])))/(γc (Sqrt[4-((2 γf (-2 ((γc)^2) (-2+γf)+γc (4+(((-2+γf))^2) γf t-2 (Sqrt[4-2 γc (((-2+γf))^2) t]))+2 γf (-2+(Sqrt[4-2 γc (((-2+γf))^2) t]))))/((γc)^2))])))]))

(1/(3 ((γc)^2) (((-2+γf))^2) ((γf)^3)))(-2 ((γf)^2) (Sqrt[((γc)^2) (((-1+γf))^2)-2 γc γf+((γf)^2)])+2 ((γc)^3) (1+γf (3+2 (-3+γf) γf))+2 γc γf (3 γf+2 (Sqrt[((γc)^2) (((-1+γf))^2)-2 γc γf+((γf)^2)]))+((γc)^2) (-2 (Sqrt[((γc)^2) (((-1+γf))^2)-2 γc γf+((γf)^2)])+2 γf (-3-4 (Sqrt[((γc)^2) (((-1+γf))^2)-2 γc γf+((γf)^2)])+γf (3+2 (Sqrt[((γc)^2) (((-1+γf))^2)-2 γc γf+((γf)^2)]))))+6 ((γc)^2) (γc-γf) (-2+γf) γf (ArcCoth[1-γf]+ArcTanh[((-γc+γf)/(Sqrt[((γc)^2) (((-1+γf))^2)-2 γc γf+((γf)^2)]))]))

Simplify[TrigToExp[(1/(3 ((γc)^2) (((-2+γf))^2) ((γf)^3)))(-2 ((γf)^2) (Sqrt[((γc)^2) (((-1+γf))^2)-2 γc γf+((γf)^2)])+2 ((γc)^3) (1+γf (3+2 (-3+γf) γf))+2 γc γf (3 γf+2 (Sqrt[((γc)^2) (((-1+γf))^2)-2 γc γf+((γf)^2)]))+((γc)^2) (-2 (Sqrt[((γc)^2) (((-1+γf))^2)-2 γc γf+((γf)^2)])+2 γf (-3-4 (Sqrt[((γc)^2) (((-1+γf))^2)-2 γc γf+((γf)^2)])+γf (3+2 (Sqrt[((γc)^2) (((-1+γf))^2)-2 γc γf+((γf)^2)]))))+6 ((γc)^2) (γc-γf) (-2+γf) γf (ArcCoth[1-γf]+ArcTanh[((-γc+γf)/(Sqrt[((γc)^2) (((-1+γf))^2)-2 γc γf+((γf)^2)]))]))]]

(1/(3 ((γc)^2) (((-2+γf))^2) ((γf)^3)))(-2 ((γf)^2) (Sqrt[((γc)^2) (((-1+γf))^2)-2 γc γf+((γf)^2)])+2 ((γc)^3) (1+γf (3+2 (-3+γf) γf))+2 γc γf (3 γf+2 (Sqrt[((γc)^2) (((-1+γf))^2)-2 γc γf+((γf)^2)]))+((γc)^2) (-2 (Sqrt[((γc)^2) (((-1+γf))^2)-2 γc γf+((γf)^2)])+2 γf (-3-4 (Sqrt[((γc)^2) (((-1+γf))^2)-2 γc γf+((γf)^2)])+γf (3+2 (Sqrt[((γc)^2) (((-1+γf))^2)-2 γc γf+((γf)^2)]))))+3 ((γc)^2) (γc-γf) (-2+γf) γf (Log[1+(1/(1-γf))]-Log[((γf)/(-1+γf))]-Log[1+((γc-γf)/(Sqrt[((γc)^2) (((-1+γf))^2)-2 γc γf+((γf)^2)]))]+Log[1+((-γc+γf)/(Sqrt[((γc)^2) (((-1+γf))^2)-2 γc γf+((γf)^2)]))]))

FullSimplify[(1+(1/(1-γf)))*(1+((-γc+γf)/(Sqrt[((γc)^2) (((-1+γf))^2)-2 γc γf+((γf)^2)])))/((((γf)/(-1+γf)))*(1+((γc-γf)/(Sqrt[((γc)^2) (((-1+γf))^2)-2 γc γf+((γf)^2)]))))]

(((-2+γf) (-γc+γf+(Sqrt[((γc)^2) (((-1+γf))^2)-2 γc γf+((γf)^2)])))/(γf (γc-γf+(Sqrt[((γc)^2) (((-1+γf))^2)-2 γc γf+((γf)^2)]))))

IndefIntPress2[γf_,γc_,t_]=FullSimplify[Integrate[rho[γf,γc,t],t,Assumptions->{{γf,γc,t} ∈ℝ &&0<γf<1&& γc<γf/(1-Sqrt[γf*(2-γf)])},PrincipalValue->True]]
IntPress2[γf_,γc_]=Evaluate[IndefIntPress2[γf,γc,tmax2[γf,γc]]-IndefIntPress2[γf,γc,0]]

(1/(γf))(t+((γc (((4-((2 γf (-2 ((γc)^2) (-2+γf)+γc (4+(((-2+γf))^2) γf t-2 (Sqrt[4-2 γc (((-2+γf))^2) t]))+2 γf (-2+(Sqrt[4-2 γc (((-2+γf))^2) t]))))/((γc)^2))))^(3/2)))/(6 (((-2+γf))^2) ((γf)^2)))-(1/(γc (((-2+γf))^2) ((γf)^2)))(γc-γf) (((Sqrt[2]) γc+γf (-(Sqrt[2])+(Sqrt[2-γc (((-2+γf))^2) t]))) (Sqrt[2-((γf (-2 ((γc)^2) (-2+γf)+γc (4+(((-2+γf))^2) γf t-2 (Sqrt[4-2 γc (((-2+γf))^2) t]))+2 γf (-2+(Sqrt[4-2 γc (((-2+γf))^2) t]))))/((γc)^2))])+2 γc (-2+γf) γf ArcTanh[((2 γc+γf (-2+(Sqrt[4-2 γc (((-2+γf))^2) t])))/(γc (Sqrt[4-((2 γf (-2 ((γc)^2) (-2+γf)+γc (4+(((-2+γf))^2) γf t-2 (Sqrt[4-2 γc (((-2+γf))^2) t]))+2 γf (-2+(Sqrt[4-2 γc (((-2+γf))^2) t]))))/((γc)^2))])))]))

-((((γc (((4+4 (-2+γf) γf))^(3/2)))/(6 (((-2+γf))^2) ((γf)^2)))-(((γc-γf) ((Sqrt[2]) γc (Sqrt[2+2 (-2+γf) γf])+2 γc (-2+γf) γf ArcTanh[(2/(Sqrt[4+4 (-2+γf) γf]))]))/(γc (((-2+γf))^2) ((γf)^2))))/(γf))+(1/(γf))(((2 (γc (((-1+γf))^2)+2 (-γc+γf)-2 (Sqrt[(2-γf) γf]) (-γc+γf)))/((((2-γf))^2) ((γf)^2)))+((γc (((4-((2 γf (-2 ((γc)^2) (-2+γf)+γc (4+((2 (((-2+γf))^2) (γc (((-1+γf))^2)+2 (-γc+γf)-2 (Sqrt[(2-γf) γf]) (-γc+γf)))/((((2-γf))^2) γf))-2 (Sqrt[4-((4 γc (((-2+γf))^2) (γc (((-1+γf))^2)+2 (-γc+γf)-2 (Sqrt[(2-γf) γf]) (-γc+γf)))/((((2-γf))^2) ((γf)^2)))]))+2 γf (-2+(Sqrt[4-((4 γc (((-2+γf))^2) (γc (((-1+γf))^2)+2 (-γc+γf)-2 (Sqrt[(2-γf) γf]) (-γc+γf)))/((((2-γf))^2) ((γf)^2)))]))))/((γc)^2))))^(3/2)))/(6 (((-2+γf))^2) ((γf)^2)))-(1/(γc (((-2+γf))^2) ((γf)^2)))(γc-γf) (((Sqrt[2]) γc+γf (-(Sqrt[2])+(Sqrt[2-((2 γc (((-2+γf))^2) (γc (((-1+γf))^2)+2 (-γc+γf)-2 (Sqrt[(2-γf) γf]) (-γc+γf)))/((((2-γf))^2) ((γf)^2)))]))) √(2-(1/((γc)^2))γf (-2 ((γc)^2) (-2+γf)+γc (4+((2 (((-2+γf))^2) (γc (((-1+γf))^2)+2 (-γc+γf)-2 (Sqrt[(2-γf) γf]) (-γc+γf)))/((((2-γf))^2) γf))-2 (Sqrt[4-((4 γc (((-2+γf))^2) (γc (((-1+γf))^2)+2 (-γc+γf)-2 (Sqrt[(2-γf) γf]) (-γc+γf)))/((((2-γf))^2) ((γf)^2)))]))+2 γf (-2+(Sqrt[4-((4 γc (((-2+γf))^2) (γc (((-1+γf))^2)+2 (-γc+γf)-2 (Sqrt[(2-γf) γf]) (-γc+γf)))/((((2-γf))^2) ((γf)^2)))]))))+2 γc (-2+γf) γf ArcTanh[((2 γc+γf (-2+(Sqrt[4-((4 γc (((-2+γf))^2) (γc (((-1+γf))^2)+2 (-γc+γf)-2 (Sqrt[(2-γf) γf]) (-γc+γf)))/((((2-γf))^2) ((γf)^2)))])))/(γc (Sqrt[4-((2 γf (-2 ((γc)^2) (-2+γf)+γc (4+((2 (((-2+γf))^2) (γc (((-1+γf))^2)+2 (-γc+γf)-2 (Sqrt[(2-γf) γf]) (-γc+γf)))/((((2-γf))^2) γf))-2 (Sqrt[4-((4 γc (((-2+γf))^2) (γc (((-1+γf))^2)+2 (-γc+γf)-2 (Sqrt[(2-γf) γf]) (-γc+γf)))/((((2-γf))^2) ((γf)^2)))]))+2 γf (-2+(Sqrt[4-((4 γc (((-2+γf))^2) (γc (((-1+γf))^2)+2 (-γc+γf)-2 (Sqrt[(2-γf) γf]) (-γc+γf)))/((((2-γf))^2) ((γf)^2)))]))))/((γc)^2))])))]))

IntPress2Simpler[γf_,γc_]=FullSimplify[Evaluate[IndefIntPress2[γf,γc,tmax2[γf,γc]]-IndefIntPress2[γf,γc,0]],Assumptions->{0<γf<1,γc>0}]

Power1[γf_,γc_]:=IntPress1[γf,γc]*q[γf]
Throughput1[γf_,γc_]:=tmax1[γf,γc]*q[γf]

Power2[γf_,γc_]:=IntPress2[γf,γc]*q[γf]
Throughput2[γf_,γc_]:=tmax2[γf,γc]*q[γf]

γcc=0.5;
SetOptions[ParametricPlot,BaseStyle->{FontFamily->"Times New Roman",FontSize->20}, FrameTicks->{{0,0.2,0.4,0.6,0.8,1},{0.0,0.5,1.0,1.5,2.0}}, FrameTicksStyle->Directive[Black,20], Ticks->None, Frame->{True,True,False,False}, PlotStyle->Directive[Thickness[0.0015]]]
graph1=ParametricPlot[{Throughput1[γf,γcc],Power1[γf,γcc]},{γf,0.01,γcc*(1+γcc-Sqrt[2*γcc])/(1+γcc^2)},PlotStyle->Blue]
graph2 = ParametricPlot[{Throughput2[γf,γcc],Power2[γf,γcc]},{γf,γcc*(1+γcc-Sqrt[2*γcc])/(1+γcc^2),1},Exclusions->None,PlotPoints->10000,PlotStyle->Blue]

γcc=1.0 ;
SetOptions[ParametricPlot,BaseStyle->{FontFamily->"Times New Roman",FontSize->20}, FrameTicks->{{0,0.2,0.4,0.6,0.8,1},{0.0,0.5,1.0,1.5,2.0}}, FrameTicksStyle->Directive[Black,20], Ticks->None, Frame->{True,True,False,False}, PlotStyle->Directive[Thickness[0.0015]]]
graph1Medium=ParametricPlot[{Throughput1[γf,γcc],Power1[γf,γcc]},{γf,0.01,γcc*(1+γcc-Sqrt[2*γcc])/(1+γcc^2)},PlotStyle->Red]
graph2Medium = ParametricPlot[{Throughput2[γf,γcc],Power2[γf,γcc]},{γf,γcc*(1+γcc-Sqrt[2*γcc])/(1+γcc^2),1},Exclusions->None,PlotPoints->10000, PlotStyle->Red]

γcc=2.0 ;
SetOptions[ParametricPlot,BaseStyle->{FontFamily->"Times New Roman",FontSize->20}, FrameTicks->{{0,0.2,0.4,0.6,0.8,1},{0.0,0.5,1.0,1.5,2.0}}, FrameTicksStyle->Directive[Black,20], Ticks->None, Frame->{True,True,False,False}, PlotStyle->Directive[Thickness[0.0015]]]
graph1High=ParametricPlot[{Throughput1[γf,γcc],Power1[γf,γcc]},{γf,0.01,γcc*(1+γcc-Sqrt[2*γcc])/(1+γcc^2)},PlotStyle->Green]
graph2High = ParametricPlot[{Throughput2[γf,γcc],Power2[γf,γcc]},{γf,γcc*(1+γcc-Sqrt[2*γcc])/(1+γcc^2),1},Exclusions->None,PlotPoints->10000, PlotStyle->Green]

Show[graph2,graph1,graph1Medium,graph2Medium, graph1High, graph2High,PlotRange->All]

SetDirectory["C:\\Users\\admin\\Desktop"]
Export["ThroughputVsPower.eps",Show[graph2,graph1,graph1Medium,graph2Medium, graph1High, graph2High,PlotRange->All]]