2015-02-13

Sistemas de ecuaciones lineales en R

Title 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

No hay comentarios:

Publicar un comentario

Nube de datos