Defining the vectors

b <- matrix(c(-1, 0, 4), nrow = 3, byrow = TRUE) c <- matrix(c(3, 4, 2), nrow = 3, byrow = TRUE) a_coefficient <- matrix(c(1, 2, -1), nrow = 3, byrow = TRUE)

Calculate the right-hand side using matrix operations

rhs <- b - 4 * c

Solve for ‘a’ using element-wise division

a <- rhs / a_coefficient

Print the solution

print(“Solution for a:”) print(a)