python3 <<'EOF' V=[[6.54,1.63,4.71,6.64,5.51,3.29,6.94,3.11],[6.05,0.70,1.41,4.56,4.18,5.37,5.30,2.32],[6.76,3.83,1.47,6.65,3.21,6.57,0.97,3.16]] p=[1.0,0.88,0.73,0.61] rf=[0.98,0.96,0.94,0.92] pr=[0.25,0.45,0.30] pg=[0.75,0.40,0.15] A=[3,2,4,1] B=[1,3,1,2] g=[0,3,6,4] r=[0,3,5,4] for s in range(3): wg=pr[s]*pg[s] wr=pr[s]*(1-pg[s]) lg=0 lr=0 for k in range(4): if g[k]==r[k]: lg+=wg*p[k]*V[s][g[k]] lr+=wr*p[k]*V[s][r[k]] else: lg+=wg*p[k]*rf[k]*V[s][g[k]] lr+=wr*p[k]*rf[k]*V[s][r[k]] ptg=[1]*8 ptr=[1]*8 for k in range(4): if g[k]!=r[k]: ptg[g[k]]+=A[k] ptg[r[k]]+=B[k] ptr[r[k]]+=A[k] ptr[g[k]]+=B[k] for k in range(4): if g[k]==r[k]: ptg[g[k]]=0 ptr[r[k]]=0 print(s,lg,lr,ptg,ptr,wg*0.4*sum(ptg[j]*V[s][j] for j in range(8))/sum(ptg),wr*0.4*sum(ptr[j]*V[s][j] for j in range(8))/sum(ptr)) EOF