En R
base la función solve
permite obtener la inversa de una matriz cuadrada (filas = columnas)
la cual usa la forma Ax=b
. Con ello se puede resolver sistemas de ecuaciones lineales mediante la inversa de A

La implementación de las funciones en R
así como la asignación de valores en las variables requieren del operador <-
en lugar del =
Opción 1:
linear_solver_1 <- function(A,b){
# Ax + b = 0
solve(A) %*% (-b) # solve(A): la inversa de la matriz cuadrada A
}
linear_solver_2 <- function(A,b){
# Ax = b
solve(A) %*% (b) # solve(A): la inversa de la matriz cuadrada A
}
Su empleo:
# 2x + 3 = 0
print(linear_solver_1(2,3)) # Rpta: -1.5
# 2x = -3
print(linear_solver_3(2,-3)) # Rpta: -1.5
# x + y - 4 = 0
# x - y - 2 = 0
a <- rbind(c(1,1), c(1,-1))
# [,1] [,2]
# [1,] 1 1
# [2,] 1 -1
b <- c(-4,-2)
print(linear_solver_1(a,b)) # Rpta: 3, 1 (x, y)
# x + y = 4
# x - y = 2
a <- rbind(c(1,1), c(1,-1))
b <- c(4,2)
print(linear_solver_2(a,b)) # Rpta: 3, 1 (x, y)
Opción 2:
También si se tiene que a %*% x = b
basta con hacer lo siguiente solve(a, b)
linear_solver_1 <- function(A,b){
# Ax + b = 0
solve(A, -b)
}
linear_solver_2 <- function(A,b){
# Ax = b
solve(A, b)
}
Qué es %*%
Este es el operador de multiplicación para matrices que emplea R
.