Sistemas de ecuaciones lineales en R

En R se pueden resolver fácilmente sistemas de ecuaciones lineales con la función solve.

Ejemplo

2x + 3y + 3z = 20

x + 4y + 3z = 15

5x + 3y + 4z = 30

Solución

a <- rbind(c(2, 3, 3), 
           c(1, 4, 3), 
           c(5, 3, 4))
b <- c(20, 15, 30)
solve(a, b)

[1] -1.25 -6.25 13.75

Si queremos el resultado en fracciones, utilizamos la función fractions del paquete MASS.

library(MASS)
fractions(solve(a, b))

[1]  -5/4 -25/4  55/4