(* Demand Compatible and Ramsey Pricing *)



<<Pricing`

demandsystem1 = CESDemandSystem[{1, 1}, {r, r}, r, y]

demandsystem2 = 

  LinearDemandSystem[1000, {10, 20, 30, 40}, 

    1.1 IdentityMatrix[4] - Table[.1, {4}, {4}]]

Demands[{p1, p2, p3, p4}, demandsystem2]

Demands[{1, 2, 3, 4}, demandsystem2]

demandsystem3 = 

  LinearDemandSystem[100, {10, 20, 30}, 

    1.1 IdentityMatrix[3] - Table[.1, {3}, {3}]]

costfunction6 = QuadraticCostFunction[10, {1, 1, 1}, Table[0, {3}, {3}]]



FP[z_] := N[ShapleyValuePrice[Demands[z, demandsystem3],

          costfunction6] /. x_ /; x < 0 -> 0];

svp = {1, 1, 1}

Do[svp = FP[svp]; Print[svp], {5}]



FixedPoint[FP, {1, 1, 1}]

costfunction7 = CESCostFunction[{1, 1, 1, 1}, {2, 2, 2, 2}, 3]



FP[z_] := AumannShapleyPrice[Demands[z, demandsystem2],

        costfunction7] /. x_ /; x < 0 -> 0;

svp = {1, 1, 1, 1}

Do[asp = FP[asp]; Print[asp], {5}]



FixedPoint[FP, {1, 1, 1, 1}]