(* Computing the Post-Merger Equilibrium *)

<<Merger`



merging = {1, 2};

mergedProfits = Map[Plus, Part[profits, merging]]



rules = Flatten[{

      Table[a[i] -> aa[[i]], {i, n}],

      Table[c[i] -> cc[[i]], {i, n}], p[5] -> 0}]

eqns = Table[If[

        MemberQ[merging, i],

        D[mergedProfits, p[i]],

        D[profits[[i]], p[i]]], {i, n - 1}] /. rules



startingValues = Table[{p[i], p0[[i]]}, {i, 1, n - 1}]

FindRoot @@ FlattenAt[{eqns, startingValues, MaxIterations -> 500}, {2}]

NumberForm[

  TableForm[Transpose[

      {{"firm", FirmNames}, {"P0", p0}, {"Q0", q00},

        {"P1", p1 = #[[2]] & @ PostMerger[merging = {1, 2}, aa, b, cc, p0]},

        {"Q1", q1 = #[[3]] & @ PostMerger[merging, aa, b, cc, p0]}}]], 3]