python3 <<'PY' import math a=0.25; b=1/12 P=[[a,0,b],[b,0,a],[0,b,a]] pX=[sum(P[x]) for x in range(3)] pY=[sum(P[x][y] for x in range(3)) for y in range(3)] print(pX,pY) aff_Y=sum(pY[y]**1.5 for y in range(3)) aff_X=sum(pX[x]**1.5 for x in range(3)) print(aff_Y,1-aff_Y,aff_X,1-aff_X) print("done") PY