(* Logit Choice Probabilities *)

<<Merger`

<mergerSimulationFunction`



n = 5;

rules = Flatten[{a[n] -> 10, Table[p[i] -> p0[[i]], {i, n}]}]

eqns = Take[q00 - q /. rules, {1, n - 1}]



startingValues = Table[{a[i], 20}, {i, 1, n - 1}]

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

NumberForm[TableForm[Transpose[

      {{"firms", FirmNames}, {"p0", p0}, {"choice prob.", 

          q00}, {"demand constants", 

          aa = DemandConstants[{}, 1., 9.557, p0, q00]}}]], 3]