t = 0.3 c1 = 0.1 c2 = 0.07 p1 = 92.21 p2 = 75.84 n = 10 a11 = 100*c1 - p1/n a12 = 100*c2 - p2/n a21 = 100 - (100 - p1) * t + 100 * c1 * (1-t) - t*p1/n a22 = 100 - (100 - p2) * t + 100 * c2 * (1-t) - t*p2/n A = matrix(data=c(a11, a12, a21, a22), nrow=2, ncol=2, byrow=TRUE) #print(A) b = c( -100/n, 100 - 100*t/n ) #print(b) x = solve(A, b) print(x) p0 = sum( x * c(p1, p2) ) print(p0)