A = matrix(c(1, -1, 1, 1), nrow=2, ncol=2, byrow=T) ## Solve the linear systems one at a time: ## x1 = solve(A, c(2*5*11, 4*1*1)) x2 = solve(A, c(2*1*11, 4*5*1)) x3 = solve(A, c(2*5*1, 4*1*11)) x4 = solve(A, c(2*1*1, 4*5*11)) x5 = solve(A, c(4*5*11, 2*1*1)) x6 = solve(A, c(4*1*11, 2*5*1)) x7 = solve(A, c(4*5*1, 2*1*11)) x8 = solve(A, c(4*1*1, 2*5*11)) R = cbind(x1, x2, x3, x4, x5, x6, x7, x8) rownames(R) = c('n', 'b') print(R)