\[ x + y = 5 \\ 2x -y = 3 \]
# Define the coefficients matrix (A) and the right-hand side vector (b)
# Solve the system of equations
solution <- solve(matrix(c(1,1,2,-1), nrow = 2, byrow = TRUE), c(5,3))
# Print the solution
print(solution)## [1] 2.666667 2.333333