En este artículo, discutiremos cómo resolver un sistema de ecuaciones en el lenguaje de programación R.
La función solve() en lenguaje R se usa para resolver la ecuación. Aquí la ecuación es como a*x = b, donde b es un vector o array yx es una variable cuyo valor se va a calcular.
Sintaxis: solve(a, b)
Parámetros:
- a: coeficientes de la ecuación
- b: vector o array de la ecuación
Ejemplo 1: Resolución de la ecuación del sistema de tres ecuaciones
Given Equations: x + 2y + 3z = 20 2x + 2y + 3z = 100 3x + 2y + 8z = 200 Matrix A and B for solution using coefficient of equations: A-> 1 2 3 2 2 3 3 2 8 B-> 20 100 200
Para resolver esto usando dos arrays en R usamos el siguiente código:
R
# create matrix A and B using given equations A <- rbind(c(1, 2, 3), c(2, 2, 3), c(3, 2, 8)) B <- c(20, 100, 200) # Solve them using solve function in R solve(A, B)
Producción:
80 -36 3.99999999999999
lo que significa que x=80, y=-36 y z=4 es la solución para ecuaciones lineales.
Ejemplo 2: Resolución de la ecuación del sistema de tres ecuaciones
Para obtener soluciones en forma de fracciones, usamos la biblioteca MASS en lenguaje R y ajustamos la función de resolución en fracciones.
Given Equations: 19x + 32y + 31z = 1110 22x + 28y + 13z = 1406 31x + 12y + 81z = 3040 Matrix A and B for solution using coefficient of equations: A-> 19 32 31 22 28 13 31 12 81 B-> 1110 1406 3040
Para resolver esto usando dos arrays en R usamos el siguiente código:
R
# Load package MASS library(MASS) # create matrix A and B using given equations A <- rbind(c(19, 32, 31), c(22, 28, 31), c(31, 12, 81)) B <- c(1110, 1406, 3040) # Solve them using solve # function wrapped in fractions fractions(solve(A, B))
Producción:
[1] 159950/2243 -92039/4486 29784/2243
lo que significa que x=159950/2243, y=-92039/4486 y z=29784/2243 es la solución para la ecuación lineal anterior.
Ejemplo 3: Resolución de array inversa
R
# create matrix A and B using given equations
A <- matrix(c(4, 7, 3, 6), ncol = 2)
print(A)
print("Inverse matrix")
# Solve them using solve function in R
print(solve(A))
Producción:
[,1] [,2]
[1,] 4 3
[2,] 7 6
[1] "Inverse matrix"
[,1] [,2]
[1,] 2.000000 -1.000000
[2,] -2.333333 1.333333
Publicación traducida automáticamente
Artículo escrito por mishrapriyank17 y traducido por Barcelona Geeks. The original can be accessed here. Licence: CCBY-SA