cF = 50000;
cL=150;
cC=5;
cD=5000;
theta1=0.4;
nDmax1=50;
nDstar[nDmax_,nA_,theta_]:=nDmax(1-((nDmax-1)/nDmax)^(theta nA));
TotalCost[nA_]:=cF+nA(cL+cC) + cD nDstar[nDmax1,nA,theta1];
AverageCost[nA_,nC_]:=TotalCost[nA/nC]/(nA/nC);

GraphCosts=Plot[{AverageCost[nA,1],AverageCost[nA,2],AverageCost[nA,5]},{nA,1,5000},AxesOrigin->{0,0},GridLines->{{},{155}},PlotRange->{{0,5000},{0,5000}},PlotStyle->{{Blue},{Green,Dashed},{Red,Dotted}},PlotLegends->{"(Subscript[n, C]) = 1","(Subscript[n, C]) = 2","(Subscript[n, C])=5"},AxesLabel->{"(Subscript[n, A])",None}]

Export["U:\\Papers\\Reproducibility\\Cost Analysis\\Graph_Costs.pdf",GraphCosts,"PDF"]