r_f = 0.1
b = c(0.15, 0.2) - r_f
A = matrix(data=c(2, 1, 3, 4), nrow=2, ncol=2, byrow=TRUE)
print(solve(A, b))