The Strategic Determination of the Supply of Liquid Assets (WP 2016)

Athanasios Geromichalos and Lucas Herrenbrueck (UC Davis working paper, 2016)
This version: June 2017

Using probabilistic-TIOLI bargaining. Note that the code uses λ where the paper uses θ (the probability that the C-type makes the TIOLI offer).

Finding equilibria

The N entry choice - interactive graphics to figure things out

Asset supplies are exogenous. The entry choice of C, as well as buyer portfolios, are optimal conditional on a fixed number eN of N-types entering the A-market. Then, the incentive of an individual N-type to deviate is analyzed. Positive numbers mean that the N-individual prefers to enter the A-market (i.e., the initial guess of eN was “too low”), negative numbers mean B-market.

Manipulate[Module[{uinv,entryC},
uinv=z/.Solve[u[z]==#,z][[1]]&;
entryC[eN_]:=(* ex-ante choice *)Quiet@FindRoot[If[eN<1,If[eN>0,{
i==ℓ(u'[q0A]-1)+λ ℓ δA (((1+(ℓ/(1-ℓ))(eC/eN)))^(ρ-1)) ((((1-ℓ)eN))^ρ)(u'[q1A]-u'[q0A]),
i==ℓ(u'[q0B]-1)+λ ℓ δB (((1+(ℓ/(1-ℓ))((1-eC)/(1-eN))))^(ρ-1)) ((((1-ℓ)(1-eN)))^ρ)(u'[q1B]-u'[q0B]),
q1A==Min[1,(1+A/M)q0A+((A(1-eC))/(M eC))q0B],
q1B==Min[1,(1+B/M)q0B+((B eC)/(M(1-eC)))q0A],
-i q0A-λ ℓ δA (((1+(ℓ/(1-ℓ))(eC/eN)))^(ρ-1)) ((((1-ℓ)eN))^ρ)(u'[q1A]-1)(q1A-q0A)+ℓ(1-λ δA (((1+(ℓ/(1-ℓ))(eC/eN)))^(ρ-1)) ((((1-ℓ)eN))^ρ))(u[q0A]-q0A)+ℓ λ δA (((1+(ℓ/(1-ℓ))(eC/eN)))^(ρ-1)) ((((1-ℓ)eN))^ρ)(u[q1A]-q1A)-(-i q0B-λ ℓ δB (((1+(ℓ/(1-ℓ))((1-eC)/(1-eN))))^(ρ-1)) ((((1-ℓ)(1-eN)))^ρ)(u'[q1B]-1)(q1B-q0B)+ℓ(1-λ δB (((1+(ℓ/(1-ℓ))((1-eC)/(1-eN))))^(ρ-1)) ((((1-ℓ)(1-eN)))^ρ))(u[q0B]-q0B)+ℓ λ δB (((1+(ℓ/(1-ℓ))((1-eC)/(1-eN))))^(ρ-1)) ((((1-ℓ)(1-eN)))^ρ)(u[q1B]-q1B))(*)eC(1-eC)*)==0
}(* middle *),
{i==ℓ(u'[q0A]-1),i==ℓ(u'[q0B]-1)+λ ℓ(1-ℓ)δB(u'[q1B]-u'[q0B]),q1A==q0A,q1B==Min[1,(1+B/M)q0B],eC==0}(* B corner *)],
{i==ℓ(u'[q0A]-1)+λ ℓ(1-ℓ)δA (u'[q1A]-u'[q0A]),i==ℓ(u'[q0B]-1),q1A==Min[1,(1+A/M)q0A],
q1B==q0B,eC==1}(* A corner *)],
{{q0A,.8,0,1},{q1A,.9},{q0B,.8,0,1},{q1B,.9},{eC,(*((δA A)/(δA A+δB B))((δA)/(δA+δB))*)((δA)/(δA+δB))}}];
Quiet@Plot[If[eC<1,If[eC>0,
(*N's surplus from entering A*)((q0A+Min[(A/(eC M))(eC q0A+(1-eC)q0B),u[1]-u[q0A]]-Min[1,uinv[u[q0A]+(A/(eC M))(eC q0A+(1-eC)q0B)]])δA (((1+((1-ℓ)/ℓ)(eN/eC)))^(ρ-1))(((ℓ eC))^ρ)-(*N's surplus from entering B*)(q0B+Min[(B/((1-eC) M))(eC q0A+(1-eC)q0B),u[1]-u[q0B]]-Min[1,uinv[u[q0B]+(B/((1-eC) M))(eC q0A+(1-eC)q0B)]])δB (((1+((1-ℓ)/ℓ)((1-eN)/(1-eC))))^(ρ-1))(((ℓ(1-eC)))^ρ))
/(*divided by sum*)((q0A+Min[(A/(eC M))(eC q0A+(1-eC)q0B),u[1]-u[q0A]]-Min[1,uinv[u[q0A]+(A/(eC M))(eC q0A+(1-eC)q0B)]])δA (((1+((1-ℓ)/ℓ)(eN/eC)))^(ρ-1))(((ℓ eC))^ρ)+(q0B+Min[(B/((1-eC) M))(eC q0A+(1-eC)q0B),u[1]-u[q0B]]-Min[1,uinv[u[q0B]+(B/((1-eC) M))(eC q0A+(1-eC)q0B)]])δB (((1+((1-ℓ)/ℓ)((1-eN)/(1-eC))))^(ρ-1))(((ℓ(1-eC)))^ρ))
,-1],1]/.entryC[eN],{eN,.001,.999},ImageSize->500,AspectRatio->.8,MaxRecursion->2,PlotRange->{{0,1},{-1,1}},AxesLabel->{"eN",None},PlotLabel->"equation (6): difference divided by sum",BaseStyle->{Bold,Medium},PlotStyle->{{Thick,Blue}}]
],{{u,2Sqrt[#]&}},Delimiter,{{ℓ,.5},0,1},{{λ,.5},0,1},{{δA,1},0,1},{{δB,.5},0,1},{{ρ,0},0,1},Delimiter,{{M,1}},{{A,.1},0,.2},{{B,.1},0,.2},{{i,.1},0,.2}]

For deeper understanding: some implied trading frequencies and the C-entry choice, for a fixed N-entry choice.

Manipulate[Module[{uinv,entryC},
uinv=z/.Solve[u[z]==#,z][[1]]&;
entryC[eN_]:=(* ex-ante choice *)Quiet@FindRoot[If[eN<1,If[eN>0,{
i==ℓ(u'[q0A]-1)+λ ℓ δA (((1+(ℓ/(1-ℓ))(eC/eN)))^(ρ-1)) ((((1-ℓ)eN))^ρ)(u'[q1A]-u'[q0A]),
i==ℓ(u'[q0B]-1)+λ ℓ δB (((1+(ℓ/(1-ℓ))((1-eC)/(1-eN))))^(ρ-1)) ((((1-ℓ)(1-eN)))^ρ)(u'[q1B]-u'[q0B]),
q1A==Min[1,(1+A/M)q0A+((A(1-eC))/(M eC))q0B],
q1B==Min[1,(1+B/M)q0B+((B eC)/(M(1-eC)))q0A],
-i q0A-λ ℓ δA (((1+(ℓ/(1-ℓ))(eC/eN)))^(ρ-1)) ((((1-ℓ)eN))^ρ)(u'[q1A]-1)(q1A-q0A)+ℓ(1-λ δA (((1+(ℓ/(1-ℓ))(eC/eN)))^(ρ-1)) ((((1-ℓ)eN))^ρ))(u[q0A]-q0A)+ℓ λ δA (((1+(ℓ/(1-ℓ))(eC/eN)))^(ρ-1)) ((((1-ℓ)eN))^ρ)(u[q1A]-q1A)-(-i q0B-λ ℓ δB (((1+(ℓ/(1-ℓ))((1-eC)/(1-eN))))^(ρ-1)) ((((1-ℓ)(1-eN)))^ρ)(u'[q1B]-1)(q1B-q0B)+ℓ(1-λ δB (((1+(ℓ/(1-ℓ))((1-eC)/(1-eN))))^(ρ-1)) ((((1-ℓ)(1-eN)))^ρ))(u[q0B]-q0B)+ℓ λ δB (((1+(ℓ/(1-ℓ))((1-eC)/(1-eN))))^(ρ-1)) ((((1-ℓ)(1-eN)))^ρ)(u[q1B]-q1B))(*)eC(1-eC)*)==0
}(* middle *),
{i==ℓ(u'[q0A]-1),i==ℓ(u'[q0B]-1)+λ ℓ(1-ℓ)δB(u'[q1B]-u'[q0B]),q1A==q0A,q1B==Min[1,(1+B/M)q0B],eC==0}(* B corner *)],
{i==ℓ(u'[q0A]-1)+λ ℓ(1-ℓ)δA (u'[q1A]-u'[q0A]),i==ℓ(u'[q0B]-1),q1A==Min[1,(1+A/M)q0A],
q1B==q0B,eC==1}(* A corner *)],
{{q0A,.8,0,1},{q1A,.9},{q0B,.8,0,1},{q1B,.9},{eC,(*((δA A)/(δA A+δB B))((δA)/(δA+δB))*)((δA)/(δA+δB))}}];
Quiet@Plot[{δA (((1+(ℓ/(1-ℓ))(eC/eN)))^(ρ-1)) ((((1-ℓ)eN))^ρ)/.entryC[eN],δB (((1+(ℓ/(1-ℓ))((1-eC)/(1-eN))))^(ρ-1)) ((((1-ℓ)(1-eN)))^ρ)/.entryC[eN],eC/.entryC[eN]}
,{eN,.001,.999},ImageSize->500,AspectRatio->.8,MaxRecursion->2,PlotRange->{{0,1},{0,1}},AxesLabel->{"eN",None},PlotLabel->"blue: αCA, red: αCB, yellow: eC",BaseStyle->{Bold,Medium},PlotStyle->{{Thick,Blue},{Thick,Red},{Thick,Yellow}}]
],{{u,Log(*2Sqrt[#]&*)}},Delimiter,{{ℓ,.5},0,1},{{λ,.5},0,1},{{δA,1},0,1},{{δB,1},0,1},{{ρ,0},0,1},Delimiter,{{M,1}},{{A,.1},0,.2},{{B,.1},0,.2},{{i,.1},0,.2}]

The N entry choice - figures for the paper

Module[{u=Log(*2Sqrt[#]&*),uinv=Exp,upinv,ℓ=.4,λ=.5,δA=1,δB=1,ρ=0,M=1,A=.05,B=.28,i=.1,entryC},
entryC[eN_]:=(* ex-ante choice *)Quiet@FindRoot[If[eN<1,If[eN>0,{
i==ℓ(u'[q0A]-1)+λ ℓ δA (((1+(ℓ/(1-ℓ))(eC/eN)))^(ρ-1)) ((((1-ℓ)eN))^ρ)(u'[q1A]-u'[q0A]),
i==ℓ(u'[q0B]-1)+λ ℓ δB (((1+(ℓ/(1-ℓ))((1-eC)/(1-eN))))^(ρ-1)) ((((1-ℓ)(1-eN)))^ρ)(u'[q1B]-u'[q0B]),
q1A==Min[1,(1+A/M)q0A+((A(1-eC))/(M eC))q0B],
q1B==Min[1,(1+B/M)q0B+((B eC)/(M(1-eC)))q0A],
-i q0A-λ ℓ δA (((1+(ℓ/(1-ℓ))(eC/eN)))^(ρ-1)) ((((1-ℓ)eN))^ρ)(u'[q1A]-1)(q1A-q0A)+ℓ(1-λ δA (((1+(ℓ/(1-ℓ))(eC/eN)))^(ρ-1)) ((((1-ℓ)eN))^ρ))(u[q0A]-q0A)+ℓ λ δA (((1+(ℓ/(1-ℓ))(eC/eN)))^(ρ-1)) ((((1-ℓ)eN))^ρ)(u[q1A]-q1A)-(-i q0B-λ ℓ δB (((1+(ℓ/(1-ℓ))((1-eC)/(1-eN))))^(ρ-1)) ((((1-ℓ)(1-eN)))^ρ)(u'[q1B]-1)(q1B-q0B)+ℓ(1-λ δB (((1+(ℓ/(1-ℓ))((1-eC)/(1-eN))))^(ρ-1)) ((((1-ℓ)(1-eN)))^ρ))(u[q0B]-q0B)+ℓ λ δB (((1+(ℓ/(1-ℓ))((1-eC)/(1-eN))))^(ρ-1)) ((((1-ℓ)(1-eN)))^ρ)(u[q1B]-q1B))(*)eC(1-eC)*)==0
}(* middle *),
{i==ℓ(u'[q0A]-1),i==ℓ(u'[q0B]-1)+λ ℓ(1-ℓ)δB(u'[q1B]-u'[q0B]),q1A==q0A,q1B==Min[1,(1+B/M)q0B],eC==0}(* B corner *)],
{i==ℓ(u'[q0A]-1)+λ ℓ(1-ℓ)δA (u'[q1A]-u'[q0A]),i==ℓ(u'[q0B]-1),q1A==Min[1,(1+A/M)q0A],
q1B==q0B,eC==1}(* A corner *)],
{{q0A,.8,0,1},{q1A,.9},{q0B,.8,0,1},{q1B,.9},{eC,(*((δA A)/(δA A+δB B))((δA)/(δA+δB))*)((δA)/(δA+δB))}}];
Quiet@Plot[If[eC<1,If[eC>0,
(*N's surplus from entering A*)((q0A+Min[(A/(eC M))(eC q0A+(1-eC)q0B),u[1]-u[q0A]]-Min[1,uinv[u[q0A]+(A/(eC M))(eC q0A+(1-eC)q0B)]])δA (((1+((1-ℓ)/ℓ)(eN/eC)))^(ρ-1))(((ℓ eC))^ρ)-(*N's surplus from entering B*)(q0B+Min[(B/((1-eC) M))(eC q0A+(1-eC)q0B),u[1]-u[q0B]]-Min[1,uinv[u[q0B]+(B/((1-eC) M))(eC q0A+(1-eC)q0B)]])δB (((1+((1-ℓ)/ℓ)((1-eN)/(1-eC))))^(ρ-1))(((ℓ(1-eC)))^ρ))
/(*divided by sum*)((q0A+Min[(A/(eC M))(eC q0A+(1-eC)q0B),u[1]-u[q0A]]-Min[1,uinv[u[q0A]+(A/(eC M))(eC q0A+(1-eC)q0B)]])δA (((1+((1-ℓ)/ℓ)(eN/eC)))^(ρ-1))(((ℓ eC))^ρ)+(q0B+Min[(B/((1-eC) M))(eC q0A+(1-eC)q0B),u[1]-u[q0B]]-Min[1,uinv[u[q0B]+(B/((1-eC) M))(eC q0A+(1-eC)q0B)]])δB (((1+((1-ℓ)/ℓ)((1-eN)/(1-eC))))^(ρ-1))(((ℓ(1-eC)))^ρ))
,-1],1]/.entryC[eN],{eN,.001,.999},ImageSize->400,AspectRatio->1,MaxRecursion->3,PlotRange->{{-.02,1.02},{-1.02,1.02}},AxesLabel->{"(Subscript[e, N])",None},(*PlotLabel->"equation (6): difference divided by sum",*)BaseStyle->{Bold,FontSize->24},PlotStyle->{{Thick,DotDashed,Blue}},AxesStyle->Opacity[0],Epilog->{PointSize@Large,Blue,Point[{0,-1}],Point[{1,1}],(*Text["(((Subscript[δ, NA])(Subscript[S, NA]) - (Subscript[δ, NB])(Subscript[S, NB]))/((Subscript[δ, NA])(Subscript[S, NA]) + (Subscript[δ, NB])(Subscript[S, NB])))",{.4,-.4}],Black,Text[Style["ρ = "<>ToString[ρ]],{.5,.95}],*)
Blue,FontSize->20,Text["high B, low A",{.3,-.5}]}]
]

Module[{u=Log(*2Sqrt[#]&*),uinv=Exp,upinv,ℓ=.4,λ=.5,δA=1,δB=1,ρ=0,M=1,A=.28,B=.05,i=.1,entryC},
entryC[eN_]:=(* ex-ante choice *)Quiet@FindRoot[If[eN<1,If[eN>0,{
i==ℓ(u'[q0A]-1)+λ ℓ δA (((1+(ℓ/(1-ℓ))(eC/eN)))^(ρ-1)) ((((1-ℓ)eN))^ρ)(u'[q1A]-u'[q0A]),
i==ℓ(u'[q0B]-1)+λ ℓ δB (((1+(ℓ/(1-ℓ))((1-eC)/(1-eN))))^(ρ-1)) ((((1-ℓ)(1-eN)))^ρ)(u'[q1B]-u'[q0B]),
q1A==Min[1,(1+A/M)q0A+((A(1-eC))/(M eC))q0B],
q1B==Min[1,(1+B/M)q0B+((B eC)/(M(1-eC)))q0A],
-i q0A-λ ℓ δA (((1+(ℓ/(1-ℓ))(eC/eN)))^(ρ-1)) ((((1-ℓ)eN))^ρ)(u'[q1A]-1)(q1A-q0A)+ℓ(1-λ δA (((1+(ℓ/(1-ℓ))(eC/eN)))^(ρ-1)) ((((1-ℓ)eN))^ρ))(u[q0A]-q0A)+ℓ λ δA (((1+(ℓ/(1-ℓ))(eC/eN)))^(ρ-1)) ((((1-ℓ)eN))^ρ)(u[q1A]-q1A)-(-i q0B-λ ℓ δB (((1+(ℓ/(1-ℓ))((1-eC)/(1-eN))))^(ρ-1)) ((((1-ℓ)(1-eN)))^ρ)(u'[q1B]-1)(q1B-q0B)+ℓ(1-λ δB (((1+(ℓ/(1-ℓ))((1-eC)/(1-eN))))^(ρ-1)) ((((1-ℓ)(1-eN)))^ρ))(u[q0B]-q0B)+ℓ λ δB (((1+(ℓ/(1-ℓ))((1-eC)/(1-eN))))^(ρ-1)) ((((1-ℓ)(1-eN)))^ρ)(u[q1B]-q1B))(*)eC(1-eC)*)==0
}(* middle *),
{i==ℓ(u'[q0A]-1),i==ℓ(u'[q0B]-1)+λ ℓ(1-ℓ)δB(u'[q1B]-u'[q0B]),q1A==q0A,q1B==Min[1,(1+B/M)q0B],eC==0}(* B corner *)],
{i==ℓ(u'[q0A]-1)+λ ℓ(1-ℓ)δA (u'[q1A]-u'[q0A]),i==ℓ(u'[q0B]-1),q1A==Min[1,(1+A/M)q0A],
q1B==q0B,eC==1}(* A corner *)],
{{q0A,.8,0,1},{q1A,.9},{q0B,.8,0,1},{q1B,.9},{eC,(*((δA A)/(δA A+δB B))((δA)/(δA+δB))*)((δA)/(δA+δB))}}];
Quiet@Plot[If[eC<1,If[eC>0,
(*N's surplus from entering A*)((q0A+Min[(A/(eC M))(eC q0A+(1-eC)q0B),u[1]-u[q0A]]-Min[1,uinv[u[q0A]+(A/(eC M))(eC q0A+(1-eC)q0B)]])δA (((1+((1-ℓ)/ℓ)(eN/eC)))^(ρ-1))(((ℓ eC))^ρ)-(*N's surplus from entering B*)(q0B+Min[(B/((1-eC) M))(eC q0A+(1-eC)q0B),u[1]-u[q0B]]-Min[1,uinv[u[q0B]+(B/((1-eC) M))(eC q0A+(1-eC)q0B)]])δB (((1+((1-ℓ)/ℓ)((1-eN)/(1-eC))))^(ρ-1))(((ℓ(1-eC)))^ρ))
/(*divided by sum*)((q0A+Min[(A/(eC M))(eC q0A+(1-eC)q0B),u[1]-u[q0A]]-Min[1,uinv[u[q0A]+(A/(eC M))(eC q0A+(1-eC)q0B)]])δA (((1+((1-ℓ)/ℓ)(eN/eC)))^(ρ-1))(((ℓ eC))^ρ)+(q0B+Min[(B/((1-eC) M))(eC q0A+(1-eC)q0B),u[1]-u[q0B]]-Min[1,uinv[u[q0B]+(B/((1-eC) M))(eC q0A+(1-eC)q0B)]])δB (((1+((1-ℓ)/ℓ)((1-eN)/(1-eC))))^(ρ-1))(((ℓ(1-eC)))^ρ))
,-1],1]/.entryC[eN],{eN,.001,.999},ImageSize->400,AspectRatio->1,MaxRecursion->3,PlotRange->{{-.02,1.02},{-1.02,1.02}},AxesLabel->{"(Subscript[e, N])",None},(*PlotLabel->"equation (6): difference divided by sum",*)BaseStyle->{Bold,FontSize->24},PlotStyle->{{Thick,Dashed,Blue}},AxesStyle->Opacity[0],Epilog->{PointSize@Large,Blue,Point[{0,-1}],Point[{1,1}],(*Text["(((Subscript[δ, NA])(Subscript[S, NA]) - (Subscript[δ, NB])(Subscript[S, NB]))/((Subscript[δ, NA])(Subscript[S, NA]) + (Subscript[δ, NB])(Subscript[S, NB])))",{.4,-.4}],Black,Text[Style["ρ = "<>ToString[ρ]],{.5,.95}],*)
Blue,FontSize->20,Text["high A, low B",{.7,.5}](*Text["higher B",{.45,-.4}]*)}]
]

Module[{u=Log(*2Sqrt[#]&*),uinv=Exp,upinv,ℓ=.4,λ=.5,δA=1,δB=1,ρ=0,M=1,A=.05,B=.05,i=.1,entryC},
entryC[eN_]:=(* ex-ante choice *)Quiet@FindRoot[If[eN<1,If[eN>0,{
i==ℓ(u'[q0A]-1)+λ ℓ δA (((1+(ℓ/(1-ℓ))(eC/eN)))^(ρ-1)) ((((1-ℓ)eN))^ρ)(u'[q1A]-u'[q0A]),
i==ℓ(u'[q0B]-1)+λ ℓ δB (((1+(ℓ/(1-ℓ))((1-eC)/(1-eN))))^(ρ-1)) ((((1-ℓ)(1-eN)))^ρ)(u'[q1B]-u'[q0B]),
q1A==Min[1,(1+A/M)q0A+((A(1-eC))/(M eC))q0B],
q1B==Min[1,(1+B/M)q0B+((B eC)/(M(1-eC)))q0A],
-i q0A-λ ℓ δA (((1+(ℓ/(1-ℓ))(eC/eN)))^(ρ-1)) ((((1-ℓ)eN))^ρ)(u'[q1A]-1)(q1A-q0A)+ℓ(1-λ δA (((1+(ℓ/(1-ℓ))(eC/eN)))^(ρ-1)) ((((1-ℓ)eN))^ρ))(u[q0A]-q0A)+ℓ λ δA (((1+(ℓ/(1-ℓ))(eC/eN)))^(ρ-1)) ((((1-ℓ)eN))^ρ)(u[q1A]-q1A)-(-i q0B-λ ℓ δB (((1+(ℓ/(1-ℓ))((1-eC)/(1-eN))))^(ρ-1)) ((((1-ℓ)(1-eN)))^ρ)(u'[q1B]-1)(q1B-q0B)+ℓ(1-λ δB (((1+(ℓ/(1-ℓ))((1-eC)/(1-eN))))^(ρ-1)) ((((1-ℓ)(1-eN)))^ρ))(u[q0B]-q0B)+ℓ λ δB (((1+(ℓ/(1-ℓ))((1-eC)/(1-eN))))^(ρ-1)) ((((1-ℓ)(1-eN)))^ρ)(u[q1B]-q1B))(*)eC(1-eC)*)==0
}(* middle *),
{i==ℓ(u'[q0A]-1),i==ℓ(u'[q0B]-1)+λ ℓ(1-ℓ)δB(u'[q1B]-u'[q0B]),q1A==q0A,q1B==Min[1,(1+B/M)q0B],eC==0}(* B corner *)],
{i==ℓ(u'[q0A]-1)+λ ℓ(1-ℓ)δA (u'[q1A]-u'[q0A]),i==ℓ(u'[q0B]-1),q1A==Min[1,(1+A/M)q0A],
q1B==q0B,eC==1}(* A corner *)],
{{q0A,.8,0,1},{q1A,.9},{q0B,.8,0,1},{q1B,.9},{eC,(*((δA A)/(δA A+δB B))((δA)/(δA+δB))*)((δA)/(δA+δB))}}];
Quiet@Plot[If[eC<1,If[eC>0,
(*N's surplus from entering A*)((q0A+Min[(A/(eC M))(eC q0A+(1-eC)q0B),u[1]-u[q0A]]-Min[1,uinv[u[q0A]+(A/(eC M))(eC q0A+(1-eC)q0B)]])δA (((1+((1-ℓ)/ℓ)(eN/eC)))^(ρ-1))(((ℓ eC))^ρ)-(*N's surplus from entering B*)(q0B+Min[(B/((1-eC) M))(eC q0A+(1-eC)q0B),u[1]-u[q0B]]-Min[1,uinv[u[q0B]+(B/((1-eC) M))(eC q0A+(1-eC)q0B)]])δB (((1+((1-ℓ)/ℓ)((1-eN)/(1-eC))))^(ρ-1))(((ℓ(1-eC)))^ρ))
/(*divided by sum*)((q0A+Min[(A/(eC M))(eC q0A+(1-eC)q0B),u[1]-u[q0A]]-Min[1,uinv[u[q0A]+(A/(eC M))(eC q0A+(1-eC)q0B)]])δA (((1+((1-ℓ)/ℓ)(eN/eC)))^(ρ-1))(((ℓ eC))^ρ)+(q0B+Min[(B/((1-eC) M))(eC q0A+(1-eC)q0B),u[1]-u[q0B]]-Min[1,uinv[u[q0B]+(B/((1-eC) M))(eC q0A+(1-eC)q0B)]])δB (((1+((1-ℓ)/ℓ)((1-eN)/(1-eC))))^(ρ-1))(((ℓ(1-eC)))^ρ))
,-1],1]/.entryC[eN],{eN,.001,.999},ImageSize->400,AspectRatio->1,MaxRecursion->3,PlotRange->{{-.02,1.02},{-1.02,1.02}},AxesLabel->{"(Subscript[e, N])",None},(*PlotLabel->"equation (6): difference divided by sum",*)BaseStyle->{Bold,FontSize->24},PlotStyle->{{Thick,Thick,Blue}},Epilog->{PointSize@Large,Blue,Point[{0,-1}],Point[{1,1}],(*Text["(((Subscript[δ, NA])(Subscript[S, NA]) - (Subscript[δ, NB])(Subscript[S, NB]))/((Subscript[δ, NA])(Subscript[S, NA]) + (Subscript[δ, NB])(Subscript[S, NB])))",{.4,-.4}],*)Black,Text[Style["ρ = "<>ToString[ρ]],{.5,.95}],
Blue,FontSize->20,Text["low A,B",{.52,.16}](*Text["higher A",{.55,.4}]*)(*Text["higher B",{.45,-.4}]*)}]
]

Module[{u=Log(*2Sqrt[#]&*),uinv=Exp,upinv,ℓ=.4,λ=.5,δA=1,δB=1,ρ=0,M=1,A=.28,B=.28,i=.1,entryC},
entryC[eN_]:=(* ex-ante choice *)Quiet@FindRoot[If[eN<1,If[eN>0,{
i==ℓ(u'[q0A]-1)+λ ℓ δA (((1+(ℓ/(1-ℓ))(eC/eN)))^(ρ-1)) ((((1-ℓ)eN))^ρ)(u'[q1A]-u'[q0A]),
i==ℓ(u'[q0B]-1)+λ ℓ δB (((1+(ℓ/(1-ℓ))((1-eC)/(1-eN))))^(ρ-1)) ((((1-ℓ)(1-eN)))^ρ)(u'[q1B]-u'[q0B]),
q1A==Min[1,(1+A/M)q0A+((A(1-eC))/(M eC))q0B],
q1B==Min[1,(1+B/M)q0B+((B eC)/(M(1-eC)))q0A],
-i q0A-λ ℓ δA (((1+(ℓ/(1-ℓ))(eC/eN)))^(ρ-1)) ((((1-ℓ)eN))^ρ)(u'[q1A]-1)(q1A-q0A)+ℓ(1-λ δA (((1+(ℓ/(1-ℓ))(eC/eN)))^(ρ-1)) ((((1-ℓ)eN))^ρ))(u[q0A]-q0A)+ℓ λ δA (((1+(ℓ/(1-ℓ))(eC/eN)))^(ρ-1)) ((((1-ℓ)eN))^ρ)(u[q1A]-q1A)-(-i q0B-λ ℓ δB (((1+(ℓ/(1-ℓ))((1-eC)/(1-eN))))^(ρ-1)) ((((1-ℓ)(1-eN)))^ρ)(u'[q1B]-1)(q1B-q0B)+ℓ(1-λ δB (((1+(ℓ/(1-ℓ))((1-eC)/(1-eN))))^(ρ-1)) ((((1-ℓ)(1-eN)))^ρ))(u[q0B]-q0B)+ℓ λ δB (((1+(ℓ/(1-ℓ))((1-eC)/(1-eN))))^(ρ-1)) ((((1-ℓ)(1-eN)))^ρ)(u[q1B]-q1B))(*)eC(1-eC)*)==0
}(* middle *),
{i==ℓ(u'[q0A]-1),i==ℓ(u'[q0B]-1)+λ ℓ(1-ℓ)δB(u'[q1B]-u'[q0B]),q1A==q0A,q1B==Min[1,(1+B/M)q0B],eC==0}(* B corner *)],
{i==ℓ(u'[q0A]-1)+λ ℓ(1-ℓ)δA (u'[q1A]-u'[q0A]),i==ℓ(u'[q0B]-1),q1A==Min[1,(1+A/M)q0A],
q1B==q0B,eC==1}(* A corner *)],
{{q0A,.8,0,1},{q1A,.9},{q0B,.8,0,1},{q1B,.9},{eC,(*((δA A)/(δA A+δB B))((δA)/(δA+δB))*)((δA)/(δA+δB))}}];
Quiet@Plot[If[eC<1,If[eC>0,
(*N's surplus from entering A*)((q0A+Min[(A/(eC M))(eC q0A+(1-eC)q0B),u[1]-u[q0A]]-Min[1,uinv[u[q0A]+(A/(eC M))(eC q0A+(1-eC)q0B)]])δA (((1+((1-ℓ)/ℓ)(eN/eC)))^(ρ-1))(((ℓ eC))^ρ)-(*N's surplus from entering B*)(q0B+Min[(B/((1-eC) M))(eC q0A+(1-eC)q0B),u[1]-u[q0B]]-Min[1,uinv[u[q0B]+(B/((1-eC) M))(eC q0A+(1-eC)q0B)]])δB (((1+((1-ℓ)/ℓ)((1-eN)/(1-eC))))^(ρ-1))(((ℓ(1-eC)))^ρ))
/(*divided by sum*)((q0A+Min[(A/(eC M))(eC q0A+(1-eC)q0B),u[1]-u[q0A]]-Min[1,uinv[u[q0A]+(A/(eC M))(eC q0A+(1-eC)q0B)]])δA (((1+((1-ℓ)/ℓ)(eN/eC)))^(ρ-1))(((ℓ eC))^ρ)+(q0B+Min[(B/((1-eC) M))(eC q0A+(1-eC)q0B),u[1]-u[q0B]]-Min[1,uinv[u[q0B]+(B/((1-eC) M))(eC q0A+(1-eC)q0B)]])δB (((1+((1-ℓ)/ℓ)((1-eN)/(1-eC))))^(ρ-1))(((ℓ(1-eC)))^ρ))
,-1],1]/.entryC[eN],{eN,.001,.999},ImageSize->400,AspectRatio->1,MaxRecursion->3,PlotRange->{{-.02,1.02},{-1.02,1.02}},AxesLabel->{"(Subscript[e, N])",None},(*PlotLabel->"equation (6): difference divided by sum",*)BaseStyle->{Bold,FontSize->24},PlotStyle->{{Thick,Dashed,Blue}},AxesStyle->Opacity[0],Epilog->{PointSize@Large,Blue,Point[{0,-1}],Point[{1,1}],(*Text["(((Subscript[δ, NA])(Subscript[S, NA]) - (Subscript[δ, NB])(Subscript[S, NB]))/((Subscript[δ, NA])(Subscript[S, NA]) + (Subscript[δ, NB])(Subscript[S, NB])))",{.4,-.4}],Black,Text[Style["ρ = "<>ToString[ρ]],{.5,.95}],*)
Blue,FontSize->20,Text["high A, B",{.89,-.25}](*Text["higher A",{.55,.4}]*)(*Text["higher B",{.45,-.4}]*)}]
]

Module[{u=Log(*2Sqrt[#]&*),uinv=Exp,upinv,ℓ=.5,λ=.5,δA=1,δB=1,ρ=1,M=1,A=.05,B=.05,i=.1,entryC},
entryC[eN_]:=(* ex-ante choice *)Quiet@FindRoot[If[eN<1,If[eN>0,{
i==ℓ(u'[q0A]-1)+λ ℓ δA (((1+(ℓ/(1-ℓ))(eC/eN)))^(ρ-1)) ((((1-ℓ)eN))^ρ)(u'[q1A]-u'[q0A]),
i==ℓ(u'[q0B]-1)+λ ℓ δB (((1+(ℓ/(1-ℓ))((1-eC)/(1-eN))))^(ρ-1)) ((((1-ℓ)(1-eN)))^ρ)(u'[q1B]-u'[q0B]),
q1A==Min[1,(1+A/M)q0A+((A(1-eC))/(M eC))q0B],
q1B==Min[1,(1+B/M)q0B+((B eC)/(M(1-eC)))q0A],
-i q0A-λ ℓ δA (((1+(ℓ/(1-ℓ))(eC/eN)))^(ρ-1)) ((((1-ℓ)eN))^ρ)(u'[q1A]-1)(q1A-q0A)+ℓ(1-λ δA (((1+(ℓ/(1-ℓ))(eC/eN)))^(ρ-1)) ((((1-ℓ)eN))^ρ))(u[q0A]-q0A)+ℓ λ δA (((1+(ℓ/(1-ℓ))(eC/eN)))^(ρ-1)) ((((1-ℓ)eN))^ρ)(u[q1A]-q1A)-(-i q0B-λ ℓ δB (((1+(ℓ/(1-ℓ))((1-eC)/(1-eN))))^(ρ-1)) ((((1-ℓ)(1-eN)))^ρ)(u'[q1B]-1)(q1B-q0B)+ℓ(1-λ δB (((1+(ℓ/(1-ℓ))((1-eC)/(1-eN))))^(ρ-1)) ((((1-ℓ)(1-eN)))^ρ))(u[q0B]-q0B)+ℓ λ δB (((1+(ℓ/(1-ℓ))((1-eC)/(1-eN))))^(ρ-1)) ((((1-ℓ)(1-eN)))^ρ)(u[q1B]-q1B))(*)eC(1-eC)*)==0
}(* middle *),
{i==ℓ(u'[q0A]-1),i==ℓ(u'[q0B]-1)+λ ℓ(1-ℓ)δB(u'[q1B]-u'[q0B]),q1A==q0A,q1B==Min[1,(1+B/M)q0B],eC==0}(* B corner *)],
{i==ℓ(u'[q0A]-1)+λ ℓ(1-ℓ)δA (u'[q1A]-u'[q0A]),i==ℓ(u'[q0B]-1),q1A==Min[1,(1+A/M)q0A],
q1B==q0B,eC==1}(* A corner *)],
{{q0A,.8,0,1},{q1A,.9},{q0B,.8,0,1},{q1B,.9},{eC,(*((δA A)/(δA A+δB B))((δA)/(δA+δB))*)((δA)/(δA+δB))}}];
Quiet@Plot[If[eC<1,If[eC>0,
(*N's surplus from entering A*)((q0A+Min[(A/(eC M))(eC q0A+(1-eC)q0B),u[1]-u[q0A]]-Min[1,uinv[u[q0A]+(A/(eC M))(eC q0A+(1-eC)q0B)]])δA (((1+((1-ℓ)/ℓ)(eN/eC)))^(ρ-1))(((ℓ eC))^ρ)-(*N's surplus from entering B*)(q0B+Min[(B/((1-eC) M))(eC q0A+(1-eC)q0B),u[1]-u[q0B]]-Min[1,uinv[u[q0B]+(B/((1-eC) M))(eC q0A+(1-eC)q0B)]])δB (((1+((1-ℓ)/ℓ)((1-eN)/(1-eC))))^(ρ-1))(((ℓ(1-eC)))^ρ))
/(*divided by sum*)((q0A+Min[(A/(eC M))(eC q0A+(1-eC)q0B),u[1]-u[q0A]]-Min[1,uinv[u[q0A]+(A/(eC M))(eC q0A+(1-eC)q0B)]])δA (((1+((1-ℓ)/ℓ)(eN/eC)))^(ρ-1))(((ℓ eC))^ρ)+(q0B+Min[(B/((1-eC) M))(eC q0A+(1-eC)q0B),u[1]-u[q0B]]-Min[1,uinv[u[q0B]+(B/((1-eC) M))(eC q0A+(1-eC)q0B)]])δB (((1+((1-ℓ)/ℓ)((1-eN)/(1-eC))))^(ρ-1))(((ℓ(1-eC)))^ρ))
,-1],1]/.entryC[eN],{eN,.001,.999},ImageSize->400,AspectRatio->1,MaxRecursion->2,PlotRange->{{-.02,1.02},{-1.02,1.02}},AxesLabel->{"(Subscript[e, N])",None},(*PlotLabel->"equation (6): difference divided by sum",*)BaseStyle->{Bold,FontSize->24},PlotStyle->{{Thick,Thick,Blue}},Epilog->{PointSize@Large,Blue,Point[{0,-1}],Point[{1,1}],(*Text["(((Subscript[δ, NA])(Subscript[S, NA]) - (Subscript[δ, NB])(Subscript[S, NB]))/((Subscript[δ, NA])(Subscript[S, NA]) + (Subscript[δ, NB])(Subscript[S, NB])))",{.4,-.4}],*)Black,Text[Style["ρ = "<>ToString[ρ]],{.5,.95}]}]
]

General equilibrium conditional on the asset supplies

The inverse of the utility function (uinv) needs to be provided as an input to the GE-function, to speed up computation later.

Clear[ge]
ge[{u_,uinv_,ℓ_,λ_,δA_,δB_,ρ_,M_,A_,B_,i_,prec_}]:=
Module[{entryC,nextN,equil},
entryC[eN_]:=(* ex-ante choice *)Quiet@FindRoot[If[eN<1,If[eN>0,{
i==ℓ(u'[q0A]-1)+λ ℓ (δA (((1+(ℓ/(1-ℓ))(eC/eN)))^(ρ-1)) ((((1-ℓ)eN))^ρ))(u'[q1A]-u'[q0A]),
i==ℓ(u'[q0B]-1)+λ ℓ (δB (((1+(ℓ/(1-ℓ))((1-eC)/(1-eN))))^(ρ-1)) ((((1-ℓ)(1-eN)))^ρ))(u'[q1B]-u'[q0B]),
q1A==Min[1,(1+A/M)q0A+((A(1-eC))/(M eC))q0B],
q1B==Min[1,(1+B/M)q0B+((B eC)/(M(1-eC)))q0A],
-i q0A-λ ℓ (δA (((1+(ℓ/(1-ℓ))(eC/eN)))^(ρ-1)) ((((1-ℓ)eN))^ρ))(u'[q1A]-1)(q1A-q0A)+ℓ(1-λ (δA (((1+(ℓ/(1-ℓ))(eC/eN)))^(ρ-1)) ((((1-ℓ)eN))^ρ)))(u[q0A]-q0A)+ℓ λ (δA (((1+(ℓ/(1-ℓ))(eC/eN)))^(ρ-1)) ((((1-ℓ)eN))^ρ))(u[q1A]-q1A)-(-i q0B-λ ℓ (δB (((1+(ℓ/(1-ℓ))((1-eC)/(1-eN))))^(ρ-1)) ((((1-ℓ)(1-eN)))^ρ))(u'[q1B]-1)(q1B-q0B)+ℓ(1-λ (δB (((1+(ℓ/(1-ℓ))((1-eC)/(1-eN))))^(ρ-1)) ((((1-ℓ)(1-eN)))^ρ)))(u[q0B]-q0B)+ℓ λ (δB (((1+(ℓ/(1-ℓ))((1-eC)/(1-eN))))^(ρ-1)) ((((1-ℓ)(1-eN)))^ρ))(u[q1B]-q1B))(*)eC(1-eC)*)==0
}(* middle *),
{i==ℓ(u'[q0A]-1),i==ℓ(u'[q0B]-1)+λ ℓ(1-ℓ)δB(u'[q1B]-u'[q0B]),(*q1A==q0A*)q1A==1,q1B==Min[1,(1+B/M)q0B],eC==0}(* B corner *)],
{i==ℓ(u'[q0A]-1)+λ ℓ(1-ℓ)δA (u'[q1A]-u'[q0A]),i==ℓ(u'[q0B]-1),q1A==Min[1,(1+A/M)q0A],
(*q1B==q0B*)q1B==1,eC==1}(* A corner *)],
{{q0A,.8,0,1},{q1A,.9},{q0B,.8,0,1},{q1B,.9},{eC,((δA A)/(δA A+δB B))}}];
nextN[eN_]:=With[{ec=entryC[eN]},If[eN<1&&Evaluate[eC/.ec]<1,If[eN>0&&Evaluate[eC/.ec]>0,Max[0,Min[1,eN+
.5((q0A+Min[(A/(eC M))(eC q0A+(1-eC)q0B),u[1]-u[q0A]]-Min[1,uinv[u[q0A]+(A/(eC M))(eC q0A+(1-eC)q0B)]])(δA (((1+((1-ℓ)/ℓ)(eN/eC)))^(ρ-1))(((ℓ eC))^ρ))-(q0B+Min[(B/((1-eC) M))(eC q0A+(1-eC)q0B),u[1]-u[q0B]]-Min[1,uinv[u[q0B]+(B/((1-eC) M))(eC q0A+(1-eC)q0B)]])(δB (((1+((1-ℓ)/ℓ)((1-eN)/(1-eC))))^(ρ-1))(((ℓ(1-eC)))^ρ)))/((q0A+Min[(A/(eC M))(eC q0A+(1-eC)q0B),u[1]-u[q0A]]-Min[1,uinv[u[q0A]+(A/(eC M))(eC q0A+(1-eC)q0B)]])(δA (((1+((1-ℓ)/ℓ)(eN/eC)))^(ρ-1))(((ℓ eC))^ρ))+(q0B+Min[(B/((1-eC) M))(eC q0A+(1-eC)q0B),u[1]-u[q0B]]-Min[1,uinv[u[q0B]+(B/((1-eC) M))(eC q0A+(1-eC)q0B)]])(δB (((1+((1-ℓ)/ℓ)((1-eN)/(1-eC))))^(ρ-1))(((ℓ(1-eC)))^ρ)))/.ec
]](* middle *),0(* B corner *)],1(* A corner *)]];
With[{en=FixedPoint[nextN,((δA A)/(δA A+δB B))(*((δA A)/(δA A+δB B)),((δA)/(δA+δB))*),10prec,SameTest->(Abs[#1-#2]<10^(-prec)&)]},entryC[en]~Join~{eN->en}]]

Check if it works, and how long it takes:

{AbsoluteTiming@ge[{Log,Exp,.1,.5,1,.99,0,1,.1,.1,.1,4}],AbsoluteTiming@ge[{2Sqrt[#]&,#^2/4&,.5,.5,1,.5,0,1,.01,.1,.1,4}]}

Analyzing the GE results (interpret as asset demand) - interactive

Module[{u=Log(*2Sqrt[#]&,Log*),uinv=Exp(*((#)^2)/4&,Exp*),upinv,ℓ=.5,λ=.5,δA=1,δB=.95,ρ=0.,M=1,i=.1,prec=6,gridlength=20,Abar,Bbar,Agrid,Bgrid,gridresults},
(*uinv=z/.Solve[u[z]==#,z][[1]]&;*)
upinv=z/.Solve[u'[z]==#,z][[1]]&;
Abar=M(1/upinv[1+i/(ℓ-λ ℓ(1-ℓ)δA)]-1);Bbar=(*M(1/upinv[1+i/(ℓ-λ ℓ(1-ℓ)δB)]-1)*)Abar;
Agrid=Range[gridlength]Abar/gridlength;Bgrid=Range[gridlength]Bbar/gridlength;
gridresults=Table[ge[{u,uinv,ℓ,λ,δA,δB,ρ,M,A,B,i,prec}],{A,Agrid},{B,Bgrid}];
Manipulate[With[{resultsA=gridresults[[;;,B]],resultsB=gridresults[[A,;;]]},
GraphicsGrid[{{ListPlot[Transpose[{{Agrid,100If[eN>0,λ ℓ δA (((1+(ℓ/(1-ℓ))(eC/eN)))^(ρ-1)) ((((1-ℓ)eN))^ρ)(((u'[q1A]-1)/((1-λ)u'[q1A]+λ))),0]/.resultsA},
{Agrid,100If[eN<1,λ ℓ δB (((1+(ℓ/(1-ℓ))((1-eC)/(1-eN))))^(ρ-1)) ((((1-ℓ)(1-eN)))^ρ)(((u'[q1B]-1)/((1-λ)u'[q1B]+λ))),0]/.resultsA}},{1,3,2}],Joined->True,ImageSize->300,GridLines->{{Bgrid[[B]]},None},PlotRange->{{0,1.05Abar},{0,2.5}},PlotStyle->{{Thick,Blue},{Thick,Dashed,Red}},AxesLabel->{"A","(Subscript[nlp, A]),(Subscript[nlp, B])"}],
ListPlot[Transpose[{{Bgrid,100If[eN>0,λ ℓ δA (((1+(ℓ/(1-ℓ))(eC/eN)))^(ρ-1)) ((((1-ℓ)eN))^ρ)(((u'[q1A]-1)/((1-λ)u'[q1A]+λ))),0]/.resultsB},
{Bgrid,100If[eN<1,λ ℓ δB (((1+(ℓ/(1-ℓ))((1-eC)/(1-eN))))^(ρ-1)) ((((1-ℓ)(1-eN)))^ρ)(((u'[q1B]-1)/((1-λ)u'[q1B]+λ))),0]/.resultsB}},{1,3,2}],Joined->True,ImageSize->300,PlotRange->{{0,1.05Bbar},{0,2.5}},GridLines->{{Agrid[[A]]},None},PlotStyle->{{Thick,Blue},{Thick,Dashed,Red}},AxesLabel->{"B","(Subscript[nlp, A]),(Subscript[nlp, B])"}]},
{ListPlot[Transpose[{{Agrid,eC/.resultsA},{Agrid,eN/.resultsA}},{1,3,2}],
Joined->True,ImageSize->300,GridLines->{{Bgrid[[B]]},None},PlotRange->{{0,1.05Abar},{0,1}},PlotStyle->{{Thick,Blue},{Thick,Dashed,Red}},AxesLabel->{"A","(Subscript[H, C]),(Subscript[H, N])"}],
ListPlot[Transpose[{{Bgrid,eC/.resultsB},{Bgrid,eN/.resultsB}},{1,3,2}],
Joined->True,ImageSize->300,GridLines->{{Agrid[[A]]},None},PlotRange->{{0,1.05Bbar},{0,1}},PlotStyle->{{Thick,Blue},{Thick,Dashed,Red}},AxesLabel->{"B","(Subscript[H, C]),(Subscript[H, N])"}]}}]],
{{A,Quotient[gridlength ,3]},1,gridlength,1},{{B,Quotient[gridlength ,3]},1,gridlength,1},ControlPlacement->Left(*,SaveDefinitions->True,SynchronousUpdating->False*)]]

With ρ=0 and δB=0.9, we clearly see that eC<eN everywhere.

Module[{u=Log(*2Sqrt[#]&*),uinv=Exp,upinv,ℓ=.5,λ=.5,δA=1,δB=1,ρ=0.0,M=1,i=.1,prec=6,gridlength=20,Abar,Bbar,Agrid,Bgrid,gridresults},
upinv=z/.Solve[u'[z]==#,z][[1]]&;
Abar=M(1/upinv[1+i/(ℓ-λ ℓ(1-ℓ)δA)]-1);Bbar=(*M(1/upinv[1+i/(ℓ-λ ℓ(1-ℓ)δB)]-1)*)Abar;
Agrid=Range[gridlength]Abar/gridlength;Bgrid=Range[gridlength]Bbar/gridlength;
ge[{u,uinv,ℓ,λ,δA,δB,ρ,M,.1,.05,i,prec}]]

Asset demand functions - figures for the paper

Module[{u=Log(*2Sqrt[#]&*),uinv=Exp,upinv,ℓ=.5,λ=.5,δA=1,δB=1,ρ=0.0,M=1,i=.1,prec=6,gridlength=20,Abar,Bbar,Agrid,Bgrid,gridresults},
upinv=z/.Solve[u'[z]==#,z][[1]]&;
Abar=M(1/upinv[1+i/(ℓ-λ ℓ(1-ℓ)δA)]-1);Bbar=(*M(1/upinv[1+i/(ℓ-λ ℓ(1-ℓ)δB)]-1)*)Abar;
Agrid=Range[gridlength]Abar/gridlength;Bgrid=Range[gridlength]Bbar/gridlength;
gridresults=Table[ge[{u,uinv,ℓ,λ,δA,δB,ρ,M,A,B,i,prec}],{A,Agrid},{B,Bgrid}];
With[{B=Quotient[gridlength ,3],resultsA=gridresults[[;;,Quotient[gridlength ,3]]]},
GraphicsGrid[{{ListPlot[Transpose[{{Agrid,100If[eN>0,λ ℓ δA (((1+(ℓ/(1-ℓ))(eC/eN)))^(ρ-1)) ((((1-ℓ)eN))^ρ)(((u'[q1A]-1)/((1-λ)u'[q1A]+λ))),0]/.resultsA},
{Agrid,100If[eN<1,λ ℓ δB (((1+(ℓ/(1-ℓ))((1-eC)/(1-eN))))^(ρ-1)) ((((1-ℓ)(1-eN)))^ρ)(((u'[q1B]-1)/((1-λ)u'[q1B]+λ))),0]/.resultsA}},{1,3,2}],Joined->True,ImageSize->250,GridLines->{{Bgrid[[B]]},None},Ticks->{{0,{Bgrid[[B]],Style["B",Darker@Gray]},.1,.2},{0,1,2}},PlotRange->{{0,1.01Abar},{0,2.01}},PlotStyle->{{Thick,Darker@Blue},{Thick,Dashed,Darker@Red}},BaseStyle->{Bold,FontSize->16},AxesLabel->{"A",None},Epilog->{Darker@Blue,Text["(Subscript[L, A])",{.2,1}],Darker@Red,Text["(Subscript[L, B])",{.14,1}],Black,Text["=",{.17,1}]}]},
{ListPlot[Transpose[{{Agrid,eC/.resultsA},{Agrid,eN/.resultsA}},{1,3,2}],
Joined->True,ImageSize->250,GridLines->{{Bgrid[[B]]},None},Ticks->{{0,{Bgrid[[B]],Style["B",Darker@Gray]},.1,.2},{0,1}},PlotRange->{{0,1.01Abar},{0,1}},PlotStyle->{{Thick,Darker@Brown},{Thick,Dashed,Darker@Green}},BaseStyle->{Bold,FontSize->16},AxesLabel->{"A",None},Epilog->{Darker@Brown,Text["(Subscript[e, C])",{.1,.45}],Darker@Green,Text["(Subscript[e, N])",{.16,.45}],Black,Text["=",{.13,.45}]}]}}]]]
Export["demandr0a1.pdf",%]

Module[{u=Log(*2Sqrt[#]&*),uinv=Exp,upinv,ℓ=.5,λ=.5,δA=1,δB=1,ρ=0.25,M=1,i=.1,prec=6,gridlength=20,Abar,Bbar,Agrid,Bgrid,gridresults},
upinv=z/.Solve[u'[z]==#,z][[1]]&;
Abar=M(1/upinv[1+i/(ℓ-λ ℓ(1-ℓ)δA)]-1);Bbar=(*M(1/upinv[1+i/(ℓ-λ ℓ(1-ℓ)δB)]-1)*)Abar;
Agrid=Range[gridlength]Abar/gridlength;Bgrid=Range[gridlength]Bbar/gridlength;
gridresults=Table[ge[{u,uinv,ℓ,λ,δA,δB,ρ,M,A,B,i,prec}],{A,Agrid},{B,Bgrid}];
With[{B=Quotient[gridlength ,3],resultsA=gridresults[[;;,Quotient[gridlength ,3]]]},
GraphicsGrid[{{ListPlot[Transpose[{{Agrid,100If[eN>0,λ ℓ δA (((1+(ℓ/(1-ℓ))(eC/eN)))^(ρ-1)) ((((1-ℓ)eN))^ρ)(((u'[q1A]-1)/((1-λ)u'[q1A]+λ))),0]/.resultsA},
{Agrid,100If[eN<1,λ ℓ δB (((1+(ℓ/(1-ℓ))((1-eC)/(1-eN))))^(ρ-1)) ((((1-ℓ)(1-eN)))^ρ)(((u'[q1B]-1)/((1-λ)u'[q1B]+λ))),0]/.resultsA}},{1,3,2}],Joined->True,ImageSize->250,GridLines->{{Bgrid[[B]]},None},Ticks->{{0,{Bgrid[[B]],Style["B",Darker@Gray]},.1,.2},{0,1,2}},PlotRange->{{0,1.01Abar},{0,2.01}},PlotStyle->{{Thick,Darker@Blue},{Thick,Dashed,Darker@Red}},BaseStyle->{Bold,FontSize->16},AxesLabel->{"A",None},Epilog->{Darker@Blue,Text["(Subscript[L, A])",{.2,1}],Darker@Red,Text["(Subscript[L, B])",{.05,1.6}]}]},
{ListPlot[Transpose[{{Agrid,eC/.resultsA},{Agrid,eN/.resultsA}},{1,3,2}],
Joined->True,ImageSize->250,GridLines->{{Bgrid[[B]]},None},Ticks->{{0,{Bgrid[[B]],Style["B",Darker@Gray]},.1,.2},{0,1}},PlotRange->{{0,1.01Abar},{0,1}},PlotStyle->{{Thick,Darker@Brown},{Thick,Dashed,Darker@Green}},BaseStyle->{Bold,FontSize->16},AxesLabel->{"A",None},Epilog->{Darker@Brown,Text["(Subscript[e, C])",{.1,.42}],Darker@Green,Text["(Subscript[e, N])",{.1,.78}]}]}}]]]
Export["demandr2a1.pdf",%]

Module[{u=Log(*2Sqrt[#]&*),uinv=Exp,upinv,ℓ=.5,λ=.5,δA=1,δB=.9,ρ=0.,M=1,i=.1,prec=6,gridlength=20,Abar,Bbar,Agrid,Bgrid,gridresults},
upinv=z/.Solve[u'[z]==#,z][[1]]&;
Abar=M(1/upinv[1+i/(ℓ-λ ℓ(1-ℓ)δA)]-1);Bbar=(*M(1/upinv[1+i/(ℓ-λ ℓ(1-ℓ)δB)]-1)*)Abar;
Agrid=Range[gridlength]Abar/gridlength;Bgrid=Range[gridlength]Bbar/gridlength;
gridresults=Table[ge[{u,uinv,ℓ,λ,δA,δB,ρ,M,A,B,i,prec}],{A,Agrid},{B,Bgrid}];
With[{B=Quotient[gridlength ,3],resultsA=gridresults[[;;,Quotient[gridlength ,3]]]},
GraphicsGrid[{{ListPlot[Transpose[{{Agrid,100If[eN>0,λ ℓ δA (((1+(ℓ/(1-ℓ))(eC/eN)))^(ρ-1)) ((((1-ℓ)eN))^ρ)(((u'[q1A]-1)/((1-λ)u'[q1A]+λ))),0]/.resultsA},
{Agrid,100If[eN<1,λ ℓ δB (((1+(ℓ/(1-ℓ))((1-eC)/(1-eN))))^(ρ-1)) ((((1-ℓ)(1-eN)))^ρ)(((u'[q1B]-1)/((1-λ)u'[q1B]+λ))),0]/.resultsA}},{1,3,2}],Joined->True,ImageSize->250,GridLines->{{Bgrid[[B]]},None},Ticks->{{0,{Bgrid[[B]],Style["B",Darker@Gray]},.1,.2},{0,1,2}},PlotRange->{{0,1.01Abar},{0,2.01}},PlotStyle->{{Thick,Darker@Blue},{Thick,Dashed,Darker@Red}},BaseStyle->{Bold,FontSize->16},AxesLabel->{"A",None},Epilog->{Darker@Blue,Text["(Subscript[L, A])",{.2,.55}],Darker@Red,Text["(Subscript[L, B])",{.05,.55}]}]},
{ListPlot[Transpose[{{Agrid,eC/.resultsA},{Agrid,eN/.resultsA}},{1,3,2}],
Joined->True,ImageSize->250,GridLines->{{Bgrid[[B]]},None},Ticks->{{0,{Bgrid[[B]],Style["B",Darker@Gray]},.1,.2},{0,1}},PlotRange->{{0,1.01Abar},{0,1}},PlotStyle->{{Thick,Darker@Brown},{Thick,Dashed,Darker@Green}},BaseStyle->{Bold,FontSize->16},AxesLabel->{"A",None},Epilog->{Darker@Brown,Text["(Subscript[e, C])",{.1,.45}],Darker@Green,Text["(Subscript[e, N])",{.1,.75}]}]}}]]]
Export["demandr0a9.pdf",%]

The “playing field” of the game between two issuers - figures for the paper

(* GRID SEARCH METHOD, BY POSITION *)

Balanced CRS benchmark:

Module[{u=Log(*2Sqrt[#]&*),uinv=Exp,upinv,ℓ=.5,λ=.5,δA=1,δB=1,ρ=0.0,M=1,i=.1,prec=13,gridlength=21,Abar,Bbar,Agrid,Bgrid,gridresults,Aseign,Bseign,ALP,findCNEa,findCNEb},
upinv=z/.Solve[u'[z]==#,z][[1]]&;
Abar=.98 M(1/upinv[1+i/(ℓ-λ ℓ(1-ℓ)δA)]-1);Bbar=Abar;(*M(1/upinv[1+i/(ℓ-λ ℓ(1-ℓ)δB)]-1);*)
(* using the .98 to get positive (if tiny) profits at the upper bound. looks better on the graph to have the final column filled in faintly rather than all white. *)
Agrid=Range[gridlength]Abar/gridlength;Bgrid=Range[gridlength]Bbar/gridlength;
gridresults=Table[ge[{u,uinv,ℓ,λ,δA,δB,ρ,M,A,B,i,prec}],{A,Agrid},{B,Bgrid}];
Aseign=(If[eN>0,λ ℓ δA (((1+(ℓ/(1-ℓ))(eC/eN)))^(ρ-1)) ((((1-ℓ)eN))^ρ)(u'[q1A]-1),0]/.gridresults)Agrid;
Bseign=Transpose[Transpose[If[eN<1,λ ℓ δB (((1+(ℓ/(1-ℓ))((1-eC)/(1-eN))))^(ρ-1)) ((((1-ℓ)(1-eN)))^ρ)(u'[q1B]-1),0]/.gridresults]Bgrid];
findCNEa[a0_,b0_]:=Module[{stepA,stepB},
stepA[{ai_,bi_}]:=If[Abs[Max[Aseign[[;;,bi]]]-Aseign[[ai,bi]]]<10^(-prec),ai,Ordering[Aseign[[;;,bi]],-1][[1]]];
stepB[{ai_,bi_}]:=With[{anext=stepA[{ai,bi}]},
{anext,If[Abs[Max[Bseign[[anext,;;]]]-Bseign[[anext,bi]]]<10^(-prec),bi,Ordering[Bseign[[anext,;;]],-1][[1]]]}];
FixedPoint[stepB,{a0,b0},2prec(*,SameTest->(Abs[#1-#2]<10^(-prec)&)*)]];
findCNEb[a0_,b0_]:=Module[{stepA,stepB},
stepB[{ai_,bi_}]:=If[Abs[Max[Bseign[[ai,;;]]]-Bseign[[ai,bi]]]<10^(-prec),bi,Ordering[Bseign[[ai,;;]],-1][[1]]];
stepA[{ai_,bi_}]:=With[{bnext=stepB[{ai,bi}]},{If[Abs[Max[Aseign[[;;,bnext]]]-Aseign[[ai,bnext]]]<10^(-prec),ai,Ordering[Aseign[[;;,bnext]],-1][[1]]],bnext}];
FixedPoint[stepA,{a0,b0},2prec(*,SameTest->(Abs[#1-#2]<10^(-prec)&)*)]];
With[{asol=findCNEa[1,1],bsol=findCNEb[1,1]},GraphicsGrid[
{{MatrixPlot[Reverse@Transpose@Aseign,Frame->{{True,False},{True,False}},FrameLabel->{"B","A"},FrameTicks->None,BaseStyle->{FontSize->18,Bold},ImageSize->200,ColorFunction->(Lighter[Darker@Red,(2-2#)^(4-3.5#)]&),Epilog->{PointSize->Large,Darker@Green,Point[gridlength {Agrid[[bsol[[1]]]]-.5Agrid[[1]],Bgrid[[bsol[[2]]]]-.5Bgrid[[1]]}/{Abar,Bbar}],Darker@Blue,Point[gridlength {Agrid[[asol[[1]]]]-.5Agrid[[1]],Bgrid[[asol[[2]]]]-.5Bgrid[[1]]}/{Abar,Bbar}]}]},{MatrixPlot[Reverse@Transpose@Bseign,Frame->{{True,False},{True,False}},FrameLabel->{"B","A"},FrameTicks->None,BaseStyle->{FontSize->18,Bold},ImageSize->200,ColorFunction->(Lighter[Darker@Red,(2-2#)^(4-3.5#)]&),(*ColorFunction->ColorData[{"RedGreenSplit","Reverse"}],*)Epilog->{PointSize->Large,Darker@Green,Point[gridlength {Agrid[[bsol[[1]]]]-.5Agrid[[1]],Bgrid[[bsol[[2]]]]-.5Bgrid[[1]]}/{Abar,Bbar}],Darker@Blue,Point[gridlength {Agrid[[asol[[1]]]]-.5Agrid[[1]],Bgrid[[asol[[2]]]]-.5Bgrid[[1]]}/{Abar,Bbar}]}]}}]]]
Export["gamer00a10.pdf",%]

Comparing increasing degrees of a disadvantage for B (with CRS):

Module[{u=Log(*2Sqrt[#]&*),uinv=Exp,upinv,ℓ=.5,λ=.5,δA=1,δB=.99,ρ=0.0,M=1,i=.1,prec=6,gridlength=21,Abar,Bbar,Agrid,Bgrid,gridresults,Aseign,Bseign,ALP,findCNEa,findCNEb},
upinv=z/.Solve[u'[z]==#,z][[1]]&;
Abar=.98 M(1/upinv[1+i/(ℓ-λ ℓ(1-ℓ)δA)]-1);Bbar=Abar;(*M(1/upinv[1+i/(ℓ-λ ℓ(1-ℓ)δB)]-1);*)
Agrid=Range[gridlength]Abar/gridlength;Bgrid=Range[gridlength]Bbar/gridlength;
gridresults=Table[ge[{u,uinv,ℓ,λ,δA,δB,ρ,M,A,B,i,prec}],{A,Agrid},{B,Bgrid}];
Aseign=(If[eN>0,λ ℓ δA (((1+(ℓ/(1-ℓ))(eC/eN)))^(ρ-1)) ((((1-ℓ)eN))^ρ)(u'[q1A]-1),0]/.gridresults)Agrid;
Bseign=Transpose[Transpose[If[eN<1,λ ℓ δB (((1+(ℓ/(1-ℓ))((1-eC)/(1-eN))))^(ρ-1)) ((((1-ℓ)(1-eN)))^ρ)(u'[q1B]-1),0]/.gridresults]Bgrid];
findCNEa[a0_,b0_]:=Module[{stepA,stepB},
stepA[{ai_,bi_}]:=If[Abs[Max[Aseign[[;;,bi]]]-Aseign[[ai,bi]]]<10^(-prec),ai,Ordering[Aseign[[;;,bi]],-1][[1]]];
stepB[{ai_,bi_}]:=With[{anext=stepA[{ai,bi}]},
{anext,If[Abs[Max[Bseign[[anext,;;]]]-Bseign[[anext,bi]]]<10^(-prec),bi,Ordering[Bseign[[anext,;;]],-1][[1]]]}];
FixedPoint[stepB,{a0,b0},2prec(*,SameTest->(Abs[#1-#2]<10^(-prec)&)*)]];
findCNEb[a0_,b0_]:=Module[{stepA,stepB},
stepB[{ai_,bi_}]:=If[Abs[Max[Bseign[[ai,;;]]]-Bseign[[ai,bi]]]<10^(-prec),bi,Ordering[Bseign[[ai,;;]],-1][[1]]];
stepA[{ai_,bi_}]:=With[{bnext=stepB[{ai,bi}]},{If[Abs[Max[Aseign[[;;,bnext]]]-Aseign[[ai,bnext]]]<10^(-prec),ai,Ordering[Aseign[[;;,bnext]],-1][[1]]],bnext}];
FixedPoint[stepA,{a0,b0},2prec(*,SameTest->(Abs[#1-#2]<10^(-prec)&)*)]];
With[{asol=findCNEa[1,1],bsol=findCNEb[1,1]},GraphicsGrid[
{{MatrixPlot[Reverse@Transpose@Aseign,Frame->{{True,False},{True,False}},FrameLabel->{"B","A"},FrameTicks->None,BaseStyle->{FontSize->18,Bold},ImageSize->200,ColorFunction->(Lighter[Darker@Red,(2-2#)^(4-3.5#)]&),Epilog->{PointSize->Large,Darker@Green,Point[gridlength {Agrid[[bsol[[1]]]]-.5Agrid[[1]],Bgrid[[bsol[[2]]]]-.5Bgrid[[1]]}/{Abar,Bbar}],Darker@Blue,Point[gridlength {Agrid[[asol[[1]]]]-.5Agrid[[1]],Bgrid[[asol[[2]]]]-.5Bgrid[[1]]}/{Abar,Bbar}]}]},{MatrixPlot[Reverse@Transpose@Bseign,Frame->{{True,False},{True,False}},FrameLabel->{"B","A"},FrameTicks->None,BaseStyle->{FontSize->18,Bold},ImageSize->200,ColorFunction->(Lighter[Darker@Red,(2-2#)^(4-3.5#)]&),Epilog->{PointSize->Large,Darker@Green,Point[gridlength {Agrid[[bsol[[1]]]]-.5Agrid[[1]],Bgrid[[bsol[[2]]]]-.5Bgrid[[1]]}/{Abar,Bbar}],Darker@Blue,Point[gridlength {Agrid[[asol[[1]]]]-.5Agrid[[1]],Bgrid[[asol[[2]]]]-.5Bgrid[[1]]}/{Abar,Bbar}]}]}}]]]
Export["gamer00a99.pdf",%]

Module[{u=Log(*2Sqrt[#]&*),uinv=Exp,upinv,ℓ=.5,λ=.5,δA=1,δB=.9,ρ=0.0,M=1,i=.1,prec=6,gridlength=21,Abar,Bbar,Agrid,Bgrid,gridresults,Aseign,Bseign,ALP,findCNEa,findCNEb},
upinv=z/.Solve[u'[z]==#,z][[1]]&;
Abar=.98 M(1/upinv[1+i/(ℓ-λ ℓ(1-ℓ)δA)]-1);Bbar=Abar;(*M(1/upinv[1+i/(ℓ-λ ℓ(1-ℓ)δB)]-1);*)
Agrid=Range[gridlength]Abar/gridlength;Bgrid=Range[gridlength]Bbar/gridlength;
gridresults=Table[ge[{u,uinv,ℓ,λ,δA,δB,ρ,M,A,B,i,prec}],{A,Agrid},{B,Bgrid}];
Aseign=(If[eN>0,λ ℓ δA (((1+(ℓ/(1-ℓ))(eC/eN)))^(ρ-1)) ((((1-ℓ)eN))^ρ)(u'[q1A]-1),0]/.gridresults)Agrid;
Bseign=Transpose[Transpose[If[eN<1,λ ℓ δB (((1+(ℓ/(1-ℓ))((1-eC)/(1-eN))))^(ρ-1)) ((((1-ℓ)(1-eN)))^ρ)(u'[q1B]-1),0]/.gridresults]Bgrid];
findCNEa[a0_,b0_]:=Module[{stepA,stepB},
stepA[{ai_,bi_}]:=If[Abs[Max[Aseign[[;;,bi]]]-Aseign[[ai,bi]]]<10^(-prec),ai,Ordering[Aseign[[;;,bi]],-1][[1]]];
stepB[{ai_,bi_}]:=With[{anext=stepA[{ai,bi}]},
{anext,If[Abs[Max[Bseign[[anext,;;]]]-Bseign[[anext,bi]]]<10^(-prec),bi,Ordering[Bseign[[anext,;;]],-1][[1]]]}];
FixedPoint[stepB,{a0,b0},2prec(*,SameTest->(Abs[#1-#2]<10^(-prec)&)*)]];
findCNEb[a0_,b0_]:=Module[{stepA,stepB},
stepB[{ai_,bi_}]:=If[Abs[Max[Bseign[[ai,;;]]]-Bseign[[ai,bi]]]<10^(-prec),bi,Ordering[Bseign[[ai,;;]],-1][[1]]];
stepA[{ai_,bi_}]:=With[{bnext=stepB[{ai,bi}]},{If[Abs[Max[Aseign[[;;,bnext]]]-Aseign[[ai,bnext]]]<10^(-prec),ai,Ordering[Aseign[[;;,bnext]],-1][[1]]],bnext}];
FixedPoint[stepA,{a0,b0},2prec(*,SameTest->(Abs[#1-#2]<10^(-prec)&)*)]];
With[{asol=findCNEa[1,1],bsol=findCNEb[1,1]},GraphicsGrid[
{{MatrixPlot[Reverse@Transpose@Aseign,Frame->{{True,False},{True,False}},FrameLabel->{"B","A"},FrameTicks->None,BaseStyle->{FontSize->18,Bold},ImageSize->200,ColorFunction->(Lighter[Darker@Red,(2-2#)^(4-3.5#)]&),Epilog->{PointSize->Large,Darker@Green,Point[gridlength {Agrid[[bsol[[1]]]]-.5Agrid[[1]],Bgrid[[bsol[[2]]]]-.5Bgrid[[1]]}/{Abar,Bbar}],Darker@Blue,Point[gridlength {Agrid[[asol[[1]]]]-.5Agrid[[1]],Bgrid[[asol[[2]]]]-.5Bgrid[[1]]}/{Abar,Bbar}]}]},{MatrixPlot[Reverse@Transpose@Bseign,Frame->{{True,False},{True,False}},FrameLabel->{"B","A"},FrameTicks->None,BaseStyle->{FontSize->18,Bold},ImageSize->200,ColorFunction->(Lighter[Darker@Red,(2-2#)^(4-3.5#)]&),Epilog->{PointSize->Large,Darker@Green,Point[gridlength {Agrid[[bsol[[1]]]]-.5Agrid[[1]],Bgrid[[bsol[[2]]]]-.5Bgrid[[1]]}/{Abar,Bbar}],Darker@Blue,Point[gridlength {Agrid[[asol[[1]]]]-.5Agrid[[1]],Bgrid[[asol[[2]]]]-.5Bgrid[[1]]}/{Abar,Bbar}]}]}}]]]
Export["gamer00a90.pdf",%]

Module[{u=Log(*2Sqrt[#]&*),uinv=Exp,upinv,ℓ=.5,λ=.5,δA=1,δB=.6,ρ=0.0,M=1,i=.1,prec=6,gridlength=21,Abar,Bbar,Agrid,Bgrid,gridresults,Aseign,Bseign,ALP,findCNEa,findCNEb},
upinv=z/.Solve[u'[z]==#,z][[1]]&;
Abar=.98 M(1/upinv[1+i/(ℓ-λ ℓ(1-ℓ)δA)]-1);Bbar=Abar;(*M(1/upinv[1+i/(ℓ-λ ℓ(1-ℓ)δB)]-1);*)
Agrid=Range[gridlength]Abar/gridlength;Bgrid=Range[gridlength]Bbar/gridlength;
gridresults=Table[ge[{u,uinv,ℓ,λ,δA,δB,ρ,M,A,B,i,prec}],{A,Agrid},{B,Bgrid}];
Aseign=(If[eN>0,λ ℓ δA (((1+(ℓ/(1-ℓ))(eC/eN)))^(ρ-1)) ((((1-ℓ)eN))^ρ)(u'[q1A]-1),0]/.gridresults)Agrid;
Bseign=Transpose[Transpose[If[eN<1,λ ℓ δB (((1+(ℓ/(1-ℓ))((1-eC)/(1-eN))))^(ρ-1)) ((((1-ℓ)(1-eN)))^ρ)(u'[q1B]-1),0]/.gridresults]Bgrid];
findCNEa[a0_,b0_]:=Module[{stepA,stepB},
stepA[{ai_,bi_}]:=If[Abs[Max[Aseign[[;;,bi]]]-Aseign[[ai,bi]]]<10^(-prec),ai,Ordering[Aseign[[;;,bi]],-1][[1]]];
stepB[{ai_,bi_}]:=With[{anext=stepA[{ai,bi}]},
{anext,If[Abs[Max[Bseign[[anext,;;]]]-Bseign[[anext,bi]]]<10^(-prec),bi,Ordering[Bseign[[anext,;;]],-1][[1]]]}];
FixedPoint[stepB,{a0,b0},2prec(*,SameTest->(Abs[#1-#2]<10^(-prec)&)*)]];
findCNEb[a0_,b0_]:=Module[{stepA,stepB},
stepB[{ai_,bi_}]:=If[Abs[Max[Bseign[[ai,;;]]]-Bseign[[ai,bi]]]<10^(-prec),bi,Ordering[Bseign[[ai,;;]],-1][[1]]];
stepA[{ai_,bi_}]:=With[{bnext=stepB[{ai,bi}]},{If[Abs[Max[Aseign[[;;,bnext]]]-Aseign[[ai,bnext]]]<10^(-prec),ai,Ordering[Aseign[[;;,bnext]],-1][[1]]],bnext}];
FixedPoint[stepA,{a0,b0},2prec(*,SameTest->(Abs[#1-#2]<10^(-prec)&)*)]];
With[{asol=findCNEa[1,1],bsol=findCNEb[1,1]},GraphicsGrid[
{{MatrixPlot[Reverse@Transpose@Aseign,Frame->{{True,False},{True,False}},FrameLabel->{"B","A"},FrameTicks->None,BaseStyle->{FontSize->18,Bold},ImageSize->200,ColorFunction->(Lighter[Darker@Red,(2-2#)^(4-3.5#)]&),Epilog->{PointSize->Large,Darker@Green,Point[gridlength {Agrid[[bsol[[1]]]]-.5Agrid[[1]],Bgrid[[bsol[[2]]]]-.5Bgrid[[1]]}/{Abar,Bbar}],Darker@Blue,Point[gridlength {Agrid[[asol[[1]]]]-.5Agrid[[1]],Bgrid[[asol[[2]]]]-.5Bgrid[[1]]}/{Abar,Bbar}]}]},{MatrixPlot[Reverse@Transpose@Bseign,Frame->{{True,False},{True,False}},FrameLabel->{"B","A"},FrameTicks->None,BaseStyle->{FontSize->18,Bold},ImageSize->200,ColorFunction->(Lighter[Darker@Red,(2-2#)^(4-3.5#)]&),Epilog->{PointSize->Large,Darker@Green,Point[gridlength {Agrid[[bsol[[1]]]]-.5Agrid[[1]],Bgrid[[bsol[[2]]]]-.5Bgrid[[1]]}/{Abar,Bbar}],Darker@Blue,Point[gridlength {Agrid[[asol[[1]]]]-.5Agrid[[1]],Bgrid[[asol[[2]]]]-.5Bgrid[[1]]}/{Abar,Bbar}]}]}}]]]
Export["gamer00a60.pdf",%]

Comparing increasing degrees of IRS (with balance):

Module[{u=Log(*2Sqrt[#]&*),uinv=Exp,upinv,ℓ=.5,λ=.5,δA=1,δB=1,ρ=0.01,M=1,i=.1,prec=13,gridlength=21,Abar,Bbar,Agrid,Bgrid,gridresults,Aseign,Bseign,ALP,findCNEa,findCNEb},
upinv=z/.Solve[u'[z]==#,z][[1]]&;
Abar=.98 M(1/upinv[1+i/(ℓ-λ ℓ(1-ℓ)δA)]-1);Bbar=Abar;(*M(1/upinv[1+i/(ℓ-λ ℓ(1-ℓ)δB)]-1);*)
Agrid=Range[gridlength]Abar/gridlength;Bgrid=Range[gridlength]Bbar/gridlength;
gridresults=Table[ge[{u,uinv,ℓ,λ,δA,δB,ρ,M,A,B,i,prec}],{A,Agrid},{B,Bgrid}];
Aseign=(If[eN>0,λ ℓ δA (((1+(ℓ/(1-ℓ))(eC/eN)))^(ρ-1)) ((((1-ℓ)eN))^ρ)(u'[q1A]-1),0]/.gridresults)Agrid;
Bseign=Transpose[Transpose[If[eN<1,λ ℓ δB (((1+(ℓ/(1-ℓ))((1-eC)/(1-eN))))^(ρ-1)) ((((1-ℓ)(1-eN)))^ρ)(u'[q1B]-1),0]/.gridresults]Bgrid];
findCNEa[a0_,b0_]:=Module[{stepA,stepB},
stepA[{ai_,bi_}]:=If[Abs[Max[Aseign[[;;,bi]]]-Aseign[[ai,bi]]]<10^(-prec),ai,Ordering[Aseign[[;;,bi]],-1][[1]]];
stepB[{ai_,bi_}]:=With[{anext=stepA[{ai,bi}]},
{anext,If[Abs[Max[Bseign[[anext,;;]]]-Bseign[[anext,bi]]]<10^(-prec),bi,Ordering[Bseign[[anext,;;]],-1][[1]]]}];
FixedPoint[stepB,{a0,b0},2prec(*,SameTest->(Abs[#1-#2]<10^(-prec)&)*)]];
findCNEb[a0_,b0_]:=Module[{stepA,stepB},
stepB[{ai_,bi_}]:=If[Abs[Max[Bseign[[ai,;;]]]-Bseign[[ai,bi]]]<10^(-prec),bi,Ordering[Bseign[[ai,;;]],-1][[1]]];
stepA[{ai_,bi_}]:=With[{bnext=stepB[{ai,bi}]},{If[Abs[Max[Aseign[[;;,bnext]]]-Aseign[[ai,bnext]]]<10^(-prec),ai,Ordering[Aseign[[;;,bnext]],-1][[1]]],bnext}];
FixedPoint[stepA,{a0,b0},2prec(*,SameTest->(Abs[#1-#2]<10^(-prec)&)*)]];
With[{asol=findCNEa[1,1],bsol=findCNEb[1,1]},GraphicsGrid[
{{MatrixPlot[Reverse@Transpose@Aseign,Frame->{{True,False},{True,False}},FrameLabel->{"B","A"},FrameTicks->None,BaseStyle->{FontSize->18,Bold},ImageSize->200,ColorFunction->(Lighter[Darker@Red,(2-2#)^(4-3.5#)]&),Epilog->{PointSize->Large,Darker@Green,Point[gridlength {Agrid[[bsol[[1]]]]-.5Agrid[[1]],Bgrid[[bsol[[2]]]]-.5Bgrid[[1]]}/{Abar,Bbar}],Darker@Blue,Point[gridlength {Agrid[[asol[[1]]]]-.5Agrid[[1]],Bgrid[[asol[[2]]]]-.5Bgrid[[1]]}/{Abar,Bbar}]}]},{MatrixPlot[Reverse@Transpose@Bseign,Frame->{{True,False},{True,False}},FrameLabel->{"B","A"},FrameTicks->None,BaseStyle->{FontSize->18,Bold},ImageSize->200,ColorFunction->(Lighter[Darker@Red,(2-2#)^(4-3.5#)]&),Epilog->{PointSize->Large,Darker@Green,Point[gridlength {Agrid[[bsol[[1]]]]-.5Agrid[[1]],Bgrid[[bsol[[2]]]]-.5Bgrid[[1]]}/{Abar,Bbar}],Darker@Blue,Point[gridlength {Agrid[[asol[[1]]]]-.5Agrid[[1]],Bgrid[[asol[[2]]]]-.5Bgrid[[1]]}/{Abar,Bbar}]}]}}]]]
Export["gamer01a10.pdf",%]

Module[{u=Log(*2Sqrt[#]&*),uinv=Exp,upinv,ℓ=.5,λ=.5,δA=1,δB=1,ρ=0.05,M=1,i=.1,prec=13,gridlength=21,Abar,Bbar,Agrid,Bgrid,gridresults,Aseign,Bseign,ALP,findCNEa,findCNEb},
upinv=z/.Solve[u'[z]==#,z][[1]]&;
Abar=.98 M(1/upinv[1+i/(ℓ-λ ℓ(1-ℓ)δA)]-1);Bbar=Abar;(*M(1/upinv[1+i/(ℓ-λ ℓ(1-ℓ)δB)]-1);*)
Agrid=Range[gridlength]Abar/gridlength;Bgrid=Range[gridlength]Bbar/gridlength;
gridresults=Table[ge[{u,uinv,ℓ,λ,δA,δB,ρ,M,A,B,i,prec}],{A,Agrid},{B,Bgrid}];
Aseign=(If[eN>0,λ ℓ δA (((1+(ℓ/(1-ℓ))(eC/eN)))^(ρ-1)) ((((1-ℓ)eN))^ρ)(u'[q1A]-1),0]/.gridresults)Agrid;
Bseign=Transpose[Transpose[If[eN<1,λ ℓ δB (((1+(ℓ/(1-ℓ))((1-eC)/(1-eN))))^(ρ-1)) ((((1-ℓ)(1-eN)))^ρ)(u'[q1B]-1),0]/.gridresults]Bgrid];
findCNEa[a0_,b0_]:=Module[{stepA,stepB},
stepA[{ai_,bi_}]:=If[Abs[Max[Aseign[[;;,bi]]]-Aseign[[ai,bi]]]<10^(-prec),ai,Ordering[Aseign[[;;,bi]],-1][[1]]];
stepB[{ai_,bi_}]:=With[{anext=stepA[{ai,bi}]},
{anext,If[Abs[Max[Bseign[[anext,;;]]]-Bseign[[anext,bi]]]<10^(-prec),bi,Ordering[Bseign[[anext,;;]],-1][[1]]]}];
FixedPoint[stepB,{a0,b0},2prec(*,SameTest->(Abs[#1-#2]<10^(-prec)&)*)]];
findCNEb[a0_,b0_]:=Module[{stepA,stepB},
stepB[{ai_,bi_}]:=If[Abs[Max[Bseign[[ai,;;]]]-Bseign[[ai,bi]]]<10^(-prec),bi,Ordering[Bseign[[ai,;;]],-1][[1]]];
stepA[{ai_,bi_}]:=With[{bnext=stepB[{ai,bi}]},{If[Abs[Max[Aseign[[;;,bnext]]]-Aseign[[ai,bnext]]]<10^(-prec),ai,Ordering[Aseign[[;;,bnext]],-1][[1]]],bnext}];
FixedPoint[stepA,{a0,b0},2prec(*,SameTest->(Abs[#1-#2]<10^(-prec)&)*)]];
With[{asol=findCNEa[1,1],bsol=findCNEb[1,1]},GraphicsGrid[
{{MatrixPlot[Reverse@Transpose@Aseign,Frame->{{True,False},{True,False}},FrameLabel->{"B","A"},FrameTicks->None,BaseStyle->{FontSize->18,Bold},ImageSize->200,ColorFunction->(Lighter[Darker@Red,(2-2#)^(4-3.5#)]&),Epilog->{PointSize->Large,Darker@Green,Point[gridlength {Agrid[[bsol[[1]]]]-.5Agrid[[1]],Bgrid[[bsol[[2]]]]-.5Bgrid[[1]]}/{Abar,Bbar}],Darker@Blue,Point[gridlength {Agrid[[asol[[1]]]]-.5Agrid[[1]],Bgrid[[asol[[2]]]]-.5Bgrid[[1]]}/{Abar,Bbar}]}]},{MatrixPlot[Reverse@Transpose@Bseign,Frame->{{True,False},{True,False}},FrameLabel->{"B","A"},FrameTicks->None,BaseStyle->{FontSize->18,Bold},ImageSize->200,ColorFunction->(Lighter[Darker@Red,(2-2#)^(4-3.5#)]&),Epilog->{PointSize->Large,Darker@Green,Point[gridlength {Agrid[[bsol[[1]]]]-.5Agrid[[1]],Bgrid[[bsol[[2]]]]-.5Bgrid[[1]]}/{Abar,Bbar}],Darker@Blue,Point[gridlength {Agrid[[asol[[1]]]]-.5Agrid[[1]],Bgrid[[asol[[2]]]]-.5Bgrid[[1]]}/{Abar,Bbar}]}]}}]]]
Export["gamer05a10.pdf",%]

Module[{u=Log(*2Sqrt[#]&*),uinv=Exp,upinv,ℓ=.5,λ=.5,δA=1,δB=1,ρ=0.5,M=1,i=.1,prec=13,gridlength=21,Abar,Bbar,Agrid,Bgrid,gridresults,Aseign,Bseign,ALP,findCNEa,findCNEb},
upinv=z/.Solve[u'[z]==#,z][[1]]&;
Abar=.98 M(1/upinv[1+i/(ℓ-λ ℓ(1-ℓ)δA)]-1);Bbar=Abar;(*M(1/upinv[1+i/(ℓ-λ ℓ(1-ℓ)δB)]-1);*)
Agrid=Range[gridlength]Abar/gridlength;Bgrid=Range[gridlength]Bbar/gridlength;
gridresults=Table[ge[{u,uinv,ℓ,λ,δA,δB,ρ,M,A,B,i,prec}],{A,Agrid},{B,Bgrid}];
Aseign=(If[eN>0,λ ℓ δA (((1+(ℓ/(1-ℓ))(eC/eN)))^(ρ-1)) ((((1-ℓ)eN))^ρ)(u'[q1A]-1),0]/.gridresults)Agrid;
Bseign=Transpose[Transpose[If[eN<1,λ ℓ δB (((1+(ℓ/(1-ℓ))((1-eC)/(1-eN))))^(ρ-1)) ((((1-ℓ)(1-eN)))^ρ)(u'[q1B]-1),0]/.gridresults]Bgrid];
findCNEa[a0_,b0_]:=Module[{stepA,stepB},
stepA[{ai_,bi_}]:=If[Abs[Max[Aseign[[;;,bi]]]-Aseign[[ai,bi]]]<10^(-prec),ai,Ordering[Aseign[[;;,bi]],-1][[1]]];
stepB[{ai_,bi_}]:=With[{anext=stepA[{ai,bi}]},
{anext,If[Abs[Max[Bseign[[anext,;;]]]-Bseign[[anext,bi]]]<10^(-prec),bi,Ordering[Bseign[[anext,;;]],-1][[1]]]}];
FixedPoint[stepB,{a0,b0},2prec(*,SameTest->(Abs[#1-#2]<10^(-prec)&)*)]];
findCNEb[a0_,b0_]:=Module[{stepA,stepB},
stepB[{ai_,bi_}]:=If[Abs[Max[Bseign[[ai,;;]]]-Bseign[[ai,bi]]]<10^(-prec),bi,Ordering[Bseign[[ai,;;]],-1][[1]]];
stepA[{ai_,bi_}]:=With[{bnext=stepB[{ai,bi}]},{If[Abs[Max[Aseign[[;;,bnext]]]-Aseign[[ai,bnext]]]<10^(-prec),ai,Ordering[Aseign[[;;,bnext]],-1][[1]]],bnext}];
FixedPoint[stepA,{a0,b0},2prec(*,SameTest->(Abs[#1-#2]<10^(-prec)&)*)]];
With[{asol=findCNEa[1,1],bsol=findCNEb[1,1]},GraphicsGrid[
{{MatrixPlot[Reverse@Transpose@Aseign,Frame->{{True,False},{True,False}},FrameLabel->{"B","A"},FrameTicks->None,BaseStyle->{FontSize->18,Bold},ImageSize->200,ColorFunction->(Lighter[Darker@Red,(2-2#)^(4-3.5#)]&),Epilog->{PointSize->Large,Darker@Green,Point[gridlength {Agrid[[bsol[[1]]]]-.5Agrid[[1]],Bgrid[[bsol[[2]]]]-.5Bgrid[[1]]}/{Abar,Bbar}],Darker@Blue,Point[gridlength {Agrid[[asol[[1]]]]-.5Agrid[[1]],Bgrid[[asol[[2]]]]-.5Bgrid[[1]]}/{Abar,Bbar}]}]},{MatrixPlot[Reverse@Transpose@Bseign,Frame->{{True,False},{True,False}},FrameLabel->{"B","A"},FrameTicks->None,BaseStyle->{FontSize->18,Bold},ImageSize->200,ColorFunction->(Lighter[Darker@Red,(2-2#)^(4-3.5#)]&),Epilog->{PointSize->Large,Darker@Green,Point[gridlength {Agrid[[bsol[[1]]]]-.5Agrid[[1]],Bgrid[[bsol[[2]]]]-.5Bgrid[[1]]}/{Abar,Bbar}],Darker@Blue,Point[gridlength {Agrid[[asol[[1]]]]-.5Agrid[[1]],Bgrid[[asol[[2]]]]-.5Bgrid[[1]]}/{Abar,Bbar}]}]}}]]]
Export["gamer50a10.pdf",%]

Finally, small IRS and small B-disadvantage:

Module[{u=Log(*2Sqrt[#]&*),uinv=Exp,upinv,ℓ=.5,λ=.5,δA=1,δB=.99,ρ=0.01,M=1,i=.1,prec=13,gridlength=21,Abar,Bbar,Agrid,Bgrid,gridresults,Aseign,Bseign,ALP,findCNEa,findCNEb},
upinv=z/.Solve[u'[z]==#,z][[1]]&;
Abar=.98 M(1/upinv[1+i/(ℓ-λ ℓ(1-ℓ)δA)]-1);Bbar=Abar;(*M(1/upinv[1+i/(ℓ-λ ℓ(1-ℓ)δB)]-1);*)
Agrid=Range[gridlength]Abar/gridlength;Bgrid=Range[gridlength]Bbar/gridlength;
gridresults=Table[ge[{u,uinv,ℓ,λ,δA,δB,ρ,M,A,B,i,prec}],{A,Agrid},{B,Bgrid}];
Aseign=(If[eN>0,λ ℓ δA (((1+(ℓ/(1-ℓ))(eC/eN)))^(ρ-1)) ((((1-ℓ)eN))^ρ)(u'[q1A]-1),0]/.gridresults)Agrid;
Bseign=Transpose[Transpose[If[eN<1,λ ℓ δB (((1+(ℓ/(1-ℓ))((1-eC)/(1-eN))))^(ρ-1)) ((((1-ℓ)(1-eN)))^ρ)(u'[q1B]-1),0]/.gridresults]Bgrid];
findCNEa[a0_,b0_]:=Module[{stepA,stepB},
stepA[{ai_,bi_}]:=If[Abs[Max[Aseign[[;;,bi]]]-Aseign[[ai,bi]]]<10^(-prec),ai,Ordering[Aseign[[;;,bi]],-1][[1]]];
stepB[{ai_,bi_}]:=With[{anext=stepA[{ai,bi}]},
{anext,If[Abs[Max[Bseign[[anext,;;]]]-Bseign[[anext,bi]]]<10^(-prec),bi,Ordering[Bseign[[anext,;;]],-1][[1]]]}];
FixedPoint[stepB,{a0,b0},2prec(*,SameTest->(Abs[#1-#2]<10^(-prec)&)*)]];
findCNEb[a0_,b0_]:=Module[{stepA,stepB},
stepB[{ai_,bi_}]:=If[Abs[Max[Bseign[[ai,;;]]]-Bseign[[ai,bi]]]<10^(-prec),bi,Ordering[Bseign[[ai,;;]],-1][[1]]];
stepA[{ai_,bi_}]:=With[{bnext=stepB[{ai,bi}]},{If[Abs[Max[Aseign[[;;,bnext]]]-Aseign[[ai,bnext]]]<10^(-prec),ai,Ordering[Aseign[[;;,bnext]],-1][[1]]],bnext}];
FixedPoint[stepA,{a0,b0},2prec(*,SameTest->(Abs[#1-#2]<10^(-prec)&)*)]];
With[{asol=findCNEa[1,1],bsol=findCNEb[1,1]},GraphicsGrid[
{{MatrixPlot[Reverse@Transpose@Aseign,Frame->{{True,False},{True,False}},FrameLabel->{"B","A"},FrameTicks->None,BaseStyle->{FontSize->18,Bold},ImageSize->200,ColorFunction->(Lighter[Darker@Red,(2-2#)^(4-3.5#)]&),Epilog->{PointSize->Large,Darker@Green,Point[gridlength {Agrid[[bsol[[1]]]]-.5Agrid[[1]],Bgrid[[bsol[[2]]]]-.5Bgrid[[1]]}/{Abar,Bbar}],Darker@Blue,Point[gridlength {Agrid[[asol[[1]]]]-.5Agrid[[1]],Bgrid[[asol[[2]]]]-.5Bgrid[[1]]}/{Abar,Bbar}]}]},{MatrixPlot[Reverse@Transpose@Bseign,Frame->{{True,False},{True,False}},FrameLabel->{"B","A"},FrameTicks->None,BaseStyle->{FontSize->18,Bold},ImageSize->200,ColorFunction->(Lighter[Darker@Red,(2-2#)^(4-3.5#)]&),Epilog->{PointSize->Large,Darker@Green,Point[gridlength {Agrid[[bsol[[1]]]]-.5Agrid[[1]],Bgrid[[bsol[[2]]]]-.5Bgrid[[1]]}/{Abar,Bbar}],Darker@Blue,Point[gridlength {Agrid[[asol[[1]]]]-.5Agrid[[1]],Bgrid[[asol[[2]]]]-.5Bgrid[[1]]}/{Abar,Bbar}]}]}}]]]
Export["gamer01a99.pdf",%]

Module[{u=Log(*2Sqrt[#]&*),uinv=Exp,upinv,ℓ=.5,λ=.5,δA=1,δB=.98,ρ=0.02,M=1,i=.1,prec=13,gridlength=21,Abar,Bbar,Agrid,Bgrid,gridresults,Aseign,Bseign,ALP,findCNEa,findCNEb},
upinv=z/.Solve[u'[z]==#,z][[1]]&;
Abar=.98 M(1/upinv[1+i/(ℓ-λ ℓ(1-ℓ)δA)]-1);Bbar=Abar;(*M(1/upinv[1+i/(ℓ-λ ℓ(1-ℓ)δB)]-1);*)
Agrid=Range[gridlength]Abar/gridlength;Bgrid=Range[gridlength]Bbar/gridlength;
gridresults=Table[ge[{u,uinv,ℓ,λ,δA,δB,ρ,M,A,B,i,prec}],{A,Agrid},{B,Bgrid}];
Aseign=(If[eN>0,λ ℓ δA (((1+(ℓ/(1-ℓ))(eC/eN)))^(ρ-1)) ((((1-ℓ)eN))^ρ)(u'[q1A]-1),0]/.gridresults)Agrid;
Bseign=Transpose[Transpose[If[eN<1,λ ℓ δB (((1+(ℓ/(1-ℓ))((1-eC)/(1-eN))))^(ρ-1)) ((((1-ℓ)(1-eN)))^ρ)(u'[q1B]-1),0]/.gridresults]Bgrid];
findCNEa[a0_,b0_]:=Module[{stepA,stepB},
stepA[{ai_,bi_}]:=If[Abs[Max[Aseign[[;;,bi]]]-Aseign[[ai,bi]]]<10^(-prec),ai,Ordering[Aseign[[;;,bi]],-1][[1]]];
stepB[{ai_,bi_}]:=With[{anext=stepA[{ai,bi}]},
{anext,If[Abs[Max[Bseign[[anext,;;]]]-Bseign[[anext,bi]]]<10^(-prec),bi,Ordering[Bseign[[anext,;;]],-1][[1]]]}];
FixedPoint[stepB,{a0,b0},2prec(*,SameTest->(Abs[#1-#2]<10^(-prec)&)*)]];
findCNEb[a0_,b0_]:=Module[{stepA,stepB},
stepB[{ai_,bi_}]:=If[Abs[Max[Bseign[[ai,;;]]]-Bseign[[ai,bi]]]<10^(-prec),bi,Ordering[Bseign[[ai,;;]],-1][[1]]];
stepA[{ai_,bi_}]:=With[{bnext=stepB[{ai,bi}]},{If[Abs[Max[Aseign[[;;,bnext]]]-Aseign[[ai,bnext]]]<10^(-prec),ai,Ordering[Aseign[[;;,bnext]],-1][[1]]],bnext}];
FixedPoint[stepA,{a0,b0},2prec(*,SameTest->(Abs[#1-#2]<10^(-prec)&)*)]];
With[{asol=findCNEa[1,1],bsol=findCNEb[1,1]},GraphicsGrid[
{{MatrixPlot[Reverse@Transpose@Aseign,Frame->{{True,False},{True,False}},FrameLabel->{"B","A"},FrameTicks->None,BaseStyle->{FontSize->18,Bold},ImageSize->200,ColorFunction->(Lighter[Darker@Red,(2-2#)^(4-3.5#)]&),Epilog->{PointSize->Large,Darker@Green,Point[gridlength {Agrid[[bsol[[1]]]]-.5Agrid[[1]],Bgrid[[bsol[[2]]]]-.5Bgrid[[1]]}/{Abar,Bbar}],Darker@Blue,Point[gridlength {Agrid[[asol[[1]]]]-.5Agrid[[1]],Bgrid[[asol[[2]]]]-.5Bgrid[[1]]}/{Abar,Bbar}]}]},{MatrixPlot[Reverse@Transpose@Bseign,Frame->{{True,False},{True,False}},FrameLabel->{"B","A"},FrameTicks->None,BaseStyle->{FontSize->18,Bold},ImageSize->200,ColorFunction->(Lighter[Darker@Red,(2-2#)^(4-3.5#)]&),Epilog->{PointSize->Large,Darker@Green,Point[gridlength {Agrid[[bsol[[1]]]]-.5Agrid[[1]],Bgrid[[bsol[[2]]]]-.5Bgrid[[1]]}/{Abar,Bbar}],Darker@Blue,Point[gridlength {Agrid[[asol[[1]]]]-.5Agrid[[1]],Bgrid[[asol[[2]]]]-.5Bgrid[[1]]}/{Abar,Bbar}]}]}}]]]
Export["gamer02a98.pdf",%]

This one is for experimentation, with lower precision:

Module[{u=Log(*2Sqrt[#]&*),uinv=Exp,upinv,ℓ=.9,λ=.5,δA=1,δB=.98,ρ=0.0,M=1,i=.1,prec=13,gridlength=15,Abar,Bbar,Agrid,Bgrid,gridresults,Aseign,Bseign,ALP,findCNEa,findCNEb},
upinv=z/.Solve[u'[z]==#,z][[1]]&;
Abar=M(1/upinv[1+i/(ℓ-λ ℓ(1-ℓ)δA)]-1);Bbar=Abar;(*M(1/upinv[1+i/(ℓ-λ ℓ(1-ℓ)δB)]-1);*)
Agrid=Range[gridlength]Abar/gridlength;Bgrid=Range[gridlength]Bbar/gridlength;
gridresults=Table[ge[{u,uinv,ℓ,λ,δA,δB,ρ,M,A,B,i,prec}],{A,Agrid},{B,Bgrid}];
Aseign=(If[eN>0,λ ℓ δA (((1+(ℓ/(1-ℓ))(eC/eN)))^(ρ-1)) ((((1-ℓ)eN))^ρ)(u'[q1A]-1),0]/.gridresults)Agrid;
Bseign=Transpose[Transpose[If[eN<1,λ ℓ δB (((1+(ℓ/(1-ℓ))((1-eC)/(1-eN))))^(ρ-1)) ((((1-ℓ)(1-eN)))^ρ)(u'[q1B]-1),0]/.gridresults]Bgrid];
findCNEa[a0_,b0_]:=Module[{stepA,stepB},
stepA[{ai_,bi_}]:=If[Abs[Max[Aseign[[;;,bi]]]-Aseign[[ai,bi]]]<10^(-prec),ai,Ordering[Aseign[[;;,bi]],-1][[1]]];
stepB[{ai_,bi_}]:=With[{anext=stepA[{ai,bi}]},
{anext,If[Abs[Max[Bseign[[anext,;;]]]-Bseign[[anext,bi]]]<10^(-prec),bi,Ordering[Bseign[[anext,;;]],-1][[1]]]}];
FixedPoint[stepB,{a0,b0},2prec(*,SameTest->(Abs[#1-#2]<10^(-prec)&)*)]];
findCNEb[a0_,b0_]:=Module[{stepA,stepB},
stepB[{ai_,bi_}]:=If[Abs[Max[Bseign[[ai,;;]]]-Bseign[[ai,bi]]]<10^(-prec),bi,Ordering[Bseign[[ai,;;]],-1][[1]]];
stepA[{ai_,bi_}]:=With[{bnext=stepB[{ai,bi}]},{If[Abs[Max[Aseign[[;;,bnext]]]-Aseign[[ai,bnext]]]<10^(-prec),ai,Ordering[Aseign[[;;,bnext]],-1][[1]]],bnext}];
FixedPoint[stepA,{a0,b0},2prec(*,SameTest->(Abs[#1-#2]<10^(-prec)&)*)]];
With[{asol=findCNEa[1,1],bsol=findCNEb[1,1]},GraphicsGrid[
{{MatrixPlot[Reverse@Transpose@Aseign,Frame->{{True,False},{True,False}},FrameLabel->{"B","A"},FrameTicks->None,BaseStyle->{FontSize->18,Bold},ImageSize->200,ColorFunction->(Lighter[Darker@Red,(2-2#)^(4-3.5#)]&),Epilog->{PointSize->Large,Darker@Green,Point[gridlength {Agrid[[bsol[[1]]]]-.5Agrid[[1]],Bgrid[[bsol[[2]]]]-.5Bgrid[[1]]}/{Abar,Bbar}],Darker@Blue,Point[gridlength {Agrid[[asol[[1]]]]-.5Agrid[[1]],Bgrid[[asol[[2]]]]-.5Bgrid[[1]]}/{Abar,Bbar}]}]},{MatrixPlot[Reverse@Transpose@Bseign,Frame->{{True,False},{True,False}},FrameLabel->{"B","A"},FrameTicks->None,BaseStyle->{FontSize->18,Bold},ImageSize->200,ColorFunction->(Lighter[Darker@Red,(2-2#)^(4-3.5#)]&),Epilog->{PointSize->Large,Darker@Green,Point[gridlength {Agrid[[bsol[[1]]]]-.5Agrid[[1]],Bgrid[[bsol[[2]]]]-.5Bgrid[[1]]}/{Abar,Bbar}],Darker@Blue,Point[gridlength {Agrid[[asol[[1]]]]-.5Agrid[[1]],Bgrid[[asol[[2]]]]-.5Bgrid[[1]]}/{Abar,Bbar}]}]}}]]]

Analyzing the strategic equilibria

Algorithm to find strategic equilibria

Clear[findCNEa,findCNEb,seq,seqi]
findCNEa[a0_,b0_,Aseign_,Bseign_,prec_]:=Module[{stepA,stepB},
stepA[{ai_,bi_}]:=If[Abs[Max[Aseign[[;;,bi]]]-Aseign[[ai,bi]]]<10^(-prec),ai,Ordering[Aseign[[;;,bi]],-1][[1]]];
stepB[{ai_,bi_}]:=With[{anext=stepA[{ai,bi}]},
{anext,If[Abs[Max[Bseign[[anext,;;]]]-Bseign[[anext,bi]]]<10^(-prec),bi,Ordering[Bseign[[anext,;;]],-1][[1]]]}];
FixedPoint[stepB,{a0,b0},2prec(*,SameTest->(Abs[#1-#2]<10^(-prec)&)*)]];
findCNEb[a0_,b0_,Aseign_,Bseign_,prec_]:=Module[{stepA,stepB},
stepB[{ai_,bi_}]:=If[Abs[Max[Bseign[[ai,;;]]]-Bseign[[ai,bi]]]<10^(-prec),bi,Ordering[Bseign[[ai,;;]],-1][[1]]];
stepA[{ai_,bi_}]:=With[{bnext=stepB[{ai,bi}]},{If[Abs[Max[Aseign[[;;,bnext]]]-Aseign[[ai,bnext]]]<10^(-prec),ai,Ordering[Aseign[[;;,bnext]],-1][[1]]],bnext}];
FixedPoint[stepA,{a0,b0},2prec(*,SameTest->(Abs[#1-#2]<10^(-prec)&)*)]];seq[δB_,ρ_,gridlength_]:=Module[{u=Log,uinv=Exp,upinv,ℓ=.5,λ=.5,δA=1,(*δB=.98,ρ=0.0,*)M=1,i=.1,prec=13,(*gridlength=21,*)Abar,Bbar,Agrid,Bgrid,gridresults,Aseign,Bseign,asol,bsol},
upinv=z/.Solve[u'[z]==#,z][[1]]&;
Abar=M(1/upinv[1+i/(ℓ-λ ℓ(1-ℓ)δA)]-1);Bbar=Abar;(*M(1/upinv[1+i/(ℓ-λ ℓ(1-ℓ)δB)]-1);*)
Agrid=Range[gridlength]Abar/gridlength;Bgrid=Range[gridlength]Bbar/gridlength;
gridresults=Table[ge[{u,uinv,ℓ,λ,δA,δB,ρ,M,A,B,i,prec}],{A,Agrid},{B,Bgrid}];
Aseign=(If[eN>0,λ ℓ δA (((1+(ℓ/(1-ℓ))(eC/eN)))^(ρ-1)) ((((1-ℓ)eN))^ρ)(u'[q1A]-1),0]/.gridresults)Agrid;
Bseign=Transpose[Transpose[If[eN<1,λ ℓ δB (((1+(ℓ/(1-ℓ))((1-eC)/(1-eN))))^(ρ-1)) ((((1-ℓ)(1-eN)))^ρ)(u'[q1B]-1),0]/.gridresults]Bgrid];
asol=findCNEa[1,1,Aseign,Bseign,prec];
bsol=findCNEb[1,1,Aseign,Bseign,prec];
{{A->Agrid[[asol[[1]]]],B->Bgrid[[asol[[2]]]]},{A->Agrid[[bsol[[1]]]],B->Bgrid[[bsol[[2]]]]}}]
seqi[δB_,ρ_,i_,gridlength_]:=Module[{u=Log,uinv=Exp,upinv,ℓ=.5,λ=.5,δA=1,(*δB=.98,ρ=0.0,*)M=1,prec=13,(*gridlength=21,*)Abar,Bbar,Agrid,Bgrid,gridresults,Aseign,Bseign,asol,bsol},
upinv=z/.Solve[u'[z]==#,z][[1]]&;
Abar=M(1/upinv[1+i/(ℓ-λ ℓ(1-ℓ)δA)]-1);Bbar=Abar;(*M(1/upinv[1+i/(ℓ-λ ℓ(1-ℓ)δB)]-1);*)
Agrid=Range[gridlength]Abar/gridlength;Bgrid=Range[gridlength]Bbar/gridlength;
gridresults=Table[ge[{u,uinv,ℓ,λ,δA,δB,ρ,M,A,B,i,prec}],{A,Agrid},{B,Bgrid}];
Aseign=(If[eN>0,λ ℓ δA (((1+(ℓ/(1-ℓ))(eC/eN)))^(ρ-1)) ((((1-ℓ)eN))^ρ)(u'[q1A]-1),0]/.gridresults)Agrid;
Bseign=Transpose[Transpose[If[eN<1,λ ℓ δB (((1+(ℓ/(1-ℓ))((1-eC)/(1-eN))))^(ρ-1)) ((((1-ℓ)(1-eN)))^ρ)(u'[q1B]-1),0]/.gridresults]Bgrid];
asol=findCNEa[1,1,Aseign,Bseign,prec];
{A->Agrid[[asol[[1]]]],B->Bgrid[[asol[[2]]]]}]

Examples and tests:

seq[1,.0,16]//Timing

Module[{δB=.9,ρ=.0,u=Log,uinv=Exp,upinv,ℓ=.5,λ=.5,δA=1,(*δB=.98,ρ=0.0,*)M=1,i=.1,prec=10,A,B},
A=0.10666666666666665`;
B=0.06399999999999999`;
ge[{u,uinv,ℓ,λ,δA,δB,ρ,M,A,B,i,prec}]]

Module[{δB=.9,ρ=.0,u=Log,uinv=Exp,upinv,ℓ=.5,λ=.5,δA=1,(*δB=.98,ρ=0.0,*)M=1,i=.1,prec=10},
Map[ge[{u,uinv,ℓ,λ,δA,δB,ρ,M,A,B,i,prec}/.#]&,{{A->0.10158730158730156`,B->0.06349206349206347`},{A->0.21587301587301583`,B->0.11428571428571425`}}]]

Computing results on a grid. This takes a long time, therefore the results are saved in tabs.

Basic case. CRS, δA = 1 is held fixed, but δB varies:

Clear[strategicresults]
(* list of "blue" eqa first; list of "green" eqa second *)
strategicresults=Transpose@Table[Map[{δB->aB}~Join~#~Join~ge[{Log,Exp,.5,.5,1,aB,.0,1,A,B,.1,13}/.#]&,seq[aB,.0,25]],{aB,.6,1,.02}];

Print["strategicresults saved here"]
strategicresults={{{δB->0.6`,A->0.12799999999999997`,B->0.010666666666666665`,q0A->0.8096926713947991`,q1A->0.9133333333333332`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.62`,A->0.12799999999999997`,B->0.010666666666666665`,q0A->0.8096926713947991`,q1A->0.9133333333333332`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.64`,A->0.14933333333333332`,B->0.02133333333333333`,q0A->0.8062645011600929`,q1A->0.9266666666666666`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.6599999999999999`,A->0.14933333333333332`,B->0.02133333333333333`,q0A->0.8062645011600929`,q1A->0.9266666666666666`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.6799999999999999`,A->0.15999999999999998`,B->0.02133333333333333`,q0A->0.8045977011494253`,q1A->0.9333333333333333`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.7`,A->0.15999999999999998`,B->0.02133333333333333`,q0A->0.8045311154626384`,q1A->0.9335245139161588`,q0B->0.8155176748089219`,q1B->1.`,eC->0.9979470943206361`,eN->0.9989838245910969`},{δB->0.72`,A->0.17066666666666663`,B->0.031999999999999994`,q0A->0.8029612756264237`,q1A->0.94`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.74`,A->0.17066666666666663`,B->0.031999999999999994`,q0A->0.8026847828586336`,q1A->0.9408310089188683`,q0B->0.8130529458416349`,q1B->1.`,eC->0.991747248006147`,eN->0.9955202958438667`},{δB->0.76`,A->0.1813333333333333`,B->0.04266666666666666`,q0A->0.801354401805869`,q1A->0.9466666666666667`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.78`,A->0.1813333333333333`,B->0.04266666666666666`,q0A->0.8007609405795661`,q1A->0.9485417501732951`,q0B->0.8103068639610654`,q1B->1.`,eC->0.9827695274115884`,eN->0.9897276562387343`},{δB->0.8`,A->0.19199999999999995`,B->0.05333333333333332`,q0A->0.7997762863534675`,q1A->0.9533333333333334`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.82`,A->0.19199999999999995`,B->0.05333333333333332`,q0A->0.7987632514747438`,q1A->0.9567162976899287`,q0B->0.8072498142563707`,q1B->1.`,eC->0.9712342963525282`,eN->0.9811350750822356`},{δB->0.84`,A->0.20266666666666663`,B->0.06399999999999999`,q0A->0.7979578528401103`,q1A->0.960927164513175`,q0B->0.8057808852289524`,q1B->1.`,eC->0.9924046357041049`,eN->0.9948231340922234`},{δB->0.86`,A->0.10666666666666665`,B->0.06399999999999999`,q0A->0.8028004817667391`,q1A->0.9342915472420301`,q0B->0.8080033230369135`,q1B->0.9562772028408018`,eC->0.6527043465285637`,eN->0.7090046725121508`},{δB->0.88`,A->0.10666666666666665`,B->0.06399999999999999`,q0A->0.8027902383105252`,q1A->0.9351272999318125`,q0B->0.8072696423740016`,q1B->0.9536579705417593`,eC->0.6483367809308015`,eN->0.6967062684961717`},{δB->0.8999999999999999`,A->0.10666666666666665`,B->0.06399999999999999`,q0A->0.8027916475870485`,q1A->0.9359549198706633`,q0B->0.8065410244099153`,q1B->0.9511526857942145`,eC->0.6441223245660486`,eN->0.6845091413027875`},{δB->0.9199999999999999`,A->0.09599999999999997`,B->0.07466666666666666`,q0A->0.8027731797650227`,q1A->0.936256779114442`,q0B->0.8057800622637098`,q1B->0.9481303408120164`,eC->0.5782580479484073`,eN->0.6126489332802433`},{δB->0.94`,A->0.09599999999999997`,B->0.07466666666666666`,q0A->0.8028029048605413`,q1A->0.937198125240863`,q0B->0.8050692448704865`,q1B->0.9459952406044482`,eC->0.5741404518061888`,eN->0.5998884966087034`},{δB->0.96`,A->0.09599999999999997`,B->0.07466666666666666`,q0A->0.8028443635004807`,q1A->0.938135291494743`,q0B->0.8043626097499195`,q1B->0.9439322659109488`,eC->0.5701469553487818`,eN->0.5872777265945431`},{δB->0.98`,A->0.08533333333333332`,B->0.08533333333333332`,q0A->0.8028877502118711`,q1A->0.938937352457007`,q0B->0.8036501955641515`,q1B->0.9417981313868604`,eC->0.5038263173191835`,eN->0.5125337566132142`},{δB->1.`,A->0.08533333333333332`,B->0.08533333333333332`,q0A->0.8029612756264237`,q1A->0.94`,q0B->0.8029612756264237`,q1B->0.94`,eC->0.5`,eN->0.5`}},{{δB->0.6`,A->0.13866666666666663`,B->0.10666666666666665`,q0A->0.8079625292740047`,q1A->0.92`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.62`,A->0.13866666666666663`,B->0.10666666666666665`,q0A->0.8079625292740047`,q1A->0.92`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.64`,A->0.14933333333333332`,B->0.10666666666666665`,q0A->0.8062645011600929`,q1A->0.9266666666666666`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.6599999999999999`,A->0.14933333333333332`,B->0.10666666666666665`,q0A->0.8062645011600929`,q1A->0.9266666666666666`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.6799999999999999`,A->0.15999999999999998`,B->0.10666666666666665`,q0A->0.8045977011494253`,q1A->0.9333333333333333`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.7`,A->0.15999999999999998`,B->0.10666666666666665`,q0A->0.804531115462788`,q1A->0.9335245139157292`,q0B->0.8155176748090879`,q1B->1.`,eC->0.9979470943252403`,eN->0.9989838245933738`},{δB->0.72`,A->0.17066666666666663`,B->0.10666666666666665`,q0A->0.8029612756264237`,q1A->0.94`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.74`,A->0.17066666666666663`,B->0.10666666666666665`,q0A->0.8026847828586348`,q1A->0.9408310089188647`,q0B->0.8130529458416358`,q1B->1.`,eC->0.991747248006181`,eN->0.9955202958438855`},{δB->0.76`,A->0.1813333333333333`,B->0.11733333333333332`,q0A->0.801354401805869`,q1A->0.9466666666666667`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.78`,A->0.1813333333333333`,B->0.11733333333333332`,q0A->0.800760940579566`,q1A->0.9485417501732955`,q0B->0.8103068639610668`,q1B->1.`,eC->0.9827695274115853`,eN->0.9897276562387333`},{δB->0.8`,A->0.19199999999999995`,B->0.11733333333333332`,q0A->0.7997762863534675`,q1A->0.9533333333333334`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.82`,A->0.19199999999999995`,B->0.11733333333333332`,q0A->0.7987632514747439`,q1A->0.9567162976899287`,q0B->0.8072498142563703`,q1B->1.`,eC->0.971234296352529`,eN->0.9811350750822321`},{δB->0.84`,A->0.20266666666666663`,B->0.11733333333333332`,q0A->0.797957852840112`,q1A->0.9609271645131694`,q0B->0.8057808852289564`,q1B->1.`,eC->0.9924046357041493`,eN->0.9948231340922549`},{δB->0.86`,A->0.2133333333333333`,B->0.11733333333333332`,q0A->0.7967032967032968`,q1A->0.9666666666666667`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.88`,A->0.2133333333333333`,B->0.11733333333333332`,q0A->0.7958079695646186`,q1A->0.9699723830901489`,q0B->0.8021555904329652`,q1B->1.`,eC->0.9749767410997364`,eN->0.9811394941804594`},{δB->0.8999999999999999`,A->0.22399999999999995`,B->0.11733333333333332`,q0A->0.7948479559789775`,q1A->0.9747139127139269`,q0B->0.8003167402688417`,q1B->1.`,eC->0.9899497216687255`,eN->0.9920738352370031`},{δB->0.9199999999999999`,A->0.23466666666666663`,B->0.11733333333333332`,q0A->0.7937365010799136`,q1A->0.9799999999999999`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.94`,A->0.09599999999999997`,B->0.07466666666666666`,q0A->0.8028029048605413`,q1A->0.937198125240863`,q0B->0.8050692448704865`,q1B->0.9459952406044482`,eC->0.5741404518061888`,eN->0.5998884966087034`},{δB->0.96`,A->0.08533333333333332`,B->0.08533333333333332`,q0A->0.8028262498816657`,q1A->0.9378742309218581`,q0B->0.8043432232586529`,q1B->0.943647514985758`,eC->0.50775701736308`,eN->0.5252639019956303`},{δB->0.98`,A->0.08533333333333332`,B->0.08533333333333332`,q0A->0.8028877502118711`,q1A->0.938937352457007`,q0B->0.8036501955641515`,q1B->0.9417981313868604`,eC->0.5038263173191835`,eN->0.5125337566132142`},{δB->1.`,A->0.08533333333333332`,B->0.08533333333333332`,q0A->0.8029612756264237`,q1A->0.94`,q0B->0.8029612756264237`,q1B->0.94`,eC->0.5`,eN->0.5`}}};

Investigating IRS. δA=δB=1 are held fixed, but ρ varies:

Clear[strategicresultsrho]
(* list of "blue" eqa first; list of "green" eqa second *)
strategicresultsrho=Transpose@Table[Map[{ρ->rho}~Join~#~Join~ge[{Log,Exp,.5,.5,1,1.,rho,1,A,B,.1,13}/.#]&,seq[1.,rho,16]],{rho,0,.06,.006}];

Print["strategicresultsrho saved here"]
strategicresultsrho={{{ρ->0.`,A->0.08333333333333331`,B->0.08333333333333331`,q0A->0.8035714285714285`,q1A->0.9374999999999999`,q0B->0.8035714285714285`,q1B->0.9374999999999999`,eC->0.5`,eN->0.5`},{ρ->0.006`,A->0.08333333333333331`,B->0.08333333333333331`,q0A->0.8036949478245697`,q1A->0.9376441057953313`,q0B->0.8036949478245697`,q1B->0.9376441057953313`,eC->0.5`,eN->0.5`},{ρ->0.012`,A->0.08333333333333331`,B->0.08333333333333331`,q0A->0.8038179544423127`,q1A->0.9377876135160315`,q0B->0.8038179544423127`,q1B->0.9377876135160315`,eC->0.5`,eN->0.5`},{ρ->0.018000000000000002`,A->0.08333333333333331`,B->0.08333333333333331`,q0A->0.8039404505522211`,q1A->0.9379305256442578`,q0B->0.8039404505522211`,q1B->0.9379305256442578`,eC->0.5`,eN->0.5`},{ρ->0.024`,A->0.24999999999999994`,B->0.2666666666666666`,q0A->0.8333333333293677`,q1A->0.8333333333293677`,q0B->0.7894736841605707`,q1B->0.9999999999045405`,eC->0.`,eN->0},{ρ->0.03`,A->0.24999999999999994`,B->0.2666666666666666`,q0A->0.8333333333293677`,q1A->0.8333333333293677`,q0B->0.7894736841605707`,q1B->0.9999999999045405`,eC->0.`,eN->0},{ρ->0.036000000000000004`,A->0.2666666666666666`,B->0.24999999999999994`,q0A->0.7894736841605707`,q1A->0.9999999999045405`,q0B->0.8333333333293677`,q1B->0.8333333333293677`,eC->1.`,eN->1},{ρ->0.042`,A->0.2666666666666666`,B->0.24999999999999994`,q0A->0.7894736841605707`,q1A->0.9999999999045405`,q0B->0.8333333333293677`,q1B->0.8333333333293677`,eC->1.`,eN->1},{ρ->0.048`,A->0.2666666666666666`,B->0.24999999999999994`,q0A->0.7894736841605707`,q1A->0.9999999999045405`,q0B->0.8333333333293677`,q1B->0.8333333333293677`,eC->1.`,eN->1},{ρ->0.054`,A->0.24999999999999994`,B->0.2666666666666666`,q0A->0.8333333333293677`,q1A->0.8333333333293677`,q0B->0.7894736841605707`,q1B->0.9999999999045405`,eC->0.`,eN->0},{ρ->0.06`,A->0.24999999999999994`,B->0.2666666666666666`,q0A->0.8333333333293677`,q1A->0.8333333333293677`,q0B->0.7894736841605707`,q1B->0.9999999999045405`,eC->0.`,eN->0}},{{ρ->0.`,A->0.08333333333333331`,B->0.08333333333333331`,q0A->0.8035714285714285`,q1A->0.9374999999999999`,q0B->0.8035714285714285`,q1B->0.9374999999999999`,eC->0.5`,eN->0.5`},{ρ->0.006`,A->0.08333333333333331`,B->0.08333333333333331`,q0A->0.8036949478245697`,q1A->0.9376441057953313`,q0B->0.8036949478245697`,q1B->0.9376441057953313`,eC->0.5`,eN->0.5`},{ρ->0.012`,A->0.08333333333333331`,B->0.08333333333333331`,q0A->0.8038179544423127`,q1A->0.9377876135160315`,q0B->0.8038179544423127`,q1B->0.9377876135160315`,eC->0.5`,eN->0.5`},{ρ->0.018000000000000002`,A->0.08333333333333331`,B->0.08333333333333331`,q0A->0.8039404505522211`,q1A->0.9379305256442578`,q0B->0.8039404505522211`,q1B->0.9379305256442578`,eC->0.5`,eN->0.5`},{ρ->0.024`,A->0.2666666666666666`,B->0.24999999999999994`,q0A->0.7894736841605707`,q1A->0.9999999999045405`,q0B->0.8333333333293677`,q1B->0.8333333333293677`,eC->1.`,eN->1},{ρ->0.03`,A->0.2666666666666666`,B->0.24999999999999994`,q0A->0.7894736841605707`,q1A->0.9999999999045405`,q0B->0.8333333333293677`,q1B->0.8333333333293677`,eC->1.`,eN->1},{ρ->0.036000000000000004`,A->0.24999999999999994`,B->0.2666666666666666`,q0A->0.8333333333293677`,q1A->0.8333333333293677`,q0B->0.7894736841605707`,q1B->0.9999999999045405`,eC->0.`,eN->0},{ρ->0.042`,A->0.24999999999999994`,B->0.2666666666666666`,q0A->0.8333333333293677`,q1A->0.8333333333293677`,q0B->0.7894736841605707`,q1B->0.9999999999045405`,eC->0.`,eN->0},{ρ->0.048`,A->0.24999999999999994`,B->0.2666666666666666`,q0A->0.8333333333293677`,q1A->0.8333333333293677`,q0B->0.7894736841605707`,q1B->0.9999999999045405`,eC->0.`,eN->0},{ρ->0.054`,A->0.2666666666666666`,B->0.24999999999999994`,q0A->0.7894736841605707`,q1A->0.9999999999045405`,q0B->0.8333333333293677`,q1B->0.8333333333293677`,eC->1.`,eN->1},{ρ->0.06`,A->0.2666666666666666`,B->0.24999999999999994`,q0A->0.7894736841605707`,q1A->0.9999999999045405`,q0B->0.8333333333293677`,q1B->0.8333333333293677`,eC->1.`,eN->1}}};

Investigating B-liquidity, with slight IRS. δA=1 and ρ=0.01 are held fixed, δB varies:

Clear[strategicresultsABirs1]
(* list of "blue" eqa first; list of "green" eqa second *)
strategicresultsABirs1=Transpose@Table[Map[{δB->aB}~Join~#~Join~ge[{Log,Exp,.5,.5,1,aB,.01,1,A,B,.1,13}/.#]&,seq[aB,.01,25]],{aB,.6,1,.02}];

Print["strategicresultsABirs1 saved here"]
strategicresultsABirs1={{{δB->0.6`,A->0.12799999999999997`,B->0.010666666666666665`,q0A->0.8096926713947991`,q1A->0.9133333333333332`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.62`,A->0.12799999999999997`,B->0.010666666666666665`,q0A->0.8096926713947991`,q1A->0.9133333333333332`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.64`,A->0.12799999999999997`,B->0.010666666666666665`,q0A->0.8096926713947991`,q1A->0.9133333333333332`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.6599999999999999`,A->0.14933333333333332`,B->0.02133333333333333`,q0A->0.8062645011600929`,q1A->0.9266666666666666`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.6799999999999999`,A->0.14933333333333332`,B->0.02133333333333333`,q0A->0.8062645011600929`,q1A->0.9266666666666666`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.7`,A->0.14933333333333332`,B->0.02133333333333333`,q0A->0.8054718537510969`,q1A->0.9288894070253082`,q0B->0.8168658769329378`,q1B->1.`,eC->0.9749538045686558`,eN->0.9881670195812693`},{δB->0.72`,A->0.15999999999999998`,B->0.031999999999999994`,q0A->0.8045977011494253`,q1A->0.9333333333333333`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.74`,A->0.15999999999999998`,B->0.031999999999999994`,q0A->0.8038293804682973`,q1A->0.9355814562067774`,q0B->0.8147237773759044`,q1B->1.`,eC->0.9764832342161098`,eN->0.987931066511258`},{δB->0.76`,A->0.17066666666666663`,B->0.031999999999999994`,q0A->0.8029612756264237`,q1A->0.94`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.78`,A->0.17066666666666663`,B->0.031999999999999994`,q0A->0.8019642434933203`,q1A->0.9430608498463949`,q0B->0.8121643818164073`,q1B->1.`,eC->0.9703996205238888`,eN->0.98334414616718`},{δB->0.8`,A->0.1813333333333333`,B->0.031999999999999994`,q0A->0.801354401805869`,q1A->0.9466666666666667`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.82`,A->0.19199999999999995`,B->0.05333333333333332`,q0A->0.7997762863534675`,q1A->0.9533333333333334`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.84`,A->0.19199999999999995`,B->0.06399999999999999`,q0A->0.7997762863534675`,q1A->0.9533333333333334`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.86`,A->0.20266666666666663`,B->0.06399999999999999`,q0A->0.7982261640798226`,q1A->0.96`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.88`,A->0.20266666666666663`,B->0.06399999999999999`,q0A->0.7982261640798226`,q1A->0.96`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.8999999999999999`,A->0.10666666666666665`,B->0.06399999999999999`,q0A->0.802916005513265`,q1A->0.9358396661206152`,q0B->0.8068701304479052`,q1B->0.9520528709051417`,eC->0.6454375279963576`,eN->0.6882470358166441`},{δB->0.9199999999999999`,A->0.09599999999999997`,B->0.07466666666666666`,q0A->0.8029294688871046`,q1A->0.9362668524786727`,q0B->0.8060455887707334`,q1B->0.9486944324961045`,eC->0.5790361055081359`,eN->0.6149139531872286`},{δB->0.94`,A->0.09599999999999997`,B->0.07466666666666666`,q0A->0.8029611186420532`,q1A->0.9372267201653756`,q0B->0.8053336184605697`,q1B->0.9465265874243909`,eC->0.5748389846648543`,eN->0.6019755328352169`},{δB->0.96`,A->0.09599999999999997`,B->0.07466666666666666`,q0A->0.803004694117199`,q1A->0.9381824661491958`,q0B->0.8046260080439466`,q1B->0.944433295955755`,eC->0.570768830177475`,eN->0.5891880262630927`},{δB->0.98`,A->0.08533333333333332`,B->0.08533333333333332`,q0A->0.8030932421177202`,q1A->0.9391589811219953`,q0B->0.8038602674250943`,q1B->0.9420635665693232`,eC->0.5038968309262754`,eN->0.5127210226638682`},{δB->1.`,A->0.08533333333333332`,B->0.08533333333333332`,q0A->0.8031710707536208`,q1A->0.9402456001622388`,q0B->0.8031710707536208`,q1B->0.9402456001622388`,eC->0.5`,eN->0.5`}},{{δB->0.6`,A->0.12799999999999997`,B->0.10666666666666665`,q0A->0.8096926713947991`,q1A->0.9133333333333332`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.62`,A->0.13866666666666663`,B->0.10666666666666665`,q0A->0.8079625292740047`,q1A->0.92`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.64`,A->0.13866666666666663`,B->0.10666666666666665`,q0A->0.8079625292740047`,q1A->0.92`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.6599999999999999`,A->0.14933333333333332`,B->0.10666666666666665`,q0A->0.8062645011600929`,q1A->0.9266666666666666`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.6799999999999999`,A->0.14933333333333332`,B->0.10666666666666665`,q0A->0.8062645011600929`,q1A->0.9266666666666666`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.7`,A->0.14933333333333332`,B->0.10666666666666665`,q0A->0.805471853751093`,q1A->0.9288894070253195`,q0B->0.8168658769329299`,q1B->1.`,eC->0.9749538045685332`,eN->0.9881670195812058`},{δB->0.72`,A->0.15999999999999998`,B->0.10666666666666665`,q0A->0.8045977011494253`,q1A->0.9333333333333333`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.74`,A->0.15999999999999998`,B->0.10666666666666665`,q0A->0.8038293804683009`,q1A->0.9355814562067665`,q0B->0.8147237773759097`,q1B->1.`,eC->0.9764832342162194`,eN->0.987931066511318`},{δB->0.76`,A->0.17066666666666663`,B->0.11733333333333332`,q0A->0.8029612756264237`,q1A->0.94`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.78`,A->0.17066666666666663`,B->0.11733333333333332`,q0A->0.8019642434933154`,q1A->0.9430608498464099`,q0B->0.8121643818164027`,q1B->1.`,eC->0.9703996205237487`,eN->0.983344146167099`},{δB->0.8`,A->0.1813333333333333`,B->0.11733333333333332`,q0A->0.801354401805869`,q1A->0.9466666666666667`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.82`,A->0.19199999999999995`,B->0.11733333333333332`,q0A->0.7997762863534675`,q1A->0.9533333333333334`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.84`,A->0.19199999999999995`,B->0.11733333333333332`,q0A->0.7997762863534675`,q1A->0.9533333333333334`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.86`,A->0.20266666666666663`,B->0.11733333333333332`,q0A->0.7982261640798226`,q1A->0.96`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.88`,A->0.20266666666666663`,B->0.11733333333333332`,q0A->0.7982261640798226`,q1A->0.96`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.8999999999999999`,A->0.2133333333333333`,B->0.11733333333333332`,q0A->0.7967032967032968`,q1A->0.9666666666666667`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.9199999999999999`,A->0.22399999999999995`,B->0.11733333333333332`,q0A->0.7952069716775599`,q1A->0.9733333333333333`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.94`,A->0.22399999999999995`,B->0.11733333333333332`,q0A->0.7952069716775599`,q1A->0.9733333333333333`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.96`,A->0.23466666666666663`,B->0.11733333333333332`,q0A->0.7937365010799136`,q1A->0.9799999999999999`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.98`,A->0.08533333333333332`,B->0.08533333333333332`,q0A->0.8030932421177202`,q1A->0.9391589811219953`,q0B->0.8038602674250943`,q1B->0.9420635665693232`,eC->0.5038968309262754`,eN->0.5127210226638682`},{δB->1.`,A->0.08533333333333332`,B->0.08533333333333332`,q0A->0.8031710707536208`,q1A->0.9402456001622388`,q0B->0.8031710707536208`,q1B->0.9402456001622388`,eC->0.5`,eN->0.5`}}};

Investigating B-liquidity, with slight IRS. δA=1 and ρ=0.01 are held fixed, δB varies. But this time, we double i (the cost of liquidity parameter = Friedman interest rate). Barely any difference:

Clear[strategicresultsABirs12]
(* list of "blue" eqa only *)
strategicresultsABirs12=Table[Apply[{δB->aB}~Join~#~Join~ge[{Log,Exp,.5,.5,1,aB,.01,1,A,B,.2,13}/.#]&,{seqi[aB,.01,.2,25]}],{aB,.6,1,.08}];

(* result with ρ=.01 and i=.2 saved here *)
Print["strategicresultsABirs12 saved here"]
strategicresultsABirs12={{δB->0.6`,A->0.27733333333333327`,B->0.04266666666666666`,q0A->0.6755144646585147`,q1A->0.8628571428571428`,q0B->0.7142857142857143`,q1B->1.`,eC->1.`,eN->1},{δB->0.6799999999999999`,A->0.31999999999999995`,B->0.06399999999999999`,q0A->0.670995670995671`,q1A->0.8857142857142856`,q0B->0.7142857142857143`,q1B->1.`,eC->1.`,eN->1},{δB->0.76`,A->0.34133333333333327`,B->0.08533333333333332`,q0A->0.6673150995607035`,q1A->0.903508252276342`,q0B->0.681776650756367`,q1B->1.`,eC->0.9650964721904248`,eN->0.97768570492803`},{δB->0.84`,A->0.40533333333333327`,B->0.12799999999999997`,q0A->0.6627812415288696`,q1A->0.9314285714285714`,q0B->0.7142857142857143`,q1B->1.`,eC->1.`,eN->1},{δB->0.9199999999999999`,A->0.19199999999999995`,B->0.14933333333333332`,q0A->0.6689009759732822`,q1A->0.8909118397558242`,q0B->0.6732736798678515`,q1B->0.9117970154037859`,eC->0.5800686614584466`,eN->0.6100716812757956`},{δB->1.`,A->0.17066666666666663`,B->0.17066666666666663`,q0A->0.669157966694142`,q1A->0.8975638859924094`,q0B->0.6691579666941423`,q1B->0.897563885992409`,eC->0.4999999999999991`,eN->0.5000000000000011`}};

Investigating B-liquidity, with a bit stronger IRS. δA=1 and ρ=0.02 are held fixed, δB varies:

Clear[strategicresultsABirs2]
(* list of "blue" eqa first; list of "green" eqa second *)
strategicresultsABirs2=Transpose@Table[Map[{δB->aB}~Join~#~Join~ge[{Log,Exp,.5,.5,1,aB,.02,1,A,B,.1,13}/.#]&,seq[aB,.02,25]],{aB,.6,1,.02}];

Print["strategicresultsABirs2 saved here"]
strategicresultsABirs2={{{δB->0.6`,A->0.12799999999999997`,B->0.010666666666666665`,q0A->0.8096926713947991`,q1A->0.9133333333333332`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.62`,A->0.12799999999999997`,B->0.010666666666666665`,q0A->0.8096926713947991`,q1A->0.9133333333333332`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.64`,A->0.12799999999999997`,B->0.010666666666666665`,q0A->0.8096926713947991`,q1A->0.9133333333333332`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.6599999999999999`,A->0.12799999999999997`,B->0.010666666666666665`,q0A->0.8096926713947991`,q1A->0.9133333333333332`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.6799999999999999`,A->0.13866666666666663`,B->0.02133333333333333`,q0A->0.8079625292740047`,q1A->0.92`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.7`,A->0.14933333333333332`,B->0.02133333333333333`,q0A->0.8062645011600929`,q1A->0.9266666666666666`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.72`,A->0.14933333333333332`,B->0.02133333333333333`,q0A->0.8062645011600929`,q1A->0.9266666666666666`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.74`,A->0.15999999999999998`,B->0.031999999999999994`,q0A->0.8045977011494253`,q1A->0.9333333333333333`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.76`,A->0.15999999999999998`,B->0.031999999999999994`,q0A->0.8045977011494253`,q1A->0.9333333333333333`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.78`,A->0.17066666666666663`,B->0.031999999999999994`,q0A->0.8029612756264237`,q1A->0.94`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.8`,A->0.17066666666666663`,B->0.04266666666666666`,q0A->0.8029612756264237`,q1A->0.94`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.82`,A->0.1813333333333333`,B->0.04266666666666666`,q0A->0.801354401805869`,q1A->0.9466666666666667`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.84`,A->0.1813333333333333`,B->0.06399999999999999`,q0A->0.801354401805869`,q1A->0.9466666666666667`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.86`,A->0.19199999999999995`,B->0.06399999999999999`,q0A->0.7997762863534675`,q1A->0.9533333333333334`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.88`,A->0.19199999999999995`,B->0.06399999999999999`,q0A->0.7997762863534675`,q1A->0.9533333333333334`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.8999999999999999`,A->0.20266666666666663`,B->0.06399999999999999`,q0A->0.7982261640798226`,q1A->0.96`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.9199999999999999`,A->0.2133333333333333`,B->0.07466666666666666`,q0A->0.7967032967032968`,q1A->0.9666666666666667`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.94`,A->0.2133333333333333`,B->0.06399999999999999`,q0A->0.7967032967032968`,q1A->0.9666666666666667`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.96`,A->0.22399999999999995`,B->0.08533333333333332`,q0A->0.7952069716775599`,q1A->0.9733333333333333`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.98`,A->0.08533333333333332`,B->0.08533333333333332`,q0A->0.8032971795048217`,q1A->0.9393782214270396`,q0B->0.804069033920302`,q1B->0.9423280711275869`,eC->0.503969517126294`,eN->0.5129145915314444`},{δB->1.`,A->0.2666666666666666`,B->0.25599999999999995`,q0A->0.7894736841605707`,q1A->0.9999999999045405`,q0B->0.8333333333293677`,q1B->1.`,eC->1.`,eN->1}},{{δB->0.6`,A->0.12799999999999997`,B->0.10666666666666665`,q0A->0.8096926713947991`,q1A->0.9133333333333332`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.62`,A->0.12799999999999997`,B->0.10666666666666665`,q0A->0.8096926713947991`,q1A->0.9133333333333332`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.64`,A->0.13866666666666663`,B->0.10666666666666665`,q0A->0.8079625292740047`,q1A->0.92`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.6599999999999999`,A->0.13866666666666663`,B->0.10666666666666665`,q0A->0.8079625292740047`,q1A->0.92`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.6799999999999999`,A->0.13866666666666663`,B->0.10666666666666665`,q0A->0.8079625292740047`,q1A->0.92`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.7`,A->0.14933333333333332`,B->0.10666666666666665`,q0A->0.8062645011600929`,q1A->0.9266666666666666`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.72`,A->0.14933333333333332`,B->0.10666666666666665`,q0A->0.8062645011600929`,q1A->0.9266666666666666`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.74`,A->0.15999999999999998`,B->0.10666666666666665`,q0A->0.8045977011494253`,q1A->0.9333333333333333`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.76`,A->0.15999999999999998`,B->0.11733333333333332`,q0A->0.8045977011494253`,q1A->0.9333333333333333`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.78`,A->0.17066666666666663`,B->0.11733333333333332`,q0A->0.8029612756264237`,q1A->0.94`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.8`,A->0.17066666666666663`,B->0.11733333333333332`,q0A->0.8029612756264237`,q1A->0.94`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.82`,A->0.1813333333333333`,B->0.11733333333333332`,q0A->0.801354401805869`,q1A->0.9466666666666667`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.84`,A->0.1813333333333333`,B->0.11733333333333332`,q0A->0.801354401805869`,q1A->0.9466666666666667`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.86`,A->0.19199999999999995`,B->0.11733333333333332`,q0A->0.7997762863534675`,q1A->0.9533333333333334`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.88`,A->0.19199999999999995`,B->0.11733333333333332`,q0A->0.7997762863534675`,q1A->0.9533333333333334`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.8999999999999999`,A->0.20266666666666663`,B->0.11733333333333332`,q0A->0.7982261640798226`,q1A->0.96`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.9199999999999999`,A->0.2133333333333333`,B->0.11733333333333332`,q0A->0.7967032967032968`,q1A->0.9666666666666667`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.94`,A->0.2133333333333333`,B->0.11733333333333332`,q0A->0.7967032967032968`,q1A->0.9666666666666667`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.96`,A->0.22399999999999995`,B->0.11733333333333332`,q0A->0.7952069716775599`,q1A->0.9733333333333333`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.98`,A->0.23466666666666663`,B->0.11733333333333332`,q0A->0.7937365010799136`,q1A->0.9799999999999999`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->1.`,A->0.25599999999999995`,B->0.2666666666666666`,q0A->0.8333333333293677`,q1A->1.`,q0B->0.7894736841605707`,q1B->0.9999999999045405`,eC->0.`,eN->0}}};

Figures for the paper

With ρ=0:

Module[{ρ=.0,u=Log,uinv=Exp,upinv,ℓ=.5,λ=.5,δA=1,(*δB=.98,ρ=0.0,*)M=1,i=.1,prec=10,Abar},
upinv=z/.Solve[u'[z]==#,z][[1]]&;
Abar=M(1/upinv[1+i/(ℓ-λ ℓ(1-ℓ)δA)]-1);
ListPlot[Transpose[{{δB,A 100/Abar},{δB,If[eC<1,B ,Null]100/Abar}}/.strategicresults[[1]]],PlotRange->{{.6,1},{0,98}},Joined->True,AxesLabel->{" (Subscript[δ, B])","(% of (OverBar[A]))"},BaseStyle->{FontSize->18},PlotStyle->{{Thick,Darker@Blue},{Thick,Dashed,Darker@Red}},
Epilog->{Text[Style["A",Darker@Blue,FontSize->20],{.8,.21 100/Abar}],Text[Style["B",Darker@Red,FontSize->20],{.8,.06 100/Abar}]
}]]
Export["csab-crs-assetsupply.pdf",%]

Module[{ρ=.0,u=Log,uinv=Exp,upinv,ℓ=.5,λ=.5,δA=1,(*δB=.98,ρ=0.0,*)M=1,i=.1,prec=10,Abar},
upinv=z/.Solve[u'[z]==#,z][[1]]&;
Abar=M(1/upinv[1+i/(ℓ-λ ℓ(1-ℓ)δA)]-1);
ListPlot[Transpose[{{δB,100If[eN>0,λ ℓ δA (((1+(ℓ/(1-ℓ))(eC/eN)))^(ρ-1)) ((((1-ℓ)eN))^ρ)(u'[q1A]-1),Null]},{δB,100If[eN<1,λ ℓ δB (((1+(ℓ/(1-ℓ))((1-eC)/(1-eN))))^(ρ-1)) ((((1-ℓ)(1-eN)))^ρ)(u'[q1B]-1),0]}}/.strategicresults[[1]]],PlotRange->{{.6,1},{0,100λ ℓ(1-ℓ)δA i}},Joined->True,AxesLabel->{" (Subscript[δ, B])","(%)"},BaseStyle->{FontSize->18},PlotStyle->{{Thick,Darker@Blue},{Thick,Dashed,Darker@Red}},
Epilog->{Text[Style["(Subscript[L, A])",Darker@Blue,FontSize->20],{.9,1}],Text[Style["(Subscript[L, B])",Darker@Red,FontSize->20],{.9,.3}]
}]]
Export["csab-crs-premia.pdf",%]

Module[{ρ=.0,u=Log,uinv=Exp,upinv,ℓ=.5,λ=.5,δA=1,(*δB=.98,ρ=0.0,*)M=1,i=.1,prec=10,Abar},
upinv=z/.Solve[u'[z]==#,z][[1]]&;
Abar=M(1/upinv[1+i/(ℓ-λ ℓ(1-ℓ)δA)]-1);
ListPlot[Transpose[{{δB,If[eN>0, δA (((1+(ℓ/(1-ℓ))(eC/eN)))^(ρ-1)) ((((1-ℓ)eN))^ρ),0]},{δB,If[eN<1,δB (((1+(ℓ/(1-ℓ))((1-eC)/(1-eN))))^(ρ-1)) ((((1-ℓ)(1-eN)))^ρ),0]},{δB,(1-ℓ)δB},{δB,(1-ℓ)δA}}/.strategicresults[[1]]],PlotRange->{{.6,1},{.31,.59}},Joined->True,AxesLabel->{" (Subscript[δ, B])",None},BaseStyle->{FontSize->18},PlotStyle->{{Thick,Darker@Blue},{Thick,Dashed,Darker@Red},{Thin,Dashed,Gray},{Thin,Dashed,Gray}},
Epilog->{Text[Style["(Subscript[α, CA])",Darker@Blue,FontSize->20],{.8,.53}],Text[Style["(Subscript[α, CB])",Darker@Red,FontSize->20],{.9,.38}],Text[Style["(1-ℓ)(Subscript[δ, A])",Darker@Gray,FontSize->16],{.86,.48}],Text[Style["(1-ℓ)(Subscript[δ, B])",Darker@Gray,FontSize->16],{.69,.38}]
}]]
Export["csab-crs-arrival.pdf",%]

Module[{ρ=.0,u=Log,uinv=Exp,upinv,ℓ=.5,λ=.5,δA=1,(*δB=.98,ρ=0.0,*)M=1,i=.1,prec=10,Abar,fA,fB,q1Atil,q1Btil},
upinv=z/.Solve[u'[z]==#,z][[1]]&;
Abar=M(1/upinv[1+i/(ℓ-λ ℓ(1-ℓ)δA)]-1);
fA=If[eC>0,δA((((1/(ℓ eC))+(1/((1-ℓ)eN))))^(-1)) ,0];
fB=If[eC<1,δB((((1/(ℓ(1-eC)))+(1/((1-ℓ)(1-eN)))))^(-1)),0];
q1Atil=If[eC>0,Min[1,uinv[u[q0A]+(A/(eC M))(eC q0A+(1-eC)q0B)]],1];
q1Btil=If[eC<1,Min[1,uinv[u[q0B]+(B/((1-eC) M))(eC q0A+(1-eC)q0B)]],1];
ListPlot[Transpose[{{δB,100q0A},{δB,100q1A},{δB,100If[eC<1,q0B,Null]},{δB,100If[eC<1,q1B,Null]},{δB,100(fA(λ q1A+(1-λ)q1Atil)+fB(λ q1B+(1-λ)q1Btil)+(ℓ eC-fA)q0A+(ℓ(1-eC)-fB)q0B)/ℓ}}/.strategicresults[[1]]],
PlotRange->{{.6,1},{77,101}},Joined->True,AxesLabel->{" (Subscript[δ, B])","(% of (SuperStar[q]))"},BaseStyle->{FontSize->18},PlotStyle->{{Thick,Dashed,Darker@Blue},{Thick,Darker@Blue},{Thick,Dashed,Darker@Red},{Thick,Darker@Red},{Thick,DotDashed,RGBColor[.33,0,.33]}},
Epilog->{Text[Style["(Subscript[q, 0  A])",Darker@Blue,FontSize->18],{.7,82}],Text[Style["(Subscript[q, 1  A])",Darker@Blue,FontSize->18],{.7,95.5}],Text[Style["(Subscript[q, 0  B])",Darker@Red,FontSize->18],{.9,82.5}],Text[Style["(Subscript[q, 1  B])",Darker@Red,FontSize->18],{.9,97}],Text[Style["(OverBar[q])",RGBColor[.33,0,.33],FontSize->18],{.8,89}]
}]]
Export["csab-crs-allqs.pdf",%]

With ρ=.01:

Module[{ρ=.01,u=Log,uinv=Exp,upinv,ℓ=.5,λ=.5,δA=1,(*δB=.98,ρ=0.0,*)M=1,i=.1,prec=10,Abar},
upinv=z/.Solve[u'[z]==#,z][[1]]&;
Abar=M(1/upinv[1+i/(ℓ-λ ℓ(1-ℓ)δA)]-1);
ListPlot[Transpose[{{δB,A 100/Abar},{δB,If[eC<1,B,Null]100/Abar}}/.strategicresultsABirs1[[1]]],PlotRange->{{.6,1},{0,98}},Joined->True,AxesLabel->{" (Subscript[δ, B])","(% of (OverBar[A]))"},BaseStyle->{FontSize->18},PlotStyle->{{Thick,Darker@Blue},{Thick,Dashed,Darker@Red}},
Epilog->{Text[Style["A",Darker@Blue,FontSize->20],{.9,.21 100/Abar}],Text[Style["B",Darker@Red,FontSize->20],{.9,.04 100/Abar}]
}]]
Export["csab-irs1-assetsupply.pdf",%]

Module[{ρ=.01,u=Log,uinv=Exp,upinv,ℓ=.5,λ=.5,δA=1,(*δB=.98,ρ=0.0,*)M=1,i=.1,prec=10,Abar},
upinv=z/.Solve[u'[z]==#,z][[1]]&;
Abar=M(1/upinv[1+i/(ℓ-λ ℓ(1-ℓ)δA)]-1);
ListPlot[Transpose[{{δB,100If[eN>0,λ ℓ δA (((1+(ℓ/(1-ℓ))(eC/eN)))^(ρ-1)) ((((1-ℓ)eN))^ρ)(u'[q1A]-1),0]},{δB,100If[eN<1,λ ℓ δB (((1+(ℓ/(1-ℓ))((1-eC)/(1-eN))))^(ρ-1)) ((((1-ℓ)(1-eN)))^ρ)(u'[q1B]-1),0]}}/.strategicresultsABirs1[[1]]],PlotRange->{{.6,1},{0,100λ ℓ(1-ℓ)δA i}},Joined->True,AxesLabel->{" (Subscript[δ, B])","(%)"},BaseStyle->{FontSize->18},PlotStyle->{{Thick,Darker@Blue},{Thick,Dashed,Darker@Red}},
Epilog->{Text[Style["(Subscript[L, A])",Darker@Blue,FontSize->20],{.92,1}],Text[Style["(Subscript[L, B])",Darker@Red,FontSize->20],{.92,.3}]
}]]
Export["csab-irs1-premia.pdf",%]

Module[{ρ=.01,u=Log,uinv=Exp,upinv,ℓ=.5,λ=.5,δA=1,(*δB=.98,ρ=0.0,*)M=1,i=.1,prec=10,Abar},
upinv=z/.Solve[u'[z]==#,z][[1]]&;
Abar=M(1/upinv[1+i/(ℓ-λ ℓ(1-ℓ)δA)]-1);
ListPlot[Transpose[{{δB,If[eN>0, δA (((1+(ℓ/(1-ℓ))(eC/eN)))^(ρ-1)) ((((1-ℓ)eN))^ρ),0]},{δB,If[eN<.99,δB (((1+(ℓ/(1-ℓ))((1-eC)/(1-eN))))^(ρ-1)) ((((1-ℓ)(1-eN)))^ρ),0]},{δB,(1-ℓ)δB},{δB,(1-ℓ)δA}}/.strategicresultsABirs1[[1]]],PlotRange->{{.6,1},{.31,.59}},Joined->True,AxesLabel->{" (Subscript[δ, B])",None},BaseStyle->{FontSize->18},PlotStyle->{{Thick,Darker@Blue},{Thick,Dashed,Darker@Red},{Thin,Dashed,Gray},{Thin,Dashed,Gray}},
Epilog->{Text[Style["(Subscript[α, CA])",Darker@Blue,FontSize->20],{.8,.53}],Text[Style["(Subscript[α, CB])",Darker@Red,FontSize->20],{.93,.38}],Text[Style["(1-ℓ)(Subscript[δ, A])",Darker@Gray,FontSize->16],{.88,.48}],Text[Style["(1-ℓ)(Subscript[δ, B])",Darker@Gray,FontSize->16],{.82,.38}]
}]]
Export["csab-irs1-arrival.pdf",%]

Module[{ρ=.01,u=Log,uinv=Exp,upinv,ℓ=.5,λ=.5,δA=1,(*δB=.98,ρ=0.0,*)M=1,i=.1,prec=10,Abar,fA,fB,q1Atil,q1Btil},
upinv=z/.Solve[u'[z]==#,z][[1]]&;
Abar=M(1/upinv[1+i/(ℓ-λ ℓ(1-ℓ)δA)]-1);
fA=If[eC>0,δA((((1/(ℓ eC))+(1/((1-ℓ)eN))))^(ρ-1)) (((ℓ eC(1-ℓ)eN))^ρ),0];
fB=If[eC<1,δB((((1/(ℓ(1-eC)))+(1/((1-ℓ)(1-eN)))))^(ρ-1))(((ℓ (1-eC)(1-ℓ)(1-eN)))^ρ),0];
q1Atil=If[eC>0,Min[1,uinv[u[q0A]+(A/(eC M))(eC q0A+(1-eC)q0B)]],1];
q1Btil=If[eC<1,Min[1,uinv[u[q0B]+(B/((1-eC) M))(eC q0A+(1-eC)q0B)]],1];
ListPlot[Transpose[{{δB,100q0A},{δB,100q1A},{δB,100If[eC<1,q0B,Null]},{δB,100If[eC<1,q1B,Null]},{δB,100(fA(λ q1A+(1-λ)q1Atil)+fB(λ q1B+(1-λ)q1Btil)+(ℓ eC-fA)q0A+(ℓ(1-eC)-fB)q0B)/ℓ}}/.strategicresultsABirs1[[1]]],
PlotRange->{{.6,1},{77,101}},Joined->True,AxesLabel->{" (Subscript[δ, B])","(% of (SuperStar[q]))"},BaseStyle->{FontSize->18},PlotStyle->{{Thick,Dashed,Darker@Blue},{Thick,Darker@Blue},{Thick,Dashed,Darker@Red},{Thick,Darker@Red},{Thick,DotDashed,RGBColor[.33,0,.33]}},
Epilog->{Text[Style["(Subscript[q, 0  A])",Darker@Blue,FontSize->18],{.7,82.5}],Text[Style["(Subscript[q, 1  A])",Darker@Blue,FontSize->18],{.7,95.5}],Text[Style["(Subscript[q, 0  B])",Darker@Red,FontSize->18],{.95,82.5}],Text[Style["(Subscript[q, 1  B])",Darker@Red,FontSize->18],{.95,97}],Text[Style["(OverBar[q])",RGBColor[.33,0,.33],FontSize->18],{.8,89}]
}]]
Export["csab-irs1-allqs.pdf",%]

With ρ=.02:

Module[{ρ=.02,u=Log,uinv=Exp,upinv,ℓ=.5,λ=.5,δA=1,(*δB=.98,ρ=0.0,*)M=1,i=.1,prec=10,Abar},
upinv=z/.Solve[u'[z]==#,z][[1]]&;
Abar=M(1/upinv[1+i/(ℓ-λ ℓ(1-ℓ)δA)]-1);
ListPlot[Transpose[{{δB,A 100/Abar},{δB,If[eC<1,B,Null]100/Abar}}/.strategicresultsABirs2[[2]]],PlotRange->{{.6,1},{0,98}},Joined->True,AxesLabel->{" (Subscript[δ, B])","(% of (OverBar[A]))"},BaseStyle->{FontSize->18},PlotStyle->{{Thick,Darker@Blue},{Thick,Dashed,Darker@Red}},
Epilog->{Text[Style["A",Darker@Blue,FontSize->20],{.9,.23 100/Abar}],Text[Style["(B undetermined)",Darker@Red,FontSize->16],{.9,.08 100/Abar}]
}]]
Export["csab-irs2-assetsupply.pdf",%]

Module[{ρ=.02,u=Log,uinv=Exp,upinv,ℓ=.5,λ=.5,δA=1,(*δB=.98,ρ=0.0,*)M=1,i=.1,prec=10,Abar},
upinv=z/.Solve[u'[z]==#,z][[1]]&;
Abar=M(1/upinv[1+i/(ℓ-λ ℓ(1-ℓ)δA)]-1);
ListPlot[Transpose[{{δB,100If[eN>0,λ ℓ δA (((1+(ℓ/(1-ℓ))(eC/eN)))^(ρ-1)) ((((1-ℓ)eN))^ρ)(u'[q1A]-1),0]},{δB,100If[eN<1,λ ℓ δB (((1+(ℓ/(1-ℓ))((1-eC)/(1-eN))))^(ρ-1)) ((((1-ℓ)(1-eN)))^ρ)(u'[q1B]-1),0]+.005}}/.strategicresultsABirs2[[2]]],PlotRange->{{.6,1},{0,100λ ℓ(1-ℓ)δA i}},Joined->True,AxesLabel->{" (Subscript[δ, B])","(%)"},BaseStyle->{FontSize->18},PlotStyle->{{Thick,Darker@Blue},{Thick,Dashed,Darker@Red}},
Epilog->{Text[Style["(Subscript[L, A])",Darker@Blue,FontSize->20],{.9,.8}],Text[Style["(Subscript[L, B])",Darker@Red,FontSize->20],{.9,.1}]
}]]
Export["csab-irs2-premia.pdf",%]

Module[{ρ=.02,u=Log,uinv=Exp,upinv,ℓ=.5,λ=.5,δA=1,(*δB=.98,ρ=0.0,*)M=1,i=.1,prec=10,Abar},
upinv=z/.Solve[u'[z]==#,z][[1]]&;
Abar=M(1/upinv[1+i/(ℓ-λ ℓ(1-ℓ)δA)]-1);
ListPlot[Transpose[{{δB,If[eN>0, δA (((1+(ℓ/(1-ℓ))(eC/eN)))^(ρ-1)) ((((1-ℓ)eN))^ρ),Null]},{δB,If[eN<1,δB (((1+(ℓ/(1-ℓ))((1-eC)/(1-eN))))^(ρ-1)) ((((1-ℓ)(1-eN)))^ρ),0]},{δB,(1-ℓ)δB},{δB,(1-ℓ)δA}}/.strategicresultsABirs2[[2]]],PlotRange->{{.6,1},{.31,.59}},Joined->True,AxesLabel->{" (Subscript[δ, B])",None},BaseStyle->{FontSize->18},PlotStyle->{{Thick,Darker@Blue},{Thick,Dashed,Darker@Red},{Thin,Dashed,Gray},{Thin,Dashed,Gray}},
Epilog->{Text[Style["(Subscript[α, CA])",Darker@Blue,FontSize->20],{.8,.53}],Text[Style["(Subscript[α, CB])",Darker@Red,FontSize->20],{.97,.38}],Text[Style["= (1-ℓ)(Subscript[δ, A])",Darker@Gray,FontSize->16],{.88,.528}],Text[Style["(1-ℓ)(Subscript[δ, B])",Darker@Gray,FontSize->16],{.83,.38}]
}]]
Export["csab-irs2-arrival.pdf",%]

Module[{ρ=.02,u=Log,uinv=Exp,upinv,ℓ=.5,λ=.5,δA=1,(*δB=.98,ρ=0.0,*)M=1,i=.1,prec=10,Abar,fA,fB,q1Atil,q1Btil},
upinv=z/.Solve[u'[z]==#,z][[1]]&;
Abar=M(1/upinv[1+i/(ℓ-λ ℓ(1-ℓ)δA)]-1);
fA=If[eC>0,δA((((1/(ℓ eC))+(1/((1-ℓ)eN))))^(ρ-1)) (((ℓ eC(1-ℓ)eN))^ρ),0];
fB=If[eC<1,δB((((1/(ℓ(1-eC)))+(1/((1-ℓ)(1-eN)))))^(ρ-1))(((ℓ (1-eC)(1-ℓ)(1-eN)))^ρ),0];
q1Atil=If[eC>0,Min[1,uinv[u[q0A]+(A/(eC M))(eC q0A+(1-eC)q0B)]],1];
q1Btil=If[eC<1,Min[1,uinv[u[q0B]+(B/((1-eC) M))(eC q0A+(1-eC)q0B)]],1];
ListPlot[Transpose[{{δB,100If[eC>0,q0A,Null]},{δB,100q1A},{δB,100If[eC<1,q0B,Null]},{δB,100If[eC<1,q1B,Null]},{δB,100(fA(λ q1A+(1-λ)q1Atil)+fB(λ q1B+(1-λ)q1Btil)+(ℓ eC-fA)q0A+(ℓ(1-eC)-fB)q0B)/ℓ}}/.strategicresultsABirs2[[2]]],
PlotRange->{{.6,1},{77,101}},Joined->True,AxesLabel->{" (Subscript[δ, B])","(% of (SuperStar[q]))"},BaseStyle->{FontSize->18},PlotStyle->{{Thick,Dashed,Darker@Blue},{Thick,Darker@Blue},{Thick,Dashed,Darker@Red},{Thick,Darker@Red},{Thick,DotDashed,RGBColor[.33,0,.33]}},
Epilog->{Text[Style["(Subscript[q, 0  A])",Darker@Blue,FontSize->18],{.9,82}],Text[Style["(Subscript[q, 1  A])",Darker@Blue,FontSize->18],{.9,98}],Text[Style["((Subscript[q, 0  B]), (Subscript[q, 1  B])\nundetermined)",Darker@Red,FontSize->16],{.71,98}],Text[Style["(OverBar[q])",RGBColor[.33,0,.33],FontSize->18],{.9,89}]
}]]
Export["csab-irs2-allqs.pdf",%]

Welfare:

Module[{ρ=.0,u=Log,uinv=Exp,upinv,v,ℓ=.5,λ=.5,δA=1,(*δB=.98,ρ=0.0,*)M=1,i=.1,prec=10,Abar,fA,fB,q1Atil,q1Btil},
upinv=z/.Solve[u'[z]==#,z][[1]]&;
v=(u[#]-#+1)&;
Abar=M(1/upinv[1+i/(ℓ-λ ℓ(1-ℓ)δA)]-1);
fA=If[eC>0,δA((((1/(ℓ eC))+(1/((1-ℓ)eN))))^(-1)) ,0];
fB=If[eC<1,δB((((1/(ℓ(1-eC)))+(1/((1-ℓ)(1-eN)))))^(-1)),0];
q1Atil=If[eC>0,Min[1,uinv[u[q0A]+(A/(eC M))(eC q0A+(1-eC)q0B)]],1];
q1Btil=If[eC<1,Min[1,uinv[u[q0B]+(B/((1-eC) M))(eC q0A+(1-eC)q0B)]],1];
ListPlot[{{δB,100(fA(λ v[q1A]+(1-λ)v[q1Atil])+fB(λ v[q1B]+(1-λ)v[q1Btil])+(ℓ eC-fA)v[q0A]+(ℓ(1-eC)-fB)v[q0B])/ℓ}/.strategicresults[[1]]},
PlotRange->{{.6,1},{-1.3,-1.235}},AxesOrigin->{.6,-1.3},AspectRatio->1,ImageSize->260,Joined->True,AxesLabel->{" (Subscript[δ, B])","(%)"},BaseStyle->{FontSize->18},PlotStyle->{{Thick,Darker@Green}}]]
Export["csab-crs-welfare.pdf",%]

Module[{ρ=.01,u=Log,uinv=Exp,upinv,v,ℓ=.5,λ=.5,δA=1,(*δB=.98,ρ=0.0,*)M=1,i=.1,prec=10,Abar,fA,fB,q1Atil,q1Btil},
upinv=z/.Solve[u'[z]==#,z][[1]]&;
v=(u[#]-#+1)&;
Abar=M(1/upinv[1+i/(ℓ-λ ℓ(1-ℓ)δA)]-1);
fA=If[eC>0,δA((((1/(ℓ eC))+(1/((1-ℓ)eN))))^(ρ-1)) (((ℓ eC(1-ℓ)eN))^ρ),0];
fB=If[eC<1,δB((((1/(ℓ(1-eC)))+(1/((1-ℓ)(1-eN)))))^(ρ-1))(((ℓ (1-eC)(1-ℓ)(1-eN)))^ρ),0];
q1Atil=If[eC>0,Min[1,uinv[u[q0A]+(A/(eC M))(eC q0A+(1-eC)q0B)]],1];
q1Btil=If[eC<1,Min[1,uinv[u[q0B]+(B/((1-eC) M))(eC q0A+(1-eC)q0B)]],1];
ListPlot[{{δB,100(fA(λ v[q1A]+(1-λ)v[q1Atil])+fB(λ v[q1B]+(1-λ)v[q1Btil])+(ℓ eC-fA)v[q0A]+(ℓ(1-eC)-fB)v[q0B])/ℓ}/.strategicresultsABirs1[[1]]},
PlotRange->{{.6,1},{-1.3,-1.235}},AxesOrigin->{.6,-1.3},AspectRatio->1,ImageSize->260,Joined->True,AxesLabel->{" (Subscript[δ, B])","(%)"},BaseStyle->{FontSize->18},PlotStyle->{{Thick,Darker@Green}}]]
Export["csab-irs1-welfare.pdf",%]

Module[{ρ=.02,u=Log,uinv=Exp,upinv,v,ℓ=.5,λ=.5,δA=1,(*δB=.98,ρ=0.0,*)M=1,i=.1,prec=13,Abar,fA,fB,q1Atil,q1Btil},
upinv=z/.Solve[u'[z]==#,z][[1]]&;
v=(u[#]-#+1)&;
Abar=M(1/upinv[1+i/(ℓ-λ ℓ(1-ℓ)δA)]-1);
fA=If[eC>0,δA((((1/(ℓ eC))+(1/((1-ℓ)eN))))^(ρ-1)) (((ℓ eC(1-ℓ)eN))^ρ),0];
fB=If[eC<1,δB((((1/(ℓ(1-eC)))+(1/((1-ℓ)(1-eN)))))^(ρ-1))(((ℓ (1-eC)(1-ℓ)(1-eN)))^ρ),0];
q1Atil=If[eC>0,Min[1,uinv[u[q0A]+(A/(eC M))(eC q0A+(1-eC)q0B)]],1];
q1Btil=If[eC<1,Min[1,uinv[u[q0B]+(B/((1-eC) M))(eC q0A+(1-eC)q0B)]],1];
ListPlot[{{δB,100(fA(λ v[q1A]+(1-λ)v[q1Atil])+fB(λ v[q1B]+(1-λ)v[q1Btil])+(ℓ eC-fA)v[q0A]+(ℓ(1-eC)-fB)v[q0B])/ℓ}/.strategicresultsABirs2[[2]]},
PlotRange->{{.6,1},{-1.3,-1.235}},AxesOrigin->{.6,-1.3},AspectRatio->1,ImageSize->260,Joined->True,AxesLabel->{" (Subscript[δ, B])","(%)"},BaseStyle->{FontSize->18},PlotStyle->{{Thick,Darker@Green}}(*,Epilog->{Black,Text[Style["ρ = "<>ToString[ρ]],{.8,-1.237}]}*)]]
Export["csab-irs2-welfare.pdf",%]

Analyzing semi-strategic equilibria where A is fixed

Algorithm to find semi-strategic equilibria

Clear[findBR,sseq,sseqi]
findBR[a0_,b0_,Bseign_,prec_]:=Module[{stepB},
stepB[bi_]:=If[Abs[Max[Bseign[[;;]]]-Bseign[[bi]]]<10^(-prec),bi,Ordering[Bseign[[;;]],-1][[1]]];
FixedPoint[stepB,b0,2prec(*,SameTest->(Abs[#1-#2]<10^(-prec)&)*)]];sseq[δB_,ρ_,afrac_,gridlength_]:=Module[{u=Log,uinv=Exp,upinv,ℓ=.5,λ=.5,δA=1,(*δB=.98,ρ=0.0,*)M=1,i=.1,prec=13,(*gridlength=21,*)a0,Abar,Bbar,Bgrid,gridresults,Bseign,bsol},
upinv=z/.Quiet@Solve[u'[z]==#,z][[1]]&;
Abar=M(1/upinv[1+i/(ℓ-λ ℓ(1-ℓ)δA)]-1);Bbar=Abar;(*M(1/upinv[1+i/(ℓ-λ ℓ(1-ℓ)δB)]-1);*)
a0=afrac Abar;
Bgrid=Range[gridlength]Bbar/gridlength;
gridresults=Table[ge[{u,uinv,ℓ,λ,δA,δB,ρ,M,a0,B,i,prec}],{B,Bgrid}];
Bseign=(If[eN<1,λ ℓ δB (((1+(ℓ/(1-ℓ))((1-eC)/(1-eN))))^(ρ-1)) ((((1-ℓ)(1-eN)))^ρ)(u'[q1B]-1),0]/.gridresults )Bgrid;
bsol=findBR[1,1,Bseign,prec];
{{A->a0,B->Bgrid[[bsol]]}}]
sseqi[δB_,ρ_,i_,afrac_,gridlength_]:=Module[{u=Log,uinv=Exp,upinv,ℓ=.5,λ=.5,δA=1,(*δB=.98,ρ=0.0,*)M=1,prec=13,(*gridlength=21,*)a0,Abar,Bbar,Bgrid,gridresults,Bseign,bsol},
upinv=z/.Quiet@Solve[u'[z]==#,z][[1]]&;
Abar=M(1/upinv[1+i/(ℓ-λ ℓ(1-ℓ)δA)]-1);Bbar=Abar;(*M(1/upinv[1+i/(ℓ-λ ℓ(1-ℓ)δB)]-1);*)
a0=afrac Abar;
Bgrid=Range[gridlength]Bbar/gridlength;
gridresults=Table[ge[{u,uinv,ℓ,λ,δA,δB,ρ,M,a0,B,i,prec}],{B,Bgrid}];
Bseign=Transpose[Transpose[If[eN<1,λ ℓ δB (((1+(ℓ/(1-ℓ))((1-eC)/(1-eN))))^(ρ-1)) ((((1-ℓ)(1-eN)))^ρ)(u'[q1B]-1),0]/.gridresults]Bgrid];
bsol=findBR[1,1,Bseign,prec];
{{A->a0,B->Bgrid[[bsol]]}}]

Examples and tests:

sseq[1,.0,.5,36]//Timing

Module[{δB=.9,ρ=.0,u=Log,uinv=Exp,upinv,ℓ=.5,λ=.5,δA=1,(*δB=.98,ρ=0.0,*)M=1,i=.1,prec=10,A,B},
A=0.10666666666666665`;
B=0.06399999999999999`;
ge[{u,uinv,ℓ,λ,δA,δB,ρ,M,A,B,i,prec}]]

Module[{δB=.9,ρ=.0,u=Log,uinv=Exp,upinv,ℓ=.5,λ=.5,δA=1,(*δB=.98,ρ=0.0,*)M=1,i=.1,prec=10},
Map[ge[{u,uinv,ℓ,λ,δA,δB,ρ,M,A,B,i,prec}/.#]&,{{A->0.10158730158730156`,B->0.06349206349206347`},{A->0.21587301587301583`,B->0.11428571428571425`}}]]

Results take long to compute, so they are saved here:

Clear[semistrategicresultsA5]
semistrategicresultsA5=Table[Map[{δB->aB}~Join~#~Join~ge[{Log,Exp,.5,.5,1,aB,.0,1,A,B,.1,13}/.#]&,sseq[aB,.0,.5,36]][[1]],{aB,.6,1,.02}];

Print["semistrategicresults with A=Abar/2 saved here"]
semistrategicresultsA5={{δB->0.6`,A->0.1333333333333333`,B->0.007407407407407405`,q0A->0.8088235294117646`,q1A->0.9166666666666666`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.62`,A->0.1333333333333333`,B->0.007407407407407405`,q0A->0.8088235294117646`,q1A->0.9166666666666666`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.64`,A->0.1333333333333333`,B->0.007407407407407405`,q0A->0.8088235294117646`,q1A->0.9166666666666666`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.6599999999999999`,A->0.1333333333333333`,B->0.01481481481481481`,q0A->0.806672320490298`,q1A->0.9228689372500605`,q0B->0.8165656704055115`,q1B->0.9792522153754676`,eC->0.9264753122408166`,eN->0.9573693522642425`},{δB->0.6799999999999999`,A->0.1333333333333333`,B->0.01481481481481481`,q0A->0.806639103465362`,q1A->0.9233537226900219`,q0B->0.8152734237321552`,q1B->0.9691241198001265`,eC->0.9222613010008408`,eN->0.9478470119935888`},{δB->0.7`,A->0.1333333333333333`,B->0.01481481481481481`,q0A->0.8066144326127791`,q1A->0.923577445593528`,q0B->0.8146179798484687`,q1B->0.9645536985643938`,eC->0.9202371448536059`,eN->0.9439285017926485`},{δB->0.72`,A->0.1333333333333333`,B->0.022222222222222216`,q0A->0.8054668497683719`,q1A->0.9263618808084426`,q0B->0.8148882512306482`,q1B->0.9770610024845189`,eC->0.8894857431399713`,eN->0.9349757068963337`},{δB->0.74`,A->0.1333333333333333`,B->0.02962962962962962`,q0A->0.8042996229432327`,q1A->0.9299615637484533`,q0B->0.813514790238202`,q1B->0.9779363906873023`,eC->0.8548199255397886`,eN->0.9101472066074959`},{δB->0.76`,A->0.1333333333333333`,B->0.02962962962962962`,q0A->0.8042455721184749`,q1A->0.930421367120776`,q0B->0.8127448532451196`,q1B->0.9731456598770579`,eC->0.8512041786565611`,eN->0.9030752259629775`},{δB->0.78`,A->0.1333333333333333`,B->0.03703703703703703`,q0A->0.8030626851719108`,q1A->0.9343096893261463`,q0B->0.8112755117658517`,q1B->0.9744235234247579`,eC->0.8173522332114401`,eN->0.8754887942161612`},{δB->0.8`,A->0.1333333333333333`,B->0.03703703703703703`,q0A->0.8030165557168634`,q1A->0.9348610369021412`,q0B->0.8104868805099105`,q1B->0.9702291593382285`,eC->0.8134936068295078`,eN->0.8669882507439038`},{δB->0.82`,A->0.1333333333333333`,B->0.03703703703703703`,q0A->0.8029799814979253`,q1A->0.9353981968509071`,q0B->0.8097103468812298`,q1B->0.966336195894209`,eC->0.80981806845095`,eN->0.8585213783994428`},{δB->0.84`,A->0.1333333333333333`,B->0.04444444444444443`,q0A->0.8018027357554184`,q1A->0.9396777143164593`,q0B->0.8081134280824854`,q1B->0.9680187235874522`,eC->0.7767535331920957`,eN->0.8272714667851037`},{δB->0.86`,A->0.1333333333333333`,B->0.04444444444444443`,q0A->0.801781413724912`,q1A->0.9403066430872632`,q0B->0.8073133788578538`,q1B->0.9644988307832447`,eC->0.7729398284806481`,eN->0.817527642084486`},{δB->0.88`,A->0.1333333333333333`,B->0.051851851851851836`,q0A->0.8006263797143762`,q1A->0.9448766441189886`,q0B->0.8056269914536013`,q1B->0.966314715042733`,eC->0.7412306549766112`,eN->0.7844805065811934`},{δB->0.8999999999999999`,A->0.1333333333333333`,B->0.051851851851851836`,q0A->0.8006265641767434`,q1A->0.9455992187802945`,q0B->0.8048035511553278`,q1B->0.9630819945911516`,eC->0.7373561827219067`,eN->0.773619381879687`},{δB->0.9199999999999999`,A->0.1333333333333333`,B->0.051851851851851836`,q0A->0.8006374894707746`,q1A->0.9463111284590879`,q0B->0.8039876262376371`,q1B->0.9600146234314751`,eC->0.7336306634887131`,eN->0.7628096195614931`},{δB->0.94`,A->0.1333333333333333`,B->0.05925925925925924`,q0A->0.7995374731967974`,q1A->0.9512926047931889`,q0B->0.8021866400208117`,q1B->0.9619642697082452`,eC->0.7031712326742737`,eN->0.7273375500113867`},{δB->0.96`,A->0.1333333333333333`,B->0.05925925925925924`,q0A->0.7995769751667025`,q1A->0.9521040224034659`,q0B->0.801348091269551`,q1B->0.9590920830785914`,eC->0.699425090150437`,eN->0.7155939118376178`},{δB->0.98`,A->0.1333333333333333`,B->0.05925925925925924`,q0A->0.7996277378673693`,q1A->0.9529073850941249`,q0B->0.8005159091231712`,q1B->0.9563425109806407`,eC->0.6958070235656003`,eN->0.7039183337749143`},{δB->1.`,A->0.1333333333333333`,B->0.06666666666666665`,q0A->0.7986111111111115`,q1A->0.9583333333333323`,q0B->0.7986111111111104`,q1B->0.9583333333333353`,eC->0.6666666666666722`,eN->0.6666666666666671`}};

Clear[semistrategicresultsA2]
semistrategicresultsA2=Table[Map[{δB->aB}~Join~#~Join~ge[{Log,Exp,.5,.5,1,aB,.0,1,A,B,.1,13}/.#]&,sseq[aB,.0,.2,36]][[1]],{aB,.6,1,.02}];

Print["semistrategicresults with A=Abar/5 saved here"]
semistrategicresultsA2={{δB->0.6`,A->0.05333333333333332`,B->0.07407407407407406`,q0A->0.810797372145997`,q1A->0.8959752019817699`,q0B->0.8202585641467557`,q1B->0.9436723818958463`,eC->0.5105728803952171`,eN->0.6738510881448361`},{δB->0.62`,A->0.05333333333333332`,B->0.07407407407407406`,q0A->0.8106331933187158`,q1A->0.8968911898460022`,q0B->0.8197139003334583`,q1B->0.9414487236169963`,eC->0.503999481569463`,eN->0.6585600090854645`},{δB->0.64`,A->0.05333333333333332`,B->0.07407407407407406`,q0A->0.8104841858344195`,q1A->0.8977955931816535`,q0B->0.819171938570614`,q1B->0.9393473693439892`,eC->0.497741867525425`,eN->0.6435573109697241`},{δB->0.6599999999999999`,A->0.05333333333333332`,B->0.08148148148148146`,q0A->0.8091388802265168`,q1A->0.9017201543697279`,q0B->0.8179238156319266`,q1B->0.9427566727070608`,eC->0.468809192608316`,eN->0.6102925504622198`},{δB->0.6799999999999999`,A->0.05333333333333332`,B->0.08148148148148146`,q0A->0.8090087901146777`,q1A->0.9026974985053757`,q0B->0.8173551698194819`,q1B->0.9408100404513718`,eC->0.46308821271790157`,eN->0.5953627690679814`},{δB->0.7`,A->0.05333333333333332`,B->0.08148148148148146`,q0A->0.8088933180739116`,q1A->0.9036657203395642`,q0B->0.8167892075382746`,q1B->0.9389513617342241`,eC->0.4576161563272693`,eN->0.5807498165081479`},{δB->0.72`,A->0.05333333333333332`,B->0.08148148148148146`,q0A->0.8087918310726134`,q1A->0.9046257401116939`,q0B->0.816225835919663`,q1B->0.9371718285645442`,eC->0.4523731525990226`,eN->0.5664480837114851`},{δB->0.74`,A->0.05333333333333332`,B->0.08148148148148146`,q0A->0.8087037514034526`,q1A->0.9055784024558072`,q0B->0.8156649767499529`,q1B->0.9354638239995104`,eC->0.4473414792301377`,eN->0.5524518198474814`},{δB->0.76`,A->0.05333333333333332`,B->0.08888888888888886`,q0A->0.8073931489129973`,q1A->0.9099846514423232`,q0B->0.8142561568901968`,q1B->0.9389888032475864`,eC->0.4217952314300587`,eN->0.5201607707011642`},{δB->0.78`,A->0.05333333333333332`,B->0.08888888888888886`,q0A->0.8073282183742216`,q1A->0.9110187452550322`,q0B->0.8136701372306328`,q1B->0.9373576117990301`,eC->0.41715136840282735`,eN->0.5064661214667889`},{δB->0.8`,A->0.05333333333333332`,B->0.08888888888888886`,q0A->0.807276425661426`,q1A->0.9120479749769939`,q0B->0.813086766574023`,q1B->0.9357810390978222`,eC->0.41267639612566276`,eN->0.49308611868022034`},{δB->0.82`,A->0.05333333333333332`,B->0.08888888888888886`,q0A->0.807237305900468`,q1A->0.9130730888658677`,q0B->0.8125060026066098`,q1B->0.9342545091629758`,eC->0.408358163733386`,eN->0.4800132682042657`},{δB->0.84`,A->0.05333333333333332`,B->0.08888888888888886`,q0A->0.8072104278895558`,q1A->0.9140947967732403`,q0B->0.8119278113032077`,q1B->0.9327739680140934`,eC->0.40418563034653787`,eN->0.4672401710679293`},{δB->0.86`,A->0.05333333333333332`,B->0.09629629629629627`,q0A->0.8059706426961641`,q1A->0.9190232027783395`,q0B->0.8103690100726784`,q1B->0.9362779132761558`,eC->0.3815056206818701`,eN->0.43690577408850045`},{δB->0.88`,A->0.05333333333333332`,B->0.09629629629629627`,q0A->0.8059711599468898`,q1A->0.9201392174695261`,q0B->0.8097675309735168`,q1B->0.9348333698956033`,eC->0.37761127284028445`,eN->0.42455749145546834`},{δB->0.8999999999999999`,A->0.05333333333333332`,B->0.09629629629629627`,q0A->0.8059838353886913`,q1A->0.921254836637366`,q0B->0.8091688685812356`,q1B->0.9334254633982189`,eC->0.3738335747861965`,eN->0.41250635771436533`},{δB->0.9199999999999999`,A->0.05333333333333332`,B->0.09629629629629627`,q0A->0.8060083107478139`,q1A->0.9223707835420365`,q0B->0.8085730161081334`,q1B->0.9320515622335549`,eC->0.3701645424466632`,eN->0.4007443784850348`},{δB->0.94`,A->0.05333333333333332`,B->0.09629629629629627`,q0A->0.8060442500374023`,q1A->0.9234877660747921`,q0B->0.8079799730756284`,q1B->0.9307092951775582`,eC->0.3665968097336958`,eN->0.38926377410547064`},{δB->0.96`,A->0.05333333333333332`,B->0.10370370370370367`,q0A->0.8049102361863971`,q1A->0.9290142628806134`,q0B->0.8062825846459538`,q1B->0.9341154577438248`,eC->0.3462932921869143`,eN->0.36135932632717016`},{δB->0.98`,A->0.05333333333333332`,B->0.10370370370370367`,q0A->0.8049781523453436`,q1A->0.9302452949362386`,q0B->0.8056685162045361`,q1B->0.9327854696438299`,eC->0.34291802785959596`,eN->0.3503549359394655`},{δB->1.`,A->0.05333333333333332`,B->0.10370370370370367`,q0A->0.8050576184379001`,q1A->0.9314814814814814`,q0B->0.8050576184379001`,q1B->0.9314814814814815`,eC->0.33962264150943416`,eN->0.33962264150943405`}};

Clear[semistrategicresultsA5irs1]
semistrategicresultsA5irs1=Table[Map[{δB->aB}~Join~#~Join~ge[{Log,Exp,.5,.5,1,aB,.1,1,A,B,.1,13}/.#]&,sseq[aB,.1,.5,36]][[1]],{aB,.6,1,.02}];

Print["semistrategicresultsA5irs1 (with A=Abar/2 and rho=.1) saved here"]
semistrategicresultsA5irs1={{δB->0.6`,A->0.1333333333333333`,B->0.007407407407407405`,q0A->0.8088235294117646`,q1A->0.9166666666666666`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.62`,A->0.1333333333333333`,B->0.007407407407407405`,q0A->0.8088235294117646`,q1A->0.9166666666666666`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.64`,A->0.1333333333333333`,B->0.007407407407407405`,q0A->0.8088235294117646`,q1A->0.9166666666666666`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.6599999999999999`,A->0.1333333333333333`,B->0.007407407407407405`,q0A->0.8088235294117646`,q1A->0.9166666666666666`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.6799999999999999`,A->0.1333333333333333`,B->0.007407407407407405`,q0A->0.8088235294117646`,q1A->0.9166666666666666`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.7`,A->0.1333333333333333`,B->0.007407407407407405`,q0A->0.8088235294117646`,q1A->0.9166666666666666`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.72`,A->0.1333333333333333`,B->0.007407407407407405`,q0A->0.8088235294117646`,q1A->0.9166666666666666`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.74`,A->0.1333333333333333`,B->0.007407407407407405`,q0A->0.8088235294117646`,q1A->0.9166666666666666`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.76`,A->0.1333333333333333`,B->0.007407407407407405`,q0A->0.8088235294117646`,q1A->0.9166666666666666`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.78`,A->0.1333333333333333`,B->0.007407407407407405`,q0A->0.8088235294117646`,q1A->0.9166666666666666`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.8`,A->0.1333333333333333`,B->0.007407407407407405`,q0A->0.8088235294117646`,q1A->0.9166666666666666`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.82`,A->0.1333333333333333`,B->0.007407407407407405`,q0A->0.8088235294117646`,q1A->0.9166666666666666`,q0B->0.8333333333333333`,q1B->1.`,eC->1.`,eN->1},{δB->0.84`,A->0.1333333333333333`,B->0.02962962962962962`,q0A->0.8048622051311571`,q1A->0.9285917376262631`,q0B->0.8164387515688378`,q1B->0.9987836622592229`,eC->0.8689696712449582`,eN->0.9344557061414017`},{δB->0.86`,A->0.1333333333333333`,B->0.03703703703703703`,q0A->0.8037323153671435`,q1A->0.9327630859619567`,q0B->0.8143326699930828`,q1B->0.9923051564195781`,eC->0.8323691658377038`,eN->0.9048166929442458`},{δB->0.88`,A->0.1333333333333333`,B->0.04444444444444443`,q0A->0.8026173017141801`,q1A->0.9372285331048951`,q0B->0.8122829269928269`,q1B->0.988385784064192`,eC->0.7969419063033156`,eN->0.8724707787799882`},{δB->0.8999999999999999`,A->0.1333333333333333`,B->0.051851851851851836`,q0A->0.80152899996237`,q1A->0.941997719886056`,q0B->0.8102321769093792`,q1B->0.9858778860981587`,eC->0.7627734958982711`,eN->0.8377214709009188`},{δB->0.9199999999999999`,A->0.1333333333333333`,B->0.051851851851851836`,q0A->0.8015001061747501`,q1A->0.9429132258915736`,q0B->0.8093202471114987`,q1B->0.9810996521175924`,eC->0.7574936391600812`,eN->0.8254681541974701`},{δB->0.94`,A->0.1333333333333333`,B->0.05925925925925924`,q0A->0.8004777606563215`,q1A->0.9481249621045678`,q0B->0.8072188194910467`,q1B->0.9798307068551827`,eC->0.7245511100869482`,eN->0.7872369283713505`},{δB->0.96`,A->0.1333333333333333`,B->0.21481481481481474`,q0A->0.8333333333333333`,q1A->1.`,q0B->0.7979674796747968`,q1B->0.9693827160493828`,eC->0.`,eN->0},{δB->0.98`,A->0.1333333333333333`,B->0.17777777777777773`,q0A->0.8333333333333333`,q1A->1.`,q0B->0.8025157232704403`,q1B->0.9451851851851852`,eC->0.`,eN->0},{δB->1.`,A->0.1333333333333333`,B->0.17777777777777773`,q0A->0.8333333333333333`,q1A->1.`,q0B->0.8018867924528302`,q1B->0.9444444444444445`,eC->0.`,eN->0}};

Clear[semistrategicresultsA2irs1]
semistrategicresultsA2irs1=Table[Map[{δB->aB}~Join~#~Join~ge[{Log,Exp,.5,.5,1,aB,.1,1,A,B,.1,13}/.#]&,sseq[aB,.1,.2,36]][[1]],{aB,.6,1,.02}];

Print["semistrategicresultsA2irs1 (with A=Abar/5 and rho=.1) saved here"]
semistrategicresultsA2irs1={{δB->0.6`,A->0.05333333333333332`,B->0.08148148148148146`,q0A->0.8111435785891354`,q1A->0.8989724767845375`,q0B->0.8208286566094095`,q1B->0.9526320582327451`,eC->0.49552676955724084`,eN->0.6711871896155437`},{δB->0.62`,A->0.05333333333333332`,B->0.08148148148148146`,q0A->0.8109897809932924`,q1A->0.9001953901374703`,q0B->0.8202366953783254`,q1B->0.9499779515991609`,eC->0.48769847377479303`,eN->0.6528501511635172`},{δB->0.64`,A->0.05333333333333332`,B->0.08148148148148146`,q0A->0.8108562330943676`,q1A->0.9014060740030282`,q0B->0.8196489770636504`,q1B->0.9474914125100757`,eC->0.48028119720184004`,eN->0.6348786526037103`},{δB->0.6599999999999999`,A->0.05333333333333332`,B->0.08148148148148146`,q0A->0.8107417440246132`,q1A->0.9026063369349798`,q0B->0.81906527457822`,q1B->0.9451509282112256`,eC->0.4732334988660229`,eN->0.6172712627541885`},{δB->0.6799999999999999`,A->0.05333333333333332`,B->0.08888888888888886`,q0A->0.8095882643162605`,q1A->0.907344532964529`,q0B->0.8177085405052646`,q1B->0.9478966130261902`,eC->0.44415328281821814`,eN->0.5784118186257862`},{δB->0.7`,A->0.05333333333333332`,B->0.08888888888888886`,q0A->0.8095100323075332`,q1A->0.9086567806536403`,q0B->0.8170992560626795`,q1B->0.9457532032008616`,eC->0.43774954620707973`,eN->0.5609625143579717`},{δB->0.72`,A->0.05333333333333332`,B->0.08888888888888886`,q0A->0.8094497070861235`,q1A->0.9099634534181452`,q0B->0.8164941572436578`,q1B->0.94371108032572`,eC->0.4316244652543939`,eN->0.5439208341053248`},{δB->0.74`,A->0.05333333333333332`,B->0.08888888888888886`,q0A->0.8094064207821484`,q1A->0.9112659338020701`,q0B->0.8158931542520992`,q1B->0.9417594904956639`,eC->0.4257531485322248`,eN->0.5272808477292453`},{δB->0.76`,A->0.05333333333333332`,B->0.09629629629629627`,q0A->0.808350589568639`,q1A->0.9166370327806719`,q0B->0.8144043532225045`,q1B->0.9447050203211271`,eC->0.3999187417342084`,eN->0.48954609428323226`},{δB->0.78`,A->0.05333333333333332`,B->0.09629629629629627`,q0A->0.8083476130835697`,q1A->0.9180722525844225`,q0B->0.8137806304769303`,q1B->0.9428620171866764`,eC->0.394508659379043`,eN->0.47335114107216236`},{δB->0.8`,A->0.05333333333333332`,B->0.09629629629629627`,q0A->0.8083609286708521`,q1A->0.9195088180597262`,q0B->0.8131614061971116`,q1B->0.9410857554575078`,eC->0.3892916006246651`,eN->0.45757602464711716`},{δB->0.82`,A->0.05333333333333332`,B->0.09629629629629627`,q0A->0.8083898666216153`,q1A->0.9209480674427133`,q0B->0.8125466567288526`,q1B->0.9393701872984465`,eC->0.3842512993577671`,eN->0.44221149774270746`},{δB->0.84`,A->0.05333333333333332`,B->0.09629629629629627`,q0A->0.8084338102603843`,q1A->0.9223913040855475`,q0B->0.8119363723061752`,q1B->0.9377100153335163`,eC->0.3793730766831347`,eN->0.4272482840575331`},{δB->0.86`,A->0.05333333333333332`,B->0.10370370370370367`,q0A->0.807539219005124`,q1A->0.9286577739193447`,q0B->0.8103053215998146`,q1B->0.9407066043051487`,eC->0.3563756961643189`,eN->0.3918318854296145`},{δB->0.88`,A->0.05333333333333332`,B->0.20740740740740735`,q0A->0.8333333333333333`,q1A->1.`,q0B->0.801840490797546`,q1B->0.968148148148148`,eC->0.`,eN->0},{δB->0.8999999999999999`,A->0.05333333333333332`,B->0.19999999999999993`,q0A->0.8333333333333333`,q1A->1.`,q0B->0.8020833333333334`,q1B->0.9624999999999999`,eC->0.`,eN->0},{δB->0.9199999999999999`,A->0.05333333333333332`,B->0.19259259259259254`,q0A->0.8333333333333333`,q1A->1.`,q0B->0.8023809523809524`,q1B->0.9569135802469135`,eC->0.`,eN->0},{δB->0.94`,A->0.05333333333333332`,B->0.19259259259259254`,q0A->0.8333333333333333`,q1A->1.`,q0B->0.8017080745341615`,q1B->0.9561111111111111`,eC->0.`,eN->0},{δB->0.96`,A->0.05333333333333332`,B->0.18518518518518512`,q0A->0.8333333333333333`,q1A->1.`,q0B->0.8020833333333334`,q1B->0.9506172839506173`,eC->0.`,eN->0},{δB->0.98`,A->0.05333333333333332`,B->0.17777777777777773`,q0A->0.8333333333333333`,q1A->1.`,q0B->0.8025157232704403`,q1B->0.9451851851851852`,eC->0.`,eN->0},{δB->1.`,A->0.05333333333333332`,B->0.17777777777777773`,q0A->0.8333333333333333`,q1A->1.`,q0B->0.8018867924528302`,q1B->0.9444444444444445`,eC->0.`,eN->0}};

Figures for the paper

With ρ=0, large A:

Module[{ρ=.0,u=Log,uinv=Exp,upinv,ℓ=.5,λ=.5,δA=1,(*δB=.98,ρ=0.0,*)M=1,i=.1,prec=10,Abar},
upinv=z/.Solve[u'[z]==#,z][[1]]&;
Abar=M(1/upinv[1+i/(ℓ-λ ℓ(1-ℓ)δA)]-1);
ListPlot[Transpose[{{δB,A 100/Abar},{δB,If[eC<1,B ,Null]100/Abar}}/.semistrategicresultsA5],PlotRange->{{.6,1},{0,98}},AspectRatio->1,ImageSize->250,Joined->True,AxesLabel->{" (Subscript[δ, B])","(% of (OverBar[A]))"},BaseStyle->{FontSize->18},PlotStyle->{{Thick,Darker@Blue},{Thick,Dashed,Darker@Red}},
Epilog->{Text[Style["A",Darker@Blue,FontSize->20],{.8,.16 100/Abar}],Text[Style["B",Darker@Red,FontSize->20],{.8,.06 100/Abar}]
}]]
Export["csss-crs-alarge-assetsupply.pdf",%]

Module[{ρ=.0,u=Log,uinv=Exp,upinv,ℓ=.5,λ=.5,δA=1,(*δB=.98,ρ=0.0,*)M=1,i=.1,prec=10,Abar},
upinv=z/.Solve[u'[z]==#,z][[1]]&;
Abar=M(1/upinv[1+i/(ℓ-λ ℓ(1-ℓ)δA)]-1);
ListPlot[Transpose[{{δB,100If[eN>0,λ ℓ δA (((1+(ℓ/(1-ℓ))(eC/eN)))^(ρ-1)) ((((1-ℓ)eN))^ρ)(u'[q1A]-1),Null]},{δB,100If[eN<1,λ ℓ δB (((1+(ℓ/(1-ℓ))((1-eC)/(1-eN))))^(ρ-1)) ((((1-ℓ)(1-eN)))^ρ)(u'[q1B]-1),0]}}/.semistrategicresultsA5],PlotRange->{{.6,1},{0,100λ ℓ(1-ℓ)δA i}},AspectRatio->1,ImageSize->250,Joined->True,AxesLabel->{" (Subscript[δ, B])","(%)"},BaseStyle->{FontSize->18},PlotStyle->{{Thick,Darker@Blue},{Thick,Dashed,Darker@Red}},
Epilog->{Text[Style["(Subscript[L, A])",Darker@Blue,FontSize->20],{.9,.9}],Text[Style["(Subscript[L, B])",Darker@Red,FontSize->20],{.9,.2}]
}]]
Export["csss-crs-alarge-premia.pdf",%]

Module[{ρ=.0,u=Log,uinv=Exp,upinv,ℓ=.5,λ=.5,δA=1,(*δB=.98,ρ=0.0,*)M=1,i=.1,prec=10,Abar},
upinv=z/.Solve[u'[z]==#,z][[1]]&;
Abar=M(1/upinv[1+i/(ℓ-λ ℓ(1-ℓ)δA)]-1);
ListPlot[Transpose[{{δB,If[eN>0, δA (((1+(ℓ/(1-ℓ))(eC/eN)))^(ρ-1)) ((((1-ℓ)eN))^ρ),0]},{δB,If[eN<1,δB (((1+(ℓ/(1-ℓ))((1-eC)/(1-eN))))^(ρ-1)) ((((1-ℓ)(1-eN)))^ρ),0]},{δB,(1-ℓ)δB},{δB,(1-ℓ)δA}}/.semistrategicresultsA5],PlotRange->{{.6,1},{.31,.59}},AspectRatio->1,ImageSize->250,Joined->True,AxesLabel->{" (Subscript[δ, B])",None},BaseStyle->{FontSize->18},PlotStyle->{{Thick,Darker@Blue},{Thick,Dashed,Darker@Red},{Thin,Dashed,Gray},{Thin,Dashed,Gray}},
Epilog->{Text[Style["(Subscript[α, CA])",Darker@Blue,FontSize->20],{.8,.54}],Text[Style["(Subscript[α, CB])",Darker@Red,FontSize->20],{.9,.38}],Text[Style["(1-ℓ)(Subscript[δ, A])",Darker@Gray,FontSize->16],{.86,.48}],Text[Style["(1-ℓ)(Subscript[δ, B])",Darker@Gray,FontSize->16],{.69,.38}]
}]]
Export["csss-crs-alarge-arrival.pdf",%]

Module[{ρ=.0,u=Log,uinv=Exp,upinv,ℓ=.5,λ=.5,δA=1,(*δB=.98,ρ=0.0,*)M=1,i=.1,prec=10,Abar,fA,fB,q1Atil,q1Btil},
upinv=z/.Solve[u'[z]==#,z][[1]]&;
Abar=M(1/upinv[1+i/(ℓ-λ ℓ(1-ℓ)δA)]-1);
fA=If[eC>0,δA((((1/(ℓ eC))+(1/((1-ℓ)eN))))^(-1)) (((ℓ eC+(1-ℓ)eN))^ρ),0];
fB=If[eC<1,δB((((1/(ℓ(1-eC)))+(1/((1-ℓ)(1-eN)))))^(-1))(((ℓ (1-eC)+(1-ℓ)(1-eN)))^ρ),0];
q1Atil=If[eC>0,Min[1,uinv[u[q0A]+(A/(eC M))(eC q0A+(1-eC)q0B)]],1];
q1Btil=If[eC<1,Min[1,uinv[u[q0B]+(B/((1-eC) M))(eC q0A+(1-eC)q0B)]],1];
ListPlot[Transpose[{{δB,100q0A},{δB,100q1A},{δB,100If[eC<1,q0B,Null]},{δB,100If[eC<1,q1B,Null]},{δB,100(fA(λ q1A+(1-λ)q1Atil)+fB(λ q1B+(1-λ)q1Btil)+(ℓ eC-fA)q0A+(ℓ(1-eC)-fB)q0B)/ℓ}}/.semistrategicresultsA5],
PlotRange->{{.6,1},{77,101}},AspectRatio->1,ImageSize->250,Joined->True,AxesLabel->{" (Subscript[δ, B])","(% of (SuperStar[q]))"},BaseStyle->{FontSize->18},PlotStyle->{{Thick,Dashed,Darker@Blue},{Thick,Darker@Blue},{Thick,Dashed,Darker@Red},{Thick,Darker@Red},{Thick,DotDashed,RGBColor[.33,0,.33]}},
Epilog->{Text[Style["(Subscript[q, 0  A])",Darker@Blue,FontSize->18],{.7,79}],Text[Style["(Subscript[q, 1  A])",Darker@Blue,FontSize->18],{.7,94}],Text[Style["(Subscript[q, 0  B])",Darker@Red,FontSize->18],{.9,82.5}],Text[Style["(Subscript[q, 1  B])",Darker@Red,FontSize->18],{.9,98}],Text[Style["(OverBar[q])",RGBColor[.33,0,.33],FontSize->18],{.8,88}]
}]]
Export["csss-crs-alarge-allqs.pdf",%]

With ρ=0, small A:

Module[{ρ=.0,u=Log,uinv=Exp,upinv,ℓ=.5,λ=.5,δA=1,(*δB=.98,ρ=0.0,*)M=1,i=.1,prec=10,Abar},
upinv=z/.Solve[u'[z]==#,z][[1]]&;
Abar=M(1/upinv[1+i/(ℓ-λ ℓ(1-ℓ)δA)]-1);
ListPlot[Transpose[{{δB,A 100/Abar},{δB,If[eC<1,B ,Null]100/Abar}}/.semistrategicresultsA2],PlotRange->{{.6,1},{0,98}},AspectRatio->1,ImageSize->250,Joined->True,AxesLabel->{" (Subscript[δ, B])","(% of (OverBar[A]))"},BaseStyle->{FontSize->18},PlotStyle->{{Thick,Darker@Blue},{Thick,Dashed,Darker@Red}},
Epilog->{Text[Style["A",Darker@Blue,FontSize->20],{.8,.03 100/Abar}],Text[Style["B",Darker@Red,FontSize->20],{.8,.11 100/Abar}]
}]]
Export["csss-crs-asmall-assetsupply.pdf",%]

Module[{ρ=.0,u=Log,uinv=Exp,upinv,ℓ=.5,λ=.5,δA=1,(*δB=.98,ρ=0.0,*)M=1,i=.1,prec=10,Abar},
upinv=z/.Solve[u'[z]==#,z][[1]]&;
Abar=M(1/upinv[1+i/(ℓ-λ ℓ(1-ℓ)δA)]-1);
ListPlot[Transpose[{{δB,100If[eN>0,λ ℓ δA (((1+(ℓ/(1-ℓ))(eC/eN)))^(ρ-1)) ((((1-ℓ)eN))^ρ)(u'[q1A]-1),Null]},{δB,100If[eN<1,λ ℓ δB (((1+(ℓ/(1-ℓ))((1-eC)/(1-eN))))^(ρ-1)) ((((1-ℓ)(1-eN)))^ρ)(u'[q1B]-1),0]}}/.semistrategicresultsA2],PlotRange->{{.6,1},{0,150λ ℓ(1-ℓ)δA i}},AspectRatio->1,ImageSize->250,Joined->True,AxesLabel->{" (Subscript[δ, B])","(%)"},BaseStyle->{FontSize->18},PlotStyle->{{Thick,Darker@Blue},{Thick,Dashed,Darker@Red}},
Epilog->{Text[Style["(Subscript[L, A])",Darker@Blue,FontSize->20],{.9,1.3}],Text[Style["(Subscript[L, B])",Darker@Red,FontSize->20],{.9,.6}]
}]]
Export["csss-crs-asmall-premia.pdf",%]

Module[{ρ=.0,u=Log,uinv=Exp,upinv,ℓ=.5,λ=.5,δA=1,(*δB=.98,ρ=0.0,*)M=1,i=.1,prec=10,Abar},
upinv=z/.Solve[u'[z]==#,z][[1]]&;
Abar=M(1/upinv[1+i/(ℓ-λ ℓ(1-ℓ)δA)]-1);
ListPlot[Transpose[{{δB,If[eN>0, δA (((1+(ℓ/(1-ℓ))(eC/eN)))^(ρ-1)) ((((1-ℓ)eN))^ρ),0]},{δB,If[eN<1,δB (((1+(ℓ/(1-ℓ))((1-eC)/(1-eN))))^(ρ-1)) ((((1-ℓ)(1-eN)))^ρ),0]},{δB,(1-ℓ)δB},{δB,(1-ℓ)δA}}/.semistrategicresultsA2],PlotRange->{{.6,1},{.31,.59}},AspectRatio->1,ImageSize->250,Joined->True,AxesLabel->{" (Subscript[δ, B])",None},BaseStyle->{FontSize->18},PlotStyle->{{Thick,Darker@Blue},{Thick,Dashed,Darker@Red},{Thin,Dashed,Gray},{Thin,Dashed,Gray}},
Epilog->{Text[Style["(Subscript[α, CA])",Darker@Blue,FontSize->20],{.8,.57}],Text[Style["(Subscript[α, CB])",Darker@Red,FontSize->20],{.9,.4}],Text[Style["(1-ℓ)(Subscript[δ, A])",Darker@Gray,FontSize->16],{.86,.48}],Text[Style["(1-ℓ)(Subscript[δ, B])",Darker@Gray,FontSize->16],{.69,.38}]
}]]
Export["csss-crs-asmall-arrival.pdf",%]

Module[{ρ=.0,u=Log,uinv=Exp,upinv,ℓ=.5,λ=.5,δA=1,(*δB=.98,ρ=0.0,*)M=1,i=.1,prec=10,Abar,fA,fB,q1Atil,q1Btil},
upinv=z/.Solve[u'[z]==#,z][[1]]&;
Abar=M(1/upinv[1+i/(ℓ-λ ℓ(1-ℓ)δA)]-1);
fA=If[eC>0,δA((((1/(ℓ eC))+(1/((1-ℓ)eN))))^(-1)) ,0];
fB=If[eC<1,δB((((1/(ℓ(1-eC)))+(1/((1-ℓ)(1-eN)))))^(-1)),0];
q1Atil=If[eC>0,Min[1,uinv[u[q0A]+(A/(eC M))(eC q0A+(1-eC)q0B)]],1];
q1Btil=If[eC<1,Min[1,uinv[u[q0B]+(B/((1-eC) M))(eC q0A+(1-eC)q0B)]],1];
ListPlot[Transpose[{{δB,100q0A},{δB,100q1A},{δB,100If[eC<1,q0B,Null]},{δB,100If[eC<1,q1B,Null]},{δB,100(fA(λ q1A+(1-λ)q1Atil)+fB(λ q1B+(1-λ)q1Btil)+(ℓ eC-fA)q0A+(ℓ(1-eC)-fB)q0B)/ℓ}}/.semistrategicresultsA2],
PlotRange->{{.6,1},{77,101}},AspectRatio->1,ImageSize->250,Joined->True,AxesLabel->{" (Subscript[δ, B])","(% of (SuperStar[q]))"},BaseStyle->{FontSize->18},PlotStyle->{{Thick,Dashed,Darker@Blue},{Thick,Darker@Blue},{Thick,Dashed,Darker@Red},{Thick,Darker@Red},{Thick,DotDashed,RGBColor[.33,0,.33]}},
Epilog->{Text[Style["(Subscript[q, 0  A])",Darker@Blue,FontSize->18],{.7,80}],Text[Style["(Subscript[q, 1  A])",Darker@Blue,FontSize->18],{.7,92}],Text[Style["(Subscript[q, 0  B])",Darker@Red,FontSize->18],{.9,82.5}],Text[Style["(Subscript[q, 1  B])",Darker@Red,FontSize->18],{.9,96}],Text[Style["(OverBar[q])",RGBColor[.33,0,.33],FontSize->18],{.8,87.5}]
}]]
Export["csss-crs-asmall-allqs.pdf",%]

With ρ=.1 and large A:

Module[{ρ=.1,u=Log,uinv=Exp,upinv,ℓ=.5,λ=.5,δA=1,(*δB=.98,ρ=0.0,*)M=1,i=.1,prec=10,Abar},
upinv=z/.Solve[u'[z]==#,z][[1]]&;
Abar=M(1/upinv[1+i/(ℓ-λ ℓ(1-ℓ)δA)]-1);
ListPlot[Transpose[{{δB,A 100/Abar},{δB,If[eC<1,B ,Null]100/Abar}}/.semistrategicresultsA5irs1],PlotRange->{{.6,1},{0,98}},AspectRatio->1,ImageSize->250,Joined->True,AxesLabel->{" (Subscript[δ, B])","(% of (OverBar[A]))"},BaseStyle->{FontSize->18},PlotStyle->{{Thick,Darker@Blue},{Thick,Dashed,Darker@Red}},
Epilog->{Text[Style["A",Darker@Blue,FontSize->20],{.87,.16 100/Abar}],Text[Style["B",Darker@Red,FontSize->20],{.87,.06 100/Abar}]
}]]
Export["csss-irs-alarge-assetsupply.pdf",%]

Module[{ρ=.1,u=Log,uinv=Exp,upinv,ℓ=.5,λ=.5,δA=1,(*δB=.98,ρ=0.0,*)M=1,i=.1,prec=10,Abar},
upinv=z/.Solve[u'[z]==#,z][[1]]&;
Abar=M(1/upinv[1+i/(ℓ-λ ℓ(1-ℓ)δA)]-1);
ListPlot[Transpose[{{δB,100If[eN>0,λ ℓ δA (((1+(ℓ/(1-ℓ))(eC/eN)))^(ρ-1)) ((((1-ℓ)eN))^ρ)(u'[q1A]-1),0]},{δB,100If[eN<1,λ ℓ δB (((1+(ℓ/(1-ℓ))((1-eC)/(1-eN))))^(ρ-1)) ((((1-ℓ)(1-eN)))^ρ)(u'[q1B]-1),0]}}/.semistrategicresultsA5irs1],PlotRange->{{.6,1},{0,100λ ℓ(1-ℓ)δA i}},AspectRatio->1,ImageSize->250,Joined->True,AxesLabel->{" (Subscript[δ, B])","(%)"},BaseStyle->{FontSize->18},PlotStyle->{{Thick,Darker@Blue},{Thick,Dashed,Darker@Red}},
Epilog->{Text[Style["(Subscript[L, A])",Darker@Blue,FontSize->20],{.9,.95}],Text[Style["(Subscript[L, B])",Darker@Red,FontSize->20],{.9,.25}]
}]]
Export["csss-irs-alarge-premia.pdf",%]

Module[{ρ=.1,u=Log,uinv=Exp,upinv,ℓ=.5,λ=.5,δA=1,(*δB=.98,ρ=0.0,*)M=1,i=.1,prec=10,Abar},
upinv=z/.Solve[u'[z]==#,z][[1]]&;
Abar=M(1/upinv[1+i/(ℓ-λ ℓ(1-ℓ)δA)]-1);
ListPlot[Transpose[{{δB,If[eN>0, δA (((1+(ℓ/(1-ℓ))(eC/eN)))^(ρ-1)) ((((1-ℓ)eN))^ρ),0]},{δB,If[eN<1,δB (((1+(ℓ/(1-ℓ))((1-eC)/(1-eN))))^(ρ-1)) ((((1-ℓ)(1-eN)))^ρ),0]},{δB,(1-ℓ)δB},{δB,(1-ℓ)δA}}/.semistrategicresultsA5irs1],PlotRange->{{.6,1},{.31,.59}},AspectRatio->1,ImageSize->250,Joined->True,AxesLabel->{" (Subscript[δ, B])",None},BaseStyle->{FontSize->18},PlotStyle->{{Thick,Darker@Blue},{Thick,Dashed,Darker@Red},{Thin,Dashed,Gray},{Thin,Dashed,Gray}},
Epilog->{Text[Style["(Subscript[α, CA])",Darker@Blue,FontSize->20],{.82,.54}],Text[Style["(Subscript[α, CB])",Darker@Red,FontSize->20],{.88,.34}],Text[Style["(1-ℓ)(Subscript[δ, A])",Darker@Gray,FontSize->16],{.86,.48}],Text[Style["(1-ℓ)(Subscript[δ, B])",Darker@Gray,FontSize->16],{.69,.38}]
}]]
Export["csss-irs-alarge-arrival.pdf",%]

Module[{ρ=.1,u=Log,uinv=Exp,upinv,ℓ=.5,λ=.5,δA=1,(*δB=.98,ρ=0.0,*)M=1,i=.1,prec=10,Abar,fA,fB,q1Atil,q1Btil},
upinv=z/.Solve[u'[z]==#,z][[1]]&;
Abar=M(1/upinv[1+i/(ℓ-λ ℓ(1-ℓ)δA)]-1);
fA=If[eC>0,δA((((1/(ℓ eC))+(1/((1-ℓ)eN))))^(-1)) (((ℓ eC+(1-ℓ)eN))^ρ),0];
fB=If[eC<1,δB((((1/(ℓ(1-eC)))+(1/((1-ℓ)(1-eN)))))^(-1))(((ℓ (1-eC)+(1-ℓ)(1-eN)))^ρ),0];
q1Atil=If[eC>0,Min[1,uinv[u[q0A]+(A/(eC M))(eC q0A+(1-eC)q0B)]],1];
q1Btil=If[eC<1,Min[1,uinv[u[q0B]+(B/((1-eC) M))(eC q0A+(1-eC)q0B)]],1];
ListPlot[Transpose[{{δB,100If[eC>0,q0A,Null]},{δB,100If[eC>0,q1A,Null]},{δB,100If[eC<1,q0B,Null]},{δB,100If[eC<1,q1B,Null]},{δB,100(fA(λ q1A+(1-λ)q1Atil)+fB(λ q1B+(1-λ)q1Btil)+(ℓ eC-fA)q0A+(ℓ(1-eC)-fB)q0B)/ℓ}}/.semistrategicresultsA5irs1],
PlotRange->{{.6,1},{77,101}},AspectRatio->1,ImageSize->250,Joined->True,AxesLabel->{" (Subscript[δ, B])","(% of (SuperStar[q]))"},BaseStyle->{FontSize->18},PlotStyle->{{Thick,Dashed,Darker@Blue},{Thick,Darker@Blue},{Thick,Dashed,Darker@Red},{Thick,Darker@Red},{Thick,DotDashed,RGBColor[.33,0,.33]}},
Epilog->{Text[Style["(Subscript[q, 0  A])",Darker@Blue,FontSize->18],{.7,79}],Text[Style["(Subscript[q, 1  A])",Darker@Blue,FontSize->18],{.7,94}],Text[Style["(Subscript[q, 0  B])",Darker@Red,FontSize->18],{.9,83}],Text[Style["(Subscript[q, 1  B])",Darker@Red,FontSize->18],{.9,100.5}],Text[Style["(OverBar[q])",RGBColor[.33,0,.33],FontSize->18],{.8,88}]
}]]
Export["csss-irs-alarge-allqs.pdf",%]

With ρ=.1, small A:

Module[{ρ=.1,u=Log,uinv=Exp,upinv,ℓ=.5,λ=.5,δA=1,(*δB=.98,ρ=0.0,*)M=1,i=.1,prec=10,Abar},
upinv=z/.Solve[u'[z]==#,z][[1]]&;
Abar=M(1/upinv[1+i/(ℓ-λ ℓ(1-ℓ)δA)]-1);
ListPlot[Transpose[{{δB,A 100/Abar},{δB,If[eC<1,B ,Null]100/Abar}}/.semistrategicresultsA2irs1],PlotRange->{{.6,1},{0,98}},AspectRatio->1,ImageSize->250,Joined->True,AxesLabel->{" (Subscript[δ, B])","(% of (OverBar[A]))"},BaseStyle->{FontSize->18},PlotStyle->{{Thick,Darker@Blue},{Thick,Dashed,Darker@Red}},
Epilog->{Text[Style["A",Darker@Blue,FontSize->20],{.8,.03 100/Abar}],Text[Style["B",Darker@Red,FontSize->20],{.8,.12 100/Abar}]
}]]
Export["csss-irs-asmall-assetsupply.pdf",%]

Module[{ρ=.1,u=Log,uinv=Exp,upinv,ℓ=.5,λ=.5,δA=1,(*δB=.98,ρ=0.0,*)M=1,i=.1,prec=10,Abar},
upinv=z/.Solve[u'[z]==#,z][[1]]&;
Abar=M(1/upinv[1+i/(ℓ-λ ℓ(1-ℓ)δA)]-1);
ListPlot[Transpose[{{δB,100If[eN>0,λ ℓ δA (((1+(ℓ/(1-ℓ))(eC/eN)))^(ρ-1)) ((((1-ℓ)eN))^ρ)(u'[q1A]-1),0]},{δB,100If[eN<1,λ ℓ δB (((1+(ℓ/(1-ℓ))((1-eC)/(1-eN))))^(ρ-1)) ((((1-ℓ)(1-eN)))^ρ)(u'[q1B]-1),0]}}/.semistrategicresultsA2irs1],PlotRange->{{.6,1},{0,150λ ℓ(1-ℓ)δA i}},AspectRatio->1,ImageSize->250,Joined->True,AxesLabel->{" (Subscript[δ, B])","(%)"},BaseStyle->{FontSize->18},PlotStyle->{{Thick,Darker@Blue},{Thick,Dashed,Darker@Red}},
Epilog->{Text[Style["(Subscript[L, A])",Darker@Blue,FontSize->20],{.8,1.4}],Text[Style["(Subscript[L, B])",Darker@Red,FontSize->20],{.8,.4}]
}]]
Export["csss-irs-asmall-premia.pdf",%]

Module[{ρ=.1,u=Log,uinv=Exp,upinv,ℓ=.5,λ=.5,δA=1,(*δB=.98,ρ=0.0,*)M=1,i=.1,prec=10,Abar},
upinv=z/.Solve[u'[z]==#,z][[1]]&;
Abar=M(1/upinv[1+i/(ℓ-λ ℓ(1-ℓ)δA)]-1);
ListPlot[Transpose[{{δB,If[eN>0, δA (((1+(ℓ/(1-ℓ))(eC/eN)))^(ρ-1)) ((((1-ℓ)eN))^ρ),0]},{δB,If[eN<1,δB (((1+(ℓ/(1-ℓ))((1-eC)/(1-eN))))^(ρ-1)) ((((1-ℓ)(1-eN)))^ρ),0]},{δB,(1-ℓ)δB},{δB,(1-ℓ)δA}}/.semistrategicresultsA2irs1],PlotRange->{{.6,1},{.31,.59}},AspectRatio->1,ImageSize->250,Joined->True,AxesLabel->{" (Subscript[δ, B])",None},BaseStyle->{FontSize->18},PlotStyle->{{Thick,Darker@Blue},{Thick,Dashed,Darker@Red},{Thin,Dashed,Gray},{Thin,Dashed,Gray}},
Epilog->{Text[Style["(Subscript[α, CA])",Darker@Blue,FontSize->20],{.8,.54}],Text[Style["(Subscript[α, CB])",Darker@Red,FontSize->20],{.81,.33}],Text[Style["(1-ℓ)(Subscript[δ, A])",Darker@Gray,FontSize->16],{.68,.48}],Text[Style["(1-ℓ)(Subscript[δ, B])",Darker@Gray,FontSize->16],{.68,.38}]
}]]
Export["csss-irs-asmall-arrival.pdf",%]

Module[{ρ=.1,u=Log,uinv=Exp,upinv,ℓ=.5,λ=.5,δA=1,(*δB=.98,ρ=0.0,*)M=1,i=.1,prec=10,Abar,fA,fB,q1Atil,q1Btil},
upinv=z/.Solve[u'[z]==#,z][[1]]&;
Abar=M(1/upinv[1+i/(ℓ-λ ℓ(1-ℓ)δA)]-1);
fA=If[eC>0,δA((((1/(ℓ eC))+(1/((1-ℓ)eN))))^(-1)) (((ℓ eC+(1-ℓ)eN))^ρ),0];
fB=If[eC<1,δB((((1/(ℓ(1-eC)))+(1/((1-ℓ)(1-eN)))))^(-1))(((ℓ (1-eC)+(1-ℓ)(1-eN)))^ρ),0];
q1Atil=If[eC>0,Min[1,uinv[u[q0A]+(A/(eC M))(eC q0A+(1-eC)q0B)]],1];
q1Btil=If[eC<1,Min[1,uinv[u[q0B]+(B/((1-eC) M))(eC q0A+(1-eC)q0B)]],1];
ListPlot[Transpose[{{δB,100If[eC>0,q0A,Null]},{δB,100If[eC>0,q1A,Null]},{δB,100If[eC<1,q0B,Null]},{δB,100If[eC<1,q1B,Null]},{δB,100(fA(λ q1A+(1-λ)q1Atil)+fB(λ q1B+(1-λ)q1Btil)+(ℓ eC-fA)q0A+(ℓ(1-eC)-fB)q0B)/ℓ}}/.semistrategicresultsA2irs1],
PlotRange->{{.6,1},{77,101}},AspectRatio->1,ImageSize->250,Joined->True,AxesLabel->{" (Subscript[δ, B])","(% of (SuperStar[q]))"},BaseStyle->{FontSize->18},PlotStyle->{{Thick,Dashed,Darker@Blue},{Thick,Darker@Blue},{Thick,Dashed,Darker@Red},{Thick,Darker@Red},{Thick,DotDashed,RGBColor[.33,0,.33]}},
Epilog->{Text[Style["(Subscript[q, 0  A])",Darker@Blue,FontSize->18],{.7,80}],Text[Style["(Subscript[q, 1  A])",Darker@Blue,FontSize->18],{.7,92.5}],Text[Style["(Subscript[q, 0  B])",Darker@Red,FontSize->18],{.9,82.5}],Text[Style["(Subscript[q, 1  B])",Darker@Red,FontSize->18],{.9,99}],Text[Style["(OverBar[q])",RGBColor[.33,0,.33],FontSize->18],{.8,88}]
}]]
Export["csss-irs-asmall-allqs.pdf",%]

Welfare:

Module[{ρ=.0,u=Log,uinv=Exp,upinv,v,ℓ=.5,λ=.5,δA=1,(*δB=.98,ρ=0.0,*)M=1,i=.1,prec=10,Abar,fA,fB,q1Atil,q1Btil},
upinv=z/.Solve[u'[z]==#,z][[1]]&;
v=(u[#]-#+1)&;
Abar=M(1/upinv[1+i/(ℓ-λ ℓ(1-ℓ)δA)]-1);
fA=If[eC>0,δA((((1/(ℓ eC))+(1/((1-ℓ)eN))))^(-1)) ,0];
fB=If[eC<1,δB((((1/(ℓ(1-eC)))+(1/((1-ℓ)(1-eN)))))^(-1)),0];
q1Atil=If[eC>0,Min[1,uinv[u[q0A]+(A/(eC M))(eC q0A+(1-eC)q0B)]],1];
q1Btil=If[eC<1,Min[1,uinv[u[q0B]+(B/((1-eC) M))(eC q0A+(1-eC)q0B)]],1];
ListPlot[{{δB,100(fA(λ v[q1A]+(1-λ)v[q1Atil])+fB(λ v[q1B]+(1-λ)v[q1Btil])+(ℓ eC-fA)v[q0A]+(ℓ(1-eC)-fB)v[q0B])/ℓ}/.semistrategicresultsA5},
PlotRange->{{.6,1},{-1.33,-1.225}},AxesOrigin->{.6,-1.33},AspectRatio->1,ImageSize->260,Joined->True,AxesLabel->{" (Subscript[δ, B])","(%)"},BaseStyle->{FontSize->18},PlotStyle->{{Thick,Darker@Green}}]]
Export["csss-crs-alarge-welfare.pdf",%]

Module[{ρ=.0,u=Log,uinv=Exp,upinv,v,ℓ=.5,λ=.5,δA=1,(*δB=.98,ρ=0.0,*)M=1,i=.1,prec=10,Abar,fA,fB,q1Atil,q1Btil},
upinv=z/.Solve[u'[z]==#,z][[1]]&;
v=(u[#]-#+1)&;
Abar=M(1/upinv[1+i/(ℓ-λ ℓ(1-ℓ)δA)]-1);
fA=If[eC>0,δA((((1/(ℓ eC))+(1/((1-ℓ)eN))))^(-1)) ,0];
fB=If[eC<1,δB((((1/(ℓ(1-eC)))+(1/((1-ℓ)(1-eN)))))^(-1)),0];
q1Atil=If[eC>0,Min[1,uinv[u[q0A]+(A/(eC M))(eC q0A+(1-eC)q0B)]],1];
q1Btil=If[eC<1,Min[1,uinv[u[q0B]+(B/((1-eC) M))(eC q0A+(1-eC)q0B)]],1];
ListPlot[{{δB,100(fA(λ v[q1A]+(1-λ)v[q1Atil])+fB(λ v[q1B]+(1-λ)v[q1Btil])+(ℓ eC-fA)v[q0A]+(ℓ(1-eC)-fB)v[q0B])/ℓ}/.semistrategicresultsA2},
PlotRange->{{.6,1},{-1.33,-1.225}},AxesOrigin->{.6,-1.33},AspectRatio->1,ImageSize->260,Joined->True,AxesLabel->{" (Subscript[δ, B])","(%)"},BaseStyle->{FontSize->18},PlotStyle->{{Thick,Darker@Green}}]]
Export["csss-crs-asmall-welfare.pdf",%]

Module[{ρ=.1,u=Log,uinv=Exp,upinv,v,ℓ=.5,λ=.5,δA=1,(*δB=.98,ρ=0.0,*)M=1,i=.1,prec=10,Abar,fA,fB,q1Atil,q1Btil},
upinv=z/.Solve[u'[z]==#,z][[1]]&;
v=(u[#]-#+1)&;
Abar=M(1/upinv[1+i/(ℓ-λ ℓ(1-ℓ)δA)]-1);
fA=If[eC>0,δA((((1/(ℓ eC))+(1/((1-ℓ)eN))))^(-1)) (((ℓ eC+(1-ℓ)eN))^ρ),0];
fB=If[eC<1,δB((((1/(ℓ(1-eC)))+(1/((1-ℓ)(1-eN)))))^(-1))(((ℓ (1-eC)+(1-ℓ)(1-eN)))^ρ),0];
q1Atil=If[eC>0,Min[1,uinv[u[q0A]+(A/(eC M))(eC q0A+(1-eC)q0B)]],1];
q1Btil=If[eC<1,Min[1,uinv[u[q0B]+(B/((1-eC) M))(eC q0A+(1-eC)q0B)]],1];
ListPlot[{{δB,100(fA(λ v[q1A]+(1-λ)v[q1Atil])+fB(λ v[q1B]+(1-λ)v[q1Btil])+(ℓ eC-fA)v[q0A]+(ℓ(1-eC)-fB)v[q0B])/ℓ}/.semistrategicresultsA5irs1},
PlotRange->{{.6,1},{-1.33,-1.225}},AxesOrigin->{.6,-1.33},AspectRatio->1,ImageSize->260,Joined->True,AxesLabel->{" (Subscript[δ, B])","(%)"},BaseStyle->{FontSize->18},PlotStyle->{{Thick,Darker@Green}}]]
Export["csss-irs-alarge-welfare.pdf",%]

Module[{ρ=.1,u=Log,uinv=Exp,upinv,v,ℓ=.5,λ=.5,δA=1,(*δB=.98,ρ=0.0,*)M=1,i=.1,prec=10,Abar,fA,fB,q1Atil,q1Btil},
upinv=z/.Solve[u'[z]==#,z][[1]]&;
v=(u[#]-#+1)&;
Abar=M(1/upinv[1+i/(ℓ-λ ℓ(1-ℓ)δA)]-1);
fA=If[eC>0,δA((((1/(ℓ eC))+(1/((1-ℓ)eN))))^(-1)) (((ℓ eC+(1-ℓ)eN))^ρ),0];
fB=If[eC<1,δB((((1/(ℓ(1-eC)))+(1/((1-ℓ)(1-eN)))))^(-1))(((ℓ (1-eC)+(1-ℓ)(1-eN)))^ρ),0];
q1Atil=If[eC>0,Min[1,uinv[u[q0A]+(A/(eC M))(eC q0A+(1-eC)q0B)]],1];
q1Btil=If[eC<1,Min[1,uinv[u[q0B]+(B/((1-eC) M))(eC q0A+(1-eC)q0B)]],1];
ListPlot[{{δB,100(fA(λ v[q1A]+(1-λ)v[q1Atil])+fB(λ v[q1B]+(1-λ)v[q1Btil])+(ℓ eC-fA)v[q0A]+(ℓ(1-eC)-fB)v[q0B])/ℓ}/.semistrategicresultsA2irs1},
PlotRange->{{.6,1},{-1.33,-1.225}},AxesOrigin->{.6,-1.33},AspectRatio->1,ImageSize->260,Joined->True,AxesLabel->{" (Subscript[δ, B])","(%)"},BaseStyle->{FontSize->18},PlotStyle->{{Thick,Darker@Green}}]]
Export["csss-irs-asmall-welfare.pdf",%]