(* 

Show P is diagonal.

Suppose equations 1, 2 and 3 hold.

*)
Get["basefile.m"];

ClearMonomialOrder[];

SetKnowns[Join[
Bknowns,
BknownsTP,
Cknowns,
CknownsTP,
Dvariables,
DvariablesTP,
Avariables,
AvariablesTP
]];

SetUnknowns[Join[othervariables, othervariablesTP],
AllPiiVars,
AllOffPVars
];

AllFiveEquations = Complement[AllFiveEquations,{True}];

polys= Union[AllFiveEquations, ASSUMPTIONS, matPtilSelfAdjoint,
             matDtilSPECIFIC, PiiInvertible, greekEquations ];
polys = polys/.Equal->Subtract;
polys = ExpandNonCommutativeMultiply[polys];

Banswer1 = NCProcess[polys,3,0,0,0,"Banswer1",
                 RR->False,
                 SBByCat->False,
                 DegreeCap->4,DegreeSumCap->7,NCCV->False,
                 AdditionalRegularOutputOptions->
                 {TeXReplace->{"Sigma"->"\\Sigma"},CategoryCutoff->5}];
