row1 <- c(2,-3,1,7) row2 <- c(2,8,-4,5) row3 <- c(1,3,-3,0) row4 <- c(-5,2,3,4) matrx <- rbind(row1, row2, row3, row4) solutions_vector <- c(14,-1,4,19) solve(matrx, solutions_vector)
## [1] -8.043841 -4.665971 -8.680585 3.538622