METODOS NUMERICOS

Unidad 2
Tarea 2 - Ecuaciones Lineales e Interpolación

APORTE INDIVIDUAL

Presentado a:
Ing. JAIRO LUIS GUTIERREZ
Tutor

Entregado por:

Fabio Hernández Mancipe

Código: 80090549

Grupo
100401_4

UNIVERSIDAD NACIONAL ABIERTA Y A DISTANCIA (UNAD).

ESCUELA DE CIENCIAS BÁSICAS TECNOLOGÍA E INGENIERA
Julio2019

pág. 1
 INTRODUCCION

En esta unidad veremos el sistema de ecuaciones lineales de algebra lineal (otra

vez), pero visto de otras técnicas; veremos los métodos iterativos para su
resolución, ventajas y desventajas.

pág. 2
 DESARROLLO DE EJERCICIOS

Aporte 1. Solución de Sistemas de Ecuaciones Lineales (SEL)

Teniendo en cuenta que mi documento termina en 9, me corresponde el

ejercicio impar.

Dado el siguiente sistema de ecuaciones lineales (SEL):

Para estudiantes cuyo documento de identidad termine en dígito impar:

3𝑥1 − 13𝑥2 + 9𝑥3 + 3𝑥4 = −19

−6𝑥1 + 4𝑥2 + 𝑥3 − 18𝑥4 = −34
6𝑥1 − 2𝑥2 + 2𝑥3 + 4𝑥4 = 16
12𝑥1 − 8𝑥2 + 6𝑥3 + 10𝑥4 = 26

Para cada uno de los estudiantes:

Según su último número de cédula (par o impar) y para el SEL que le
corresponde, realice los siguientes tres ejercicios:

Ejercicio 1.
Determine si el SEL tiene o no solución. ¿Si la tiene, es única? Realice una
breve explicación basándose en la teoría revisada.

Para realizar el proceso de verificación, usare la regla de cramer para hallar saber si
el sistema tiene solución.

“Para determinar si el sistema tiene solución única o no, debemos calcular su

determinante. Si este nos da diferente de cero (0), entonces el sistema tendrá única
solución y además la inversa de la matriz de coeficientes existirá (y esta será
única)” (Zuñiga, 2010, p.204)

pág. 3
Tenemos la matriz A

3 −13 9 3
−6 4 1 −18
𝐴=| |
6 −2 2 4
12 −8 6 10

Los términos independientes:

−19
−34
𝐵=| |
16
26

Las incógnitas:

𝑥1
𝑥2
𝑋 = |𝑥 |
3
𝑥4

Tomamos la matriz A y la vamos repetir para complementar hacia la derecha para

poder tomar todas las líneas superiores e inferiores en diagonal.

3 −13 9 3 3 −13 9
−6 4 1 −18 −6 4 1
𝐴=| || |
6 −2 2 4 6 −2 2
12 −8 6 10 12 −8 6

Ojo…. TOMANMOS EN DIAGONAL…. ARRANCAMOS CON LOS SUPERIORES….

REPETIMOS CON LOS INFERIORES…. LUEGO SUPERIORES MENOS INFERIORES =
RESULTADO = DETERMINANTE

|𝐴| = ((3)(4)(2)(10) + (−13)(1)(4)(12) + (9)(−18)(6)(−8) + (3)(−6)(−2)(6))

−((12)(−2)(1)(3) + (−8)(2)(−18)(3) + (6)(4)(−6)(−13) + (10)(6)(4)(9))

|𝐴| = ((240) + (−624) + (7776) + (216)) − ((−72) + (864) + (1872) + (2160))

|𝐴| = (240 − 624 + 7776 + 216) − (−72 + 864 + 1872 + 2160)

|𝐴| = (7608) − (4824)

|𝐴| = 7608 − 4824 = 2784

|𝑨| = 𝟐𝟕𝟖𝟒

Teniendo en cuenta que la determinante es diferente de cero, esta tendrá única

solución, por lo tanto, es NO SINGULAR.

pág. 4
(Nieves, 2014, P. 184)

pág. 5
Ejercicio 2.
Resuelva el SEL por los esquemas de Eliminación Gaussiana simple y Gauss
- Jordan. ¿Es necesario emplear la estrategia de Pivoteo? ¿Por qué? Valide
que las respuestas determinadas son solución del SEL.

3𝑥1 − 13𝑥2 + 9𝑥3 + 3𝑥4 = −19

−6𝑥1 + 4𝑥2 + 𝑥3 − 18𝑥4 = −34
6𝑥1 − 2𝑥2 + 2𝑥3 + 4𝑥4 = 16
12𝑥1 − 8𝑥2 + 6𝑥3 + 10𝑥4 = 26

Eliminación Gaussiana Simple

Matriz Ampliada

3 −13 9 3 −19
−6 4 1 −18 −34
| || |
6 −2 2 4 16
12 −8 6 10 26

𝑓4 − 4𝑓1

3 −13 9 3 −19
−6 4 1 −18 −34
| || |
6 −2 2 4 16
0 44 −30 −2 102

𝑓3 − 2𝑓1

3 −13 9 3 −19
−6 4 1 −18 −34
| || |
0 24 −16 −2 54
0 44 −30 −2 102

𝑓2 + 2𝑓1

3 −13 9 3 −19
0 −22 19 −12 −72
| || |
0 24 −16 −2 54
0 44 −30 −2 102

𝑓4 + 2𝑓2

pág. 6
 3 −13 9 3 −19
0 −22 19 −12 −72
| || |
0 24 −16 −2 54
0 0 8 −26 −42

𝑓3 + 1.09090909𝑓2

3 −13 9 3 −19
0 −22 19 −12 −72
| || |
0 0 4.72727272 −15.09090909 −24.54545455
0 0 8 −26 −42

𝑓4 − 1.692307692𝑓3

3 −13 9 3 −19
0 −22 19 −12 −72
| || |
0 0 4.72727272 −15.09090909 −24.54545455
0 0 0 −0.461538461 −0.461538461

Teniendo el ultimo resultado, obtenemos 𝑥4

−0.461538461𝑥4 = −0.461538461

−0.461538461
𝑥4 =
−0.461538461

𝒙𝟒 = 𝟏

Ahora podemos hallar 𝑥3

4.72727272𝑥3 − 15.09090909𝑥4 = −24.54545455

4.72727272𝑥3 − 15.09090909(1) = −24.54545455

4.72727272𝑥3 = −24.54545455 + 15.09090909

4.72727272𝑥3 = −9.45454546

−9.45454546
𝑥3 =
4.72727272

𝒙𝟑 = −𝟐

pág. 7
Ahora podemos hallar 𝑥2

−22𝑥2 + 19𝑥3 − 12𝑥4 = −72

−22𝑥2 + 19(−2) − 12(1) = −72

−22𝑥2 − 38 − 12 = −72

−22𝑥2 = −72 + 38 + 12

−22𝑥2 = −22

−22
𝑥2 =
−22

𝒙𝟐 = 𝟏

Ahora podemos hallar 𝑥1

3𝑥1 − 13𝑥2 + 9𝑥3 + 3𝑥4 = −19

3𝑥1 − 13(1) + 9(−2) + 3(1) = −19

3𝑥1 − 13 − 18 + 3 = −19

3𝑥1 − 28 = −19

3𝑥1 = −19 + 28

3𝑥1 = 9

9
𝑥1 =
3

𝒙𝟏 = 𝟑

Resultados:
𝒙𝟏 = 𝟑
𝒙𝟐 = 𝟏
𝒙𝟑 = −𝟐
𝒙𝟒 = 𝟏

pág. 8
Eliminación Gauss-Jordan

La idea es reducir el SEL a la forma:

1 0 0 0
0 1 0 0
0 0 1 0
0 0 0 1

Para luego obtener sus valores de B.

3𝑥1 − 13𝑥2 + 9𝑥3 + 3𝑥4 = −19

−6𝑥1 + 4𝑥2 + 𝑥3 − 18𝑥4 = −34
6𝑥1 − 2𝑥2 + 2𝑥3 + 4𝑥4 = 16
12𝑥1 − 8𝑥2 + 6𝑥3 + 10𝑥4 = 26

Matriz Ampliada

3 −13 9 3 −19
−6 4 1 −18 −34
| || |
6 −2 2 4 16
12 −8 6 10 26

0.33333333𝑓1

1 −4.33333333 3 1 −6.33333333
−6 4 1 −18 −34
| || |
6 −2 2 4 16
12 −8 6 10 26

𝑓2 + 6𝑓1

1 −4.33333333 3 1 −6.33333333
0 −22 19 −12 −72
| || |
6 −2 2 4 16
12 −8 6 10 26

𝑓3 − 6𝑓1

1 −4.33333333 3 1 −6.33333333
0 −22 19 −12 −72
| || |
0 24 −16 −2 54
12 −8 6 10 26

pág. 9
𝑓4 − 12𝑓1

1 −4.33333333 3 1 −6.33333333
0 −22 19 −12 −72
| || |
0 24 −16 −2 54
0 44 −30 −2 102

−0.045454545𝑓2

1 −4.33333333 3 1 −6.33333333
0 1 −1.58333333 0.545454545 3.27272727
| || |
0 24 −16 −2 54
0 44 −30 −2 102

𝑓1 + 4.33333333𝑓2

1 0 −0.74242424 3.363636364 7.84848484

0 1 −1.58333333 0.545454545 3.27272727
| || |
0 24 −16 −2 54
0 44 −30 −2 102

𝑓3 − 24𝑓2

1 0 −0.74242424 3.363636364 7.84848484

0 1 −1.58333333 0.545454545 3.27272727
| || |
0 0 4.727272727 −15.09090909 −24.54545455
0 44 −30 −2 102

𝑓4 − 44𝑓2

1 0 −0.74242424 3.363636364 7.84848484

0 1 −1.58333333 0.545454545 3.27272727
| || |
0 0 4.727272727 −15.09090909 −24.54545455
0 0 8 −26 −42

0.211538461𝑓3

1 0 −0.74242424 3.363636364 7.84848484

0 1 −1.58333333 0.545454545 3.27272727
| || |
0 0 1 −3.192307692 −5.192307692
0 0 8 −26 −42

𝑓1 + 0.74242424𝑓3

1 0 0 0.993589743 3.993589744
0 1 −1.58333333 0.545454545 3.27272727
| || |
0 0 1 −3.192307692 −5.192307692
0 0 8 −26 −42

pág. 10
𝑓2 + 1.58333333𝑓3

1 0 0 0.993589743 3.993589744
0 1 0 −2.211538462 −1.211538462
| || |
0 0 1 −3.192307692 −5.192307692
0 0 8 −26 −42

𝑓4 − 8𝑓3

1 0 0 0.993589743 3.993589744
0 1 0 −2.211538462 −1.211538462
| || |
0 0 1 −3.192307692 −5.192307692
0 0 0 −0.461538461 −0.461538461

−2.166666667𝑓4

1 0 0 0.993589743 3.993589744
0 1 0 −2.211538462 −1.211538462
| || |
0 0 1 −3.192307692 −5.192307692
0 0 0 1 1

𝑓1 − 0.993589743𝑓4

1 0 0 0 3
0 1 0 −2.211538462 −1.211538462
| || |
0 0 1 −3.192307692 −5.192307692
0 0 0 1 1

𝑓2 + 2.211538462𝑓4

1 0 0 0 3
0 1 0 0 1
| || |
0 0 1 −3.192307692 −5.192307692
0 0 0 1 1

𝑓3 + 3.192307692𝑓4

1 0 0 0 3
0 1 0 0 1
| || |
0 0 1 0 −2
0 0 0 1 1

Teniendo la matriz de forma escalonada reducida, tenemos los resultados:

Resultados:
𝒙𝟏 = 𝟑
𝒙𝟐 = 𝟏
𝒙𝟑 = −𝟐
𝒙𝟒 = 𝟏

pág. 11
Validación y confirmación que respuestas del SEL son la solución:

𝑥1 = 3
𝑥2 = 1
𝑥3 = −2
𝑥4 = 1

3𝑥1 − 13𝑥2 + 9𝑥3 + 3𝑥4 = −19

3(3) − 13(1) + 9(−2) + 3(1) = −19
9 − 13 − 18 + 3 = −19
−𝟏𝟗 = −𝟏𝟗

−6𝑥1 + 4𝑥2 + 𝑥3 − 18𝑥4 = −34

−6(3) + 4(1) + (−2) − 18(1) = −34
−18 + 4 − 2 − 18 = −34
−𝟑𝟒 = −𝟑𝟒

6𝑥1 − 2𝑥2 + 2𝑥3 + 4𝑥4 = 16

6(3) − 2(1) + 2(−2) + 4(1) = 16
18 − 2 − 4 + 4 = 16
𝟏𝟔 = 𝟏𝟔

12𝑥1 − 8𝑥2 + 6𝑥3 + 10𝑥4 = 26

12(3) − 8(1) + 6(−2) + 10(1) = 26
36 − 8 − 12 + 10 = 26
𝟐𝟔 = 𝟐𝟔

¿Es necesario emplear la estrategia de Pivoteo? ¿Por qué?

“El coeficiente de la incógnita que se eliminará en la fila pivote, se llama pivote…”

(Nieves, 2014, p. 190)

Si es necesario.
La estrategia de pivoteo fue utilizado para el método gaussiano, ya que el método
escalonado requiere establecer el punto de referencia que multiplica con la incógnita
para llevar el valor de fila a cero, es decir, la matriz llevada a su forma escalonada
necesita el pivote que multiplica con incógnita.

pág. 12
Ejercicio 3.
OBSERVACIONES GENERALES:

k está dado por la siguiente tabla:

Último dígito cédula k

1o6 5
2o7 6
3o8 7
4o9 8
0o5 9

Teniendo en cuenta que mi documento termina en 9, me corresponde el valor

k=8.

Resuelva el SEL por los esquemas iterativos: Gauss – Seidel, Jacobi y S.O.R.
Debe establecer los criterios de convergencia en cada método.
Itere la solución hasta que ésta converja a una tolerancia de
‖𝑨‖ ≤ 𝟏𝟎−𝒌
Donde ‖𝑨‖ es una norma matricial para la matriz A (escoja la que
considere más adecuada) y k está dada por la tabla que se presentará más
abajo en esta guía. Valide que las respuestas determinadas son solución
del SEL.

Realice una gráfica de la forma como va convergiendo la solución (Número

de iteraciones vs norma del error). Realice una breve explicación,
sustentándose en la teoría revisada, acerca de los resultados. ¿Cuál
considera que es el mejor método para este SEL en particular? ¿por qué?

Gauss – Seidel:
“El método de Gauss-Seidel es particularmente adecuado cuando se tiene gran
número de ecuaciones. En estos casos, los métodos de eliminación pueden estar
sujetos a errores de redondeo. Debido a que el error en el método de Gauss-Seidel
es determinado por el número de iteraciones, el error de redondeo no es un tema
que preocupe a este método. Aunque, existen ciertos ejemplos donde la técnica de
Gauss-Seidel no convergerá al resultado correcto”. (Chapra, 2007, p. 305)

pág. 13
Una de las reglas para dar inicio al método, Esque los elementos de la diagonal NO
deben ser cero.

Si observamos detalladamente, ninguno es cero en diagonal de a11 a a44(x1=3,

x2=4, x3=2, x4=10)

3𝑥1 − 13𝑥2 + 9𝑥3 + 3𝑥4 = −19

−6𝑥1 + 4𝑥2 + 𝑥3 − 18𝑥4 = −34
6𝑥1 − 2𝑥2 + 2𝑥3 + 4𝑥4 = 16
12𝑥1 − 8𝑥2 + 6𝑥3 + 10𝑥4 = 26

Estableciendo en matriz ampliada:

3 −13 9 3 −19
−6 4 1 −18 −34
| || |
6 −2 2 4 16
12 −8 6 10 26

La zona en amarillo es la matriz 𝑎, y la zona en azul es el término independiente 𝑏.

Recuerden que para interpretar el posicionamiento en una matriz, el 1er número

corresponde a la fila y el segundo número corresponde a la columna; ejemplo:

𝑎32 = 𝑀𝑎𝑡𝑟𝑖𝑧 𝑎, 𝐹𝑖𝑙𝑎 3, 𝑐𝑜𝑙𝑢𝑚𝑛𝑎 2 = −2

𝑎13 = 𝑀𝑎𝑡𝑟𝑖𝑧 𝑎, 𝐹𝑖𝑙𝑎 1, 𝑐𝑜𝑙𝑢𝑚𝑛𝑎 3 = 9
𝑏3 = 𝑇𝑒𝑟𝑚𝑖𝑛𝑜 𝑖𝑛𝑑𝑒𝑝𝑒𝑛𝑑𝑖𝑒𝑛𝑡𝑒 𝑏, 𝐹𝑖𝑙𝑎 3 = 16

Las ecuaciones utilizadas para para un sistema SEL de 4x4, es de:

𝑏1 − 𝑎12 𝑥2 − 𝑎13 𝑥3 − 𝑎14 𝑥4

𝑥1 =
𝑎11

𝑏2 − 𝑎21 𝑥1 − 𝑎23 𝑥3 − 𝑎24 𝑥4

𝑥2 =
𝑎22

𝑏3 − 𝑎31 𝑥1 − 𝑎32 𝑥2 − 𝑎34 𝑥4

𝑥3 =
𝑎33

𝑏4 − 𝑎41 𝑥1 − 𝑎42 𝑥2 − 𝑎43 𝑥3

𝑥4 =
𝑎44

Observarán que las ecuaciones están hechas para que los valores ocultos tomen un
valor de resultado; la idea es que inmediatamente sea identificado este valor, ese

pág. 14
mismo valor sea utilizado en la siguiente ecuación; teniendo en cuenta que no
tenemos valores iniciales de x, iniciaremos con Cero en todos los valores ocultos:

𝑥1 = 0
𝑥2 = 0
𝑥3 = 0
𝑥4 = 0

Iniciamos la búsqueda de los valores ocultos con las formulas; usaremos el número
de iteraciones que sean necesarios, hasta que converjan.

“Conforme un nuevo valor de x se calcula con el método de Gauss-Seidel, éste se

usa inmediatamente en la siguiente ecuación para determinar el otro valor de x. De
esta forma, si la solución es convergente, se empleará la mejor aproximación
disponible”.
(Chapra, 2007, p. 312)

Nota: He realizado la ecuación pero esto tiene una trampa ya que diverge, por
lo tanto, hay que utilizar una jugada de rearreglo para lograr convergencia.

“Si el proceso iterativo diverge, como es el caso, un rearreglo de las ecuaciones

puede originar convergencia; por ejemplo, en lugar de despejar x1 de la primera
ecuación, x2 de la segunda, etc., cabe despejar las diferentes xi de diferentes
ecuaciones, teniendo cuidado de que los coeficientes de las xi despejadas sean
distintos de cero.”
(Nieves, 2014, p. 240)

Realizando ejercicio en orden de:

6𝑥1 − 2𝑥2 + 2𝑥3 + 4𝑥4 = 16
12𝑥1 − 8𝑥2 + 6𝑥3 + 10𝑥4 = 26
3𝑥1 − 13𝑥2 + 9𝑥3 + 3𝑥4 = −19
−6𝑥1 + 4𝑥2 + 𝑥3 − 18𝑥4 = −34

pág. 15
Estableciendo en matriz ampliada:

6 −2 2 4 16
12 −8 6 10 26
| || |
3 −13 9 3 −19
−6 4 1 −18 −34

𝑥1 = 0
𝑥2 = 0
𝑥3 = 0
𝑥4 = 0

16 − (−2𝑥2 ) − (2𝑥3 ) − (4𝑥4 )

𝑥1 =
6

26 − (12𝑥1 ) − (6𝑥3 ) − (10𝑥4 )

𝑥2 =
−8

−19 − (3𝑥1 ) − (−13𝑥2 ) − (3𝑥4 )

𝑥3 =
9

−34 − (−6𝑥1 ) − (4𝑥2 ) − (1𝑥3 )

𝑥4 =
−18

Iteración 1:

16 − (−2 ∗ 0) − (2 ∗ 0) − (4 ∗ 0) 16 + 0 − 0 − 0
𝑥1 = = = 2.66666666
6 6

Ojo.. ya tenemos un valor de x1, por lo tanto debemos usarlo para la siguiente
ecuación… así vamos haciendo en todo.. vamos buscando resultados de x y los
vamos usando inmediatamente.

26 − (12 ∗ 2.66666666) − (6 ∗ 0) − (10 ∗ 0) −5.99999992 − 0 − 0

𝑥2 = = = 0.74999999
−8 −8

−19 − (3 ∗ 2.66666666) − (−13 ∗ 0.74999999) − (3 ∗ 0)

𝑥3 =
9
−26.99999998 + 9.74999987 − 0
= = −1.9166666788
9

−34 − (−6 ∗ 2.66666666) − (4 ∗ 0.74999999) − (1 ∗ −1.916666788)

𝑥4 =
−18
−18 − 3 + 1.916666788
= = 1.06018517844
−18

pág. 16
Hasta este punto, tenemos los valores de x:

𝑥1 = 2.66666666
𝑥2 = 0.74999999
𝑥3 = −1.9166666788
𝑥4 = 1.06018517844

Seguimos con la siguiente iteración y repetimos.

Iteración 2:

16 − (−2 ∗ 0.74999999) − (2 ∗ −1.9166666788) − (4 ∗ 1.06018517844)

𝑥1 =
6
17.49999998 + 308333333576 − 4.24074071376
𝑥1 = = 2.84876543731
6

Recuerden que los nuevos resultados de x se usan inmediatamente.

26 − (12 ∗ 2.84876543731) − (6 ∗ −1.9166666788) − (10 ∗ 1.06018517844)

𝑥2 =
−8
−8.1851852477 + 11.5 − 10.6018517844
𝑥2 = = 0.91087961991
−8

−19 − (3 ∗ 2.84876543731) − (−13 ∗ 0.91087961991) − (3 ∗ 1.06018517844)

𝑥3 =
9
−27.5462963119 + 11.8414350588 − 3.18055553532
𝑥3 = = −2.09837964315
9

−34 − (−6 ∗ 2.84876543731) − (4 ∗ 0.91087961991) − (1 ∗ −2.09837964315)

𝑥4 =
−18
−16.9074073762 − 3.64351847964 + 2.09837964315
𝑥4 = = 1.02514145626
−18

Hasta este punto, tenemos los valores de x:

𝑥1 = 2.84876543731
𝑥2 = 0.91087961991
𝑥3 = −2.09837964315
𝑥4 = 1.02514145626

Cada vez que se termine una iteración, debe verificarse el error porcentual para
validar si los 4 valores han convergido, esto en función de ‖𝐴‖ ≤ 10−8 ; si no han
convergido, se debe continuar con las iteraciones.

pág. 17
Tabla de resultados Gauss – Seidel:

|(𝑥)𝑎𝑐𝑡𝑢𝑎𝑙 − (𝑥)𝑎𝑛𝑡𝑒𝑟|
𝐸𝑟𝑟𝑜𝑟% = 𝐸% = ∗ 100%
(𝑥)𝑎𝑐𝑡𝑢𝑎𝑙

𝑃𝑜𝑟𝑐𝑒𝑛𝑡𝑎𝑗𝑒 𝑑𝑒 𝑒𝑟𝑟𝑜𝑟 𝑒𝑛 𝑓𝑢𝑛𝑐𝑖𝑜𝑛 𝑐𝑜𝑛 𝑠𝑢 𝑙𝑖𝑚𝑖𝑡𝑒 𝑎 𝑐𝑒𝑟𝑜.

𝐸% > 𝑇𝑜𝑙𝑒𝑟𝑎𝑛𝑐𝑖𝑎 = 𝐶𝑜𝑛𝑣𝑒𝑟𝑔𝑒; 𝑑𝑒 𝑙𝑜 𝑐𝑜𝑛𝑡𝑟𝑎𝑟𝑖𝑜 𝑁𝑂 𝑐𝑜𝑛𝑣𝑒𝑟𝑔𝑒

𝑇𝑜𝑙𝑒𝑟𝑎𝑛𝑐𝑖𝑎 𝑠𝑒𝑔ú𝑛 𝐶é𝑑𝑢𝑙𝑎 = 10−8 = 0.00000001

Iter
ació Error( Error(% Error(% Error(%)x
n x1 x2 x3 x4 %)x1 )x2 )x3 4 Rx1 Rx2 Rx3 Rx4
0.0000 0.0000 0.0000 0.0000
0 0000 0000 0000 0000
-
2.6666 0.7500 1.9166 1.0601
1 6667 0000 6667 8519
Ya
- convergi
-
2.8487 0.9108 2.0983 1.0251 0.0639 0.176620 0.086596
NO HA NO HA ó.. Ya lo NO HA
convergido. convergido.. tiene convergido
2 6543 7963 7963 4146 2199 08 80 0.03418428 . continúe continúe pequeño .. continúe
Ya
- convergi
-
2.9863 0.9371 2.0946 0.9853 0.0460 0.028011 0.001788
NO HA NO HA ó.. Ya lo NO HA
convergido. convergido.. tiene convergido
3 2545 3027 3413 2968 6330 73 14 0.04040453 . continúe continúe pequeño .. continúe
Ya
- convergi
-
3.0203 0.9412 2.0867 0.9753 0.0112 0.004365 0.003765
NO HA NO HA ó.. Ya lo NO HA
convergido. convergido.. tiene convergido
4 6835 3902 7631 3165 7111 26 53 0.01025090 . continúe continúe pequeño .. continúe
Ya
- convergi
-
3.0257 0.9427 2.0830 0.9740 0.0017 0.001611 0.001786
NO HA NO HA ó.. Ya lo NO HA
convergido. convergido.. tiene convergido
5 8401 5835 5427 7084 8984 57 82 0.00129438 . continúe continúe pequeño .. continúe
Ya
- convergi
-
3.0258 0.9441 2.0806 0.9744 0.0000 0.001456 0.001139
NO HA NO HA ó.. Ya lo NO HA
convergido. convergido.. tiene convergido
6 9032 3331 8337 7267 3513 33 48 0.00041236 . continúe continúe pequeño .. continúe
Ya
- convergi
-
3.0252 0.9455 2.0786 0.9750 0.0001 0.001460 0.000991
NO HA NO HA ó.. Ya lo NO HA
convergido. convergido.. tiene convergido
7 9045 1398 2307 9390 9828 23 19 0.00063710 . continúe continúe pequeño .. continúe
Ya
- convergi
-
3.0246 0.9468 2.0766 0.9757 0.0002 0.001437 0.000949
NO HA NO HA ó.. Ya lo NO HA
convergido. convergido.. tiene convergido
8 4975 7470 5110 1940 1183 06 59 0.00064107 . continúe continúe pequeño .. continúe
Ya
- convergi
-
3.0240 0.9482 2.0747 0.9763 0.0002 0.001402 0.000925
NO HA NO HA ó.. Ya lo NO HA
convergido. convergido.. tiene convergido
9 2900 0442 3196 2843 0527 36 00 0.00062380 . continúe continúe pequeño .. continúe
Ya
- convergi
-
3.0234 0.9495 2.0728 0.9769 0.0001 0.001365 0.000902
NO HA NO HA ó.. Ya lo NO HA
convergido. convergido.. tiene convergido
10 2651 0133 6084 2141 9927 88 68 0.00060699 . continúe continúe pequeño .. continúe
Ya
- convergi
-
3.0228 0.9507 2.0710 0.9774 0.0001 0.001329 0.000880
NO HA NO HA ó.. Ya lo NO HA
convergido. convergido.. tiene convergido
11 3978 6581 3645 9934 9410 96 90 0.00059123 . continúe continúe pequeño .. continúe

etc… etc… etc… etc… etc… etc… etc… etc… etc… etc… etc… etc… etc…

pág. 18
etc… etc… etc… etc… etc… etc… etc… etc… etc… etc… etc… etc… etc…

etc… etc… etc… etc… etc… etc… etc… etc… etc… etc… etc… etc… etc…

etc… etc… etc… etc… etc… etc… etc… etc… etc… etc… etc… etc… etc…
Ya
- Ya convergi
-
3.0000 0.9999 2.0000 0.9999 0.0000 0.000000 0.000000
convergió.. NO HA ó.. Ya lo NO HA
Ya lo tiene convergido.. tiene convergido
423 0066 9857 0206 9935 0001 04 03 0.00000002 pequeño continúe pequeño .. continúe
Ya
- Ya convergi
-
3.0000 0.9999 2.0000 0.9999 0.0000 0.000000 0.000000
convergió.. NO HA ó.. Ya lo NO HA
Ya lo tiene convergido.. tiene convergido
424 0065 9861 0201 9936 0001 04 03 0.00000002 pequeño continúe pequeño .. continúe
Ya
- Ya convergi
-
3.0000 0.9999 2.0000 0.9999 0.0000 0.000000 0.000000
convergió.. NO HA ó.. Ya lo NO HA
Ya lo tiene convergido.. tiene convergido
425 0063 9864 0196 9938 0001 03 03 0.00000002 pequeño continúe pequeño .. continúe
Ya
- Ya convergi
-
3.0000 0.9999 2.0000 0.9999 0.0000 0.000000 0.000000
convergió.. NO HA ó.. Ya lo NO HA
Ya lo tiene convergido.. tiene convergido
426 0061 9868 0191 9940 0001 03 02 0.00000002 pequeño continúe pequeño .. continúe
Ya
- Ya convergi
-
3.0000 0.9999 2.0000 0.9999 0.0000 0.000000 0.000000
convergió.. NO HA ó.. Ya lo NO HA
Ya lo tiene convergido.. tiene convergido
427 0060 9871 0186 9941 0001 03 02 0.00000002 pequeño continúe pequeño .. continúe
Ya
- Ya convergi
-
3.0000 0.9999 2.0000 0.9999 0.0000 0.000000 0.000000
convergió.. NO HA ó.. Ya lo NO HA
Ya lo tiene convergido.. tiene convergido
428 0058 9874 0182 9942 0000 03 02 0.00000001 pequeño continúe pequeño .. continúe
Ya
- Ya convergi
-
3.0000 0.9999 2.0000 0.9999 0.0000 0.000000 0.000000
convergió.. NO HA ó.. Ya lo NO HA
Ya lo tiene convergido.. tiene convergido
429 0057 9877 0177 9944 0000 03 02 0.00000001 pequeño continúe pequeño .. continúe
Ya
- Ya convergi
-
3.0000 0.9999 2.0000 0.9999 0.0000 0.000000 0.000000
convergió.. NO HA ó.. Ya lo NO HA
Ya lo tiene convergido.. tiene convergido
430 0055 9880 0173 9945 0000 03 02 0.00000001 pequeño continúe pequeño .. continúe
Ya
- Ya convergi
-
3.0000 0.9999 2.0000 0.9999 0.0000 0.000000 0.000000
convergió.. NO HA ó.. Ya lo NO HA
Ya lo tiene convergido.. tiene convergido
431 0054 9883 0168 9947 0000 03 02 0.00000001 pequeño continúe pequeño .. continúe
Ya
- Ya convergi
-
3.0000 0.9999 2.0000 0.9999 0.0000 0.000000 0.000000
convergió.. NO HA ó.. Ya lo NO HA
Ya lo tiene convergido.. tiene convergido
432 0053 9886 0164 9948 0000 03 02 0.00000001 pequeño continúe pequeño .. continúe
Ya
- Ya convergi
-
3.0000 0.9999 2.0000 0.9999 0.0000 0.000000 0.000000
convergió.. NO HA ó.. Ya lo NO HA
Ya lo tiene convergido.. tiene convergido
433 0051 9889 0160 9949 0000 03 02 0.00000001 pequeño continúe pequeño .. continúe
Ya
- Ya convergi
-
3.0000 0.9999 2.0000 0.9999 0.0000 0.000000 0.000000
convergió.. NO HA ó.. Ya lo NO HA
Ya lo tiene convergido.. tiene convergido
434 0050 9892 0156 9951 0000 03 02 0.00000001 pequeño continúe pequeño .. continúe
Ya
- Ya convergi
-
3.0000 0.9999 2.0000 0.9999 0.0000 0.000000 0.000000
convergió.. NO HA ó.. Ya lo NO HA
Ya lo tiene convergido.. tiene convergido
435 0049 9895 0152 9952 0000 03 02 0.00000001 pequeño continúe pequeño .. continúe
Ya
- Ya convergi
-
3.0000 0.9999 2.0000 0.9999 0.0000 0.000000 0.000000
convergió.. NO HA ó.. Ya lo NO HA
Ya lo tiene convergido.. tiene convergido
436 0048 9897 0148 9953 0000 03 02 0.00000001 pequeño continúe pequeño .. continúe
Ya
- Ya convergi
-
3.0000 0.9999 2.0000 0.9999 0.0000 0.000000 0.000000
convergió.. NO HA ó.. Ya lo NO HA
Ya lo tiene convergido.. tiene convergido
437 0046 9900 0145 9954 0000 03 02 0.00000001 pequeño continúe pequeño .. continúe
Ya
- Ya convergi
-
3.0000 0.9999 2.0000 0.9999 0.0000 0.000000 0.000000
convergió.. NO HA ó.. Ya lo NO HA
Ya lo tiene convergido.. tiene convergido
438 0045 9902 0141 9955 0000 03 02 0.00000001 pequeño continúe pequeño .. continúe
Ya
- Ya convergi
-
3.0000 0.9999 2.0000 0.9999 0.0000 0.000000 0.000000
convergió.. NO HA ó.. Ya lo NO HA
Ya lo tiene convergido.. tiene convergido
439 0044 9905 0137 9956 0000 02 02 0.00000001 pequeño continúe pequeño .. continúe
Ya
- Ya convergi
-
3.0000 0.9999 2.0000 0.9999 0.0000 0.000000 0.000000
convergió.. NO HA ó.. Ya lo NO HA
Ya lo tiene convergido.. tiene convergido
440 0043 9907 0134 9958 0000 02 02 0.00000001 pequeño continúe pequeño .. continúe
Ya
- Ya convergi
-
3.0000 0.9999 2.0000 0.9999 0.0000 0.000000 0.000000
convergió.. NO HA ó.. Ya lo NO HA
Ya lo tiene convergido.. tiene convergido
441 0042 9910 0131 9959 0000 02 02 0.00000001 pequeño continúe pequeño .. continúe
0.0000 0.000000 - Ya NO HA Ya NO HA
442 3.0000 0.9999 - 0.9999 0000 02 0.000000 0.00000001 convergió.. convergido.. convergi convergido

pág. 19
 0041 9912 2.0000 9960 02 Ya lo tiene continúe ó.. Ya lo .. continúe
pequeño tiene
0127 pequeño
Ya
- Ya convergi
-
3.0000 0.9999 2.0000 0.9999 0.0000 0.000000 0.000000
convergió.. NO HA ó.. Ya lo NO HA
Ya lo tiene convergido.. tiene convergido
443 0040 9914 0124 9961 0000 02 02 0.00000001 pequeño continúe pequeño .. continúe
Ya
- Ya convergi Ya
-
3.0000 0.9999 2.0000 0.9999 0.0000 0.000000 0.000000
convergió.. NO HA ó.. Ya lo convergió..
Ya lo tiene convergido.. tiene Ya lo tiene
444 0039 9916 0121 9962 0000 02 02 0.00000001 pequeño continúe pequeño pequeño
Ya
- Ya convergi Ya
-
3.0000 0.9999 2.0000 0.9999 0.0000 0.000000 0.000000
convergió.. NO HA ó.. Ya lo convergió..
Ya lo tiene convergido.. tiene Ya lo tiene
445 0038 9918 0118 9963 0000 02 02 0.00000001 pequeño continúe pequeño pequeño
Ya
- Ya convergi Ya
-
3.0000 0.9999 2.0000 0.9999 0.0000 0.000000 0.000000
convergió.. NO HA ó.. Ya lo convergió..
Ya lo tiene convergido.. tiene Ya lo tiene
446 0037 9920 0115 9964 0000 02 01 0.00000001 pequeño continúe pequeño pequeño
Ya
- Ya convergi Ya
-
3.0000 0.9999 2.0000 0.9999 0.0000 0.000000 0.000000
convergió.. NO HA ó.. Ya lo convergió..
Ya lo tiene convergido.. tiene Ya lo tiene
447 0036 9922 0112 9964 0000 02 01 0.00000001 pequeño continúe pequeño pequeño
Ya
- Ya convergi Ya
-
3.0000 0.9999 2.0000 0.9999 0.0000 0.000000 0.000000
convergió.. NO HA ó.. Ya lo convergió..
Ya lo tiene convergido.. tiene Ya lo tiene
448 0035 9924 0109 9965 0000 02 01 0.00000001 pequeño continúe pequeño pequeño
Ya
- Ya convergi Ya
-
3.0000 0.9999 2.0000 0.9999 0.0000 0.000000 0.000000
convergió.. NO HA ó.. Ya lo convergió..
Ya lo tiene convergido.. tiene Ya lo tiene
449 0034 9926 0107 9966 0000 02 01 0.00000001 pequeño continúe pequeño pequeño
Ya
- Ya convergi Ya
-
3.0000 0.9999 2.0000 0.9999 0.0000 0.000000 0.000000
convergió.. NO HA ó.. Ya lo convergió..
Ya lo tiene convergido.. tiene Ya lo tiene
450 0033 9928 0104 9967 0000 02 01 0.00000001 pequeño continúe pequeño pequeño
Ya
- Ya convergi Ya
-
3.0000 0.9999 2.0000 0.9999 0.0000 0.000000 0.000000
convergió.. NO HA ó.. Ya lo convergió..
Ya lo tiene convergido.. tiene Ya lo tiene
451 0033 9930 0101 9968 0000 02 01 0.00000001 pequeño continúe pequeño pequeño
Ya
- Ya convergi Ya
-
3.0000 0.9999 2.0000 0.9999 0.0000 0.000000 0.000000
convergió.. NO HA ó.. Ya lo convergió..
Ya lo tiene convergido.. tiene Ya lo tiene
452 0032 9932 0099 9969 0000 02 01 0.00000001 pequeño continúe pequeño pequeño
Ya
- Ya convergi Ya
-
3.0000 0.9999 2.0000 0.9999 0.0000 0.000000 0.000000
convergió.. NO HA ó.. Ya lo convergió..
Ya lo tiene convergido.. tiene Ya lo tiene
453 0031 9933 0096 9969 0000 02 01 0.00000001 pequeño continúe pequeño pequeño
Ya
- Ya convergi Ya
-
3.0000 0.9999 2.0000 0.9999 0.0000 0.000000 0.000000
convergió.. NO HA ó.. Ya lo convergió..
Ya lo tiene convergido.. tiene Ya lo tiene
454 0030 9935 0094 9970 0000 02 01 0.00000001 pequeño continúe pequeño pequeño
Ya
- Ya convergi Ya
-
3.0000 0.9999 2.0000 0.9999 0.0000 0.000000 0.000000
convergió.. NO HA ó.. Ya lo convergió..
Ya lo tiene convergido.. tiene Ya lo tiene
455 0029 9937 0092 9971 0000 02 01 0.00000001 pequeño continúe pequeño pequeño
Ya
- Ya convergi Ya
-
3.0000 0.9999 2.0000 0.9999 0.0000 0.000000 0.000000
convergió.. NO HA ó.. Ya lo convergió..
Ya lo tiene convergido.. tiene Ya lo tiene
456 0029 9938 0089 9972 0000 02 01 0.00000001 pequeño continúe pequeño pequeño
Ya
- Ya convergi Ya
-
3.0000 0.9999 2.0000 0.9999 0.0000 0.000000 0.000000
convergió.. NO HA ó.. Ya lo convergió..
Ya lo tiene convergido.. tiene Ya lo tiene
457 0028 9940 0087 9972 0000 02 01 0.00000001 pequeño continúe pequeño pequeño
Ya
- Ya convergi Ya
-
3.0000 0.9999 2.0000 0.9999 0.0000 0.000000 0.000000
convergió.. NO HA ó.. Ya lo convergió..
Ya lo tiene convergido.. tiene Ya lo tiene
458 0027 9941 0085 9973 0000 02 01 0.00000001 pequeño continúe pequeño pequeño
Ya
- Ya convergi Ya
-
3.0000 0.9999 2.0000 0.9999 0.0000 0.000000 0.000000
convergió.. NO HA ó.. Ya lo convergió..
Ya lo tiene convergido.. tiene Ya lo tiene
459 0027 9943 0083 9974 0000 01 01 0.00000001 pequeño continúe pequeño pequeño
Ya
- Ya convergi Ya
-
3.0000 0.9999 2.0000 0.9999 0.0000 0.000000 0.000000
convergió.. NO HA ó.. Ya lo convergió..
Ya lo tiene convergido.. tiene Ya lo tiene
460 0026 9944 0081 9974 0000 01 01 0.00000001 pequeño continúe pequeño pequeño
Ya
- Ya convergi Ya
-
3.0000 0.9999 2.0000 0.9999 0.0000 0.000000 0.000000
convergió.. NO HA ó.. Ya lo convergió..
Ya lo tiene convergido.. tiene Ya lo tiene
461 0025 9946 0079 9975 0000 01 01 0.00000001 pequeño continúe pequeño pequeño
Ya
- Ya convergi Ya
-
3.0000 0.9999 2.0000 0.9999 0.0000 0.000000 0.000000
convergió.. NO HA ó.. Ya lo convergió..
Ya lo tiene convergido.. tiene Ya lo tiene
462 0025 9947 0077 9976 0000 01 01 0.00000001 pequeño continúe pequeño pequeño

pág. 20
 Ya
- Ya convergi Ya
-
3.0000 0.9999 2.0000 0.9999 0.0000 0.000000 0.000000
convergió.. NO HA ó.. Ya lo convergió..
Ya lo tiene convergido.. tiene Ya lo tiene
463 0024 9948 0075 9976 0000 01 01 0.00000001 pequeño continúe pequeño pequeño
Ya
- Ya convergi Ya
-
3.0000 0.9999 2.0000 0.9999 0.0000 0.000000 0.000000
convergió.. NO HA ó.. Ya lo convergió..
Ya lo tiene convergido.. tiene Ya lo tiene
464 0023 9949 0073 9977 0000 01 01 0.00000001 pequeño continúe pequeño pequeño
Ya
- Ya convergi Ya
-
3.0000 0.9999 2.0000 0.9999 0.0000 0.000000 0.000000
convergió.. NO HA ó.. Ya lo convergió..
Ya lo tiene convergido.. tiene Ya lo tiene
465 0023 9951 0071 9977 0000 01 01 0.00000001 pequeño continúe pequeño pequeño
Ya
- Ya convergi Ya
-
3.0000 0.9999 2.0000 0.9999 0.0000 0.000000 0.000000
convergió.. NO HA ó.. Ya lo convergió..
Ya lo tiene convergido.. tiene Ya lo tiene
466 0022 9952 0069 9978 0000 01 01 0.00000001 pequeño continúe pequeño pequeño
Ya
- Ya convergi Ya
-
3.0000 0.9999 2.0000 0.9999 0.0000 0.000000 0.000000
convergió.. NO HA ó.. Ya lo convergió..
Ya lo tiene convergido.. tiene Ya lo tiene
467 0022 9953 0068 9979 0000 01 01 0.00000001 pequeño continúe pequeño pequeño
Ya
- Ya convergi Ya
-
3.0000 0.9999 2.0000 0.9999 0.0000 0.000000 0.000000
convergió.. NO HA ó.. Ya lo convergió..
Ya lo tiene convergido.. tiene Ya lo tiene
468 0021 9954 0066 9979 0000 01 01 0.00000001 pequeño continúe pequeño pequeño
Ya
- Ya convergi Ya
-
3.0000 0.9999 2.0000 0.9999 0.0000 0.000000 0.000000
convergió.. NO HA ó.. Ya lo convergió..
Ya lo tiene convergido.. tiene Ya lo tiene
469 0021 9956 0064 9980 0000 01 01 0.00000001 pequeño continúe pequeño pequeño
Ya
- Ya convergi Ya
-
3.0000 0.9999 2.0000 0.9999 0.0000 0.000000 0.000000
convergió.. NO HA ó.. Ya lo convergió..
Ya lo tiene convergido.. tiene Ya lo tiene
470 0020 9957 0063 9980 0000 01 01 0.00000001 pequeño continúe pequeño pequeño
Ya
- Ya convergi Ya
-
3.0000 0.9999 2.0000 0.9999 0.0000 0.000000 0.000000
convergió.. NO HA ó.. Ya lo convergió..
Ya lo tiene convergido.. tiene Ya lo tiene
471 0020 9958 0061 9981 0000 01 01 0.00000000 pequeño continúe pequeño pequeño
Ya
- Ya convergi Ya
-
3.0000 0.9999 2.0000 0.9999 0.0000 0.000000 0.000000
convergió.. NO HA ó.. Ya lo convergió..
Ya lo tiene convergido.. tiene Ya lo tiene
472 0019 9959 0059 9981 0000 01 01 0.00000000 pequeño continúe pequeño pequeño
Ya
- Ya convergi Ya
-
3.0000 0.9999 2.0000 0.9999 0.0000 0.000000 0.000000
convergió.. NO HA ó.. Ya lo convergió..
Ya lo tiene convergido.. tiene Ya lo tiene
473 0019 9960 0058 9982 0000 01 01 0.00000000 pequeño continúe pequeño pequeño
Ya
- Ya convergi Ya
-
3.0000 0.9999 2.0000 0.9999 0.0000 0.000000 0.000000
convergió.. NO HA ó.. Ya lo convergió..
Ya lo tiene convergido.. tiene Ya lo tiene
474 0018 9961 0057 9982 0000 01 01 0.00000000 pequeño continúe pequeño pequeño
Ya
- Ya convergi Ya
-
3.0000 0.9999 2.0000 0.9999 0.0000 0.000000 0.000000
convergió.. Ya convergió.. ó.. Ya lo convergió..
Ya lo tiene Ya lo tiene tiene Ya lo tiene
475 0018 9962 0055 9983 0000 01 01 0.00000000 pequeño pequeño pequeño pequeño

Resultados:
𝒙𝟏 = 𝟑. 𝟎𝟎𝟎𝟎𝟎𝟎𝟏𝟖 ≅ 𝟑
𝒙𝟐 = 𝟎. 𝟗𝟗𝟗𝟗𝟗𝟗𝟔𝟐 ≅ 𝟏
𝒙𝟑 = −𝟐. 𝟎𝟎𝟎𝟎𝟎𝟎𝟓𝟓 ≅ −𝟐
𝒙𝟒 = 𝟎. 𝟗𝟗𝟗𝟗𝟗𝟗𝟖𝟑 ≅ 𝟏

pág. 21
Jacobi:

El método Jacobi es iterativo, la cual permite la solución de sistemas lineales

grandes; este método es igual al anterior (gauss-seidel) con la diferencia de que el
valor de las incógnitas no se van tomando de forma inmediata, sino hasta la
siguiente iteración.

Las ecuaciones utilizadas para para un sistema SEL de 4x4, es de:

𝑏1 − 𝑎12 𝑥2 − 𝑎13 𝑥3 − 𝑎14 𝑥4

𝑥1 =
𝑎11

𝑏2 − 𝑎21 𝑥1 − 𝑎23 𝑥3 − 𝑎24 𝑥4

𝑥2 =
𝑎22

𝑏3 − 𝑎31 𝑥1 − 𝑎32 𝑥2 − 𝑎34 𝑥4

𝑥3 =
𝑎33

𝑏4 − 𝑎41 𝑥1 − 𝑎42 𝑥2 − 𝑎43 𝑥3

𝑥4 =
𝑎44

Para determinar si el método jacobi puede tener solución, el SEL debe garantizar
una matriz dominante:

“…el coeficiente diagonal en cada una de las ecuaciones debe ser mayor que la
suma del valor absoluto de los otros coeficientes de la ecuación. Este criterio es
suficiente pero no necesario para la convergencia. Es decir, el método puede
funcionar aunque no se satisfaga la convergencia se garantiza cuando la condición
se satisface. A los sistemas que cumplen con la ecuación se les conoce como
diagonalmente dominantes.“ (Chapra, 2007, p.315)

Veamos si esto se cumple… tomemos la ecuación original:

3𝑥1 − 13𝑥2 + 9𝑥3 + 3𝑥4 = −19
−6𝑥1 + 4𝑥2 + 𝑥3 − 18𝑥4 = −34
6𝑥1 − 2𝑥2 + 2𝑥3 + 4𝑥4 = 16
12𝑥1 − 8𝑥2 + 6𝑥3 + 10𝑥4 = 26

pág. 22
Estableciendo en matriz ampliada:

3 −13 9 3 −19
−6 4 1 −18 −34
| || |
6 −2 2 4 16
12 −8 6 10 26

La diagonal sería 3,4,2,10

Observamos que NO es una diagonal dominante ya que:

3 < (13 + 9 + 3) → 𝑁𝑜 𝑠𝑒 𝑐𝑢𝑚𝑝𝑙𝑒 𝑦𝑎 𝑞𝑢𝑒 3 𝑑𝑒𝑏𝑒𝑟í𝑎 𝑠𝑒𝑟 𝑚𝑎𝑦𝑜𝑟 𝑜 𝑖𝑔𝑢𝑎𝑙.

4 < (6 + 1 + 18) → 𝑁𝑜 𝑠𝑒 𝑐𝑢𝑚𝑝𝑙𝑒 𝑦𝑎 𝑞𝑢𝑒 4 𝑑𝑒𝑏𝑒𝑟í𝑎 𝑠𝑒𝑟 𝑚𝑎𝑦𝑜𝑟 𝑜 𝑖𝑔𝑢𝑎𝑙.
2 < (6 + 2 + 4) → 𝑁𝑜 𝑠𝑒 𝑐𝑢𝑚𝑝𝑙𝑒 𝑦𝑎 𝑞𝑢𝑒 2 𝑑𝑒𝑏𝑒𝑟í𝑎 𝑠𝑒𝑟 𝑚𝑎𝑦𝑜𝑟 𝑜 𝑖𝑔𝑢𝑎𝑙.
10 < (12 + 8 + 6) → 𝑁𝑜 𝑠𝑒 𝑐𝑢𝑚𝑝𝑙𝑒 𝑦𝑎 𝑞𝑢𝑒 6 𝑑𝑒𝑏𝑒𝑟í𝑎 𝑠𝑒𝑟 𝑚𝑎𝑦𝑜𝑟 𝑜 𝑖𝑔𝑢𝑎𝑙.

Nota: como pueden observar, únicamente para la verificación he omitido los

signos; esto es insignificante para ello según la teoría... se toman términos
como suma, independientemente si son + o -.

Antes de tomar una decisión, el teorema de re-arreglo de ecuaciones también aplica

para jacobi. Si re-arreglando la ecuación podemos conseguir la diagonal dominante,
podremos tener un margen de éxito en la convergencia al resultado… veamos.

Mejor combinación posible en Re-Arreglo:

12 −8 6 10 26
3 −13 9 3 −19
| || |
6 −2 2 4 16
−6 4 1 −18 −34

Volvamos a Revisar:

12 < (8 + 6 + 10) → 𝑁𝑜 𝑠𝑒 𝑐𝑢𝑚𝑝𝑙𝑒 𝑦𝑎 𝑞𝑢𝑒 12 𝑑𝑒𝑏𝑒𝑟í𝑎 𝑠𝑒𝑟 𝑚𝑎𝑦𝑜𝑟 𝑜 𝑖𝑔𝑢𝑎𝑙.

13 < (3 + 9 + 3) → 𝑁𝑜 𝑠𝑒 𝑐𝑢𝑚𝑝𝑙𝑒 𝑦𝑎 𝑞𝑢𝑒 13 𝑑𝑒𝑏𝑒𝑟í𝑎 𝑠𝑒𝑟 𝑚𝑎𝑦𝑜𝑟 𝑜 𝑖𝑔𝑢𝑎𝑙.
2 < (6 + 2 + 4) → 𝑁𝑜 𝑠𝑒 𝑐𝑢𝑚𝑝𝑙𝑒 𝑦𝑎 𝑞𝑢𝑒 2 𝑑𝑒𝑏𝑒𝑟í𝑎 𝑠𝑒𝑟 𝑚𝑎𝑦𝑜𝑟 𝑜 𝑖𝑔𝑢𝑎𝑙.
18 > (6 + 4 + 1) → 𝑆𝐼 𝑠𝑒 𝑐𝑢𝑚𝑝𝑙𝑒 𝑦𝑎 𝑞𝑢𝑒 18 𝑑𝑒𝑏𝑒𝑟í𝑎 𝑠𝑒𝑟 𝑚𝑎𝑦𝑜𝑟 𝑜 𝑖𝑔𝑢𝑎𝑙.

Definitivamente comprobamos que la ecuación NO tiene diagonal dominante; esto

obedece que solo 1 ecuación de las 4, tiene 1 termino el doble o mayor de grande
que el resto…. Pero se requieren que todos los 4 sean así; de todos modos y con el
ánimo de hacer deporte, explicaremos el teorema, aun sabiendo que lo más seguro
que el método iterativo de jacobi, no será el adecuado para que converja a la
solución de las incógnitas.

pág. 23
Teniendo en cuenta que no tenemos valores iniciales de x, iniciaremos con Cero en
todos los valores ocultos:

𝑥1 = 0
𝑥2 = 0
𝑥3 = 0
𝑥4 = 0

12 −8 6 10 26
3 −13 9 3 −19
| || |
6 −2 2 4 16
−6 4 1 −18 −34

Iniciamos la búsqueda de los valores ocultos con las formulas; usaremos el número
de iteraciones que sean necesarios, hasta que converjan (eso si tuviéramos una
diagonal dominante…. Como no es así, nos detendremos cuando me aburra).

26 − (−8𝑥2 ) − (6𝑥3 ) − (10𝑥4 )

𝑥1 =
12

−19 − (3𝑥1 ) − (9𝑥3 ) − (3𝑥4 )

𝑥2 =
−13

16 − (6𝑥1 ) − (−2𝑥2 ) − (4𝑥4 )

𝑥3 =
2

−34 − (−6𝑥1 ) − (4𝑥2 ) − (1𝑥3 )

𝑥4 =
−18

pág. 24
Iteración 1:

𝑥1 = 0
𝑥2 = 0
𝑥3 = 0
𝑥4 = 0

26 − (−8 ∗ 0) − (6 ∗ 0) − (10 ∗ 0)
𝑥1 =
12
26
𝑥1 = = 2.16666666666
12

−19 − (3 ∗ 0) − (9 ∗ 0) − (3 ∗ 0)
𝑥2 =
−13

−19
𝑥2 = = 1.46153846153
−13

16 − (6 ∗ 0) − (−2 ∗ 0) − (4 ∗ 0)
𝑥3 =
2

16
𝑥3 = =8
2

−34 − (−6 ∗ 0) − (4 ∗ 0) − (1 ∗ 0)
𝑥4 =
−18

−34
𝑥4 = = 1.88888888888
−18

Hasta este punto, tenemos los nuevos valores de x:

𝑥1 = 2.16666666666
𝑥2 = 1.46153846153
𝑥3 = 8
𝑥4 = 1.88888888888

Seguimos con la siguiente iteración y repetimos.

pág. 25
Iteración 2:

𝑥1 = 2.16666666666
𝑥2 = 1.46153846153
𝑥3 = 8
𝑥4 = 1.88888888888

26 − (−8 ∗ 1.46153846153) − (6 ∗ 8) − (10 ∗ 1.88888888888)

𝑥1 =
12

26 + 11.6923076922 − 48 − 18.8888888888
𝑥1 = = −2.43304843304843
12

−19 − (3 ∗ 2.16666666666) − (9 ∗ 8) − (3 ∗ 1.88888888888)

𝑥2 =
−13

−19 − 6.49999999999 − 72 − 5.66666666661

𝑥2 = = 7.93589743589744
−13

16 − (6 ∗ 2.16666666666) − (−2 ∗ 1.46153846153) − (4 ∗ 1.88888888888)

𝑥3 =
2

16 − 12.999999999 + 2.92307692306 − 7.55555555552

𝑥3 = = −0.816239316239316
2

−34 − (−6 ∗ 2.16666666666) − (4 ∗ 1.46153846153) − (1 ∗ 8)

𝑥4 =
−18

−34 + 12.9999999999 − 5.84615384612 − 8

𝑥4 = = 1.93589743589744
−18

Hasta este punto, tenemos los nuevos valores de x:

𝑥1 = −2.43304843304843
𝑥2 = 7.93589743589744
𝑥3 = −0.816239316239316
𝑥4 = 1.93589743589744

pág. 26
Cada vez que se termine una iteración, debe verificarse el error porcentual para
validar si los 4 valores han convergido, esto en función de ‖𝐴‖ ≤ 10−8 ; si no han
convergido, se debe continuar con las iteraciones.

Tabla de resultados Jacobi:

|(𝑥)𝑎𝑐𝑡𝑢𝑎𝑙 − (𝑥)𝑎𝑛𝑡𝑒𝑟|
𝐸𝑟𝑟𝑜𝑟% = 𝐸% = ∗ 100%
(𝑥)𝑎𝑐𝑡𝑢𝑎𝑙

𝑃𝑜𝑟𝑐𝑒𝑛𝑡𝑎𝑗𝑒 𝑑𝑒 𝑒𝑟𝑟𝑜𝑟 𝑒𝑛 𝑓𝑢𝑛𝑐𝑖𝑜𝑛 𝑐𝑜𝑛 𝑠𝑢 𝑙𝑖𝑚𝑖𝑡𝑒 𝑎 𝑐𝑒𝑟𝑜.

𝐸% > 𝑇𝑜𝑙𝑒𝑟𝑎𝑛𝑐𝑖𝑎 = 𝐶𝑜𝑛𝑣𝑒𝑟𝑔𝑒; 𝑑𝑒 𝑙𝑜 𝑐𝑜𝑛𝑡𝑟𝑎𝑟𝑖𝑜 𝑁𝑂 𝑐𝑜𝑛𝑣𝑒𝑟𝑔𝑒

𝑇𝑜𝑙𝑒𝑟𝑎𝑛𝑐𝑖𝑎 𝑠𝑒𝑔ú𝑛 𝐶é𝑑𝑢𝑙𝑎 = 10−8 = 0.00000001

Iter
acio Error( Error(% Error(%
n x1 x2 x3 x4 %)x1 )x2 )x3 Error(%)x4 Rx1 Rx2 Rx3 Rx4
0.0000 0.0000 0.000000 0.000000
0000 0000 00 00
2.1666 1.4615 8.000000 1.888888
1 6667 3846 00 89
Ya
convergi
- - - - Ya convergió.. NO HA ó.. Ya lo NO HA
2.4330 7.9358 0.816239 1.935897 1.8905 0.815831 10.80104 Ya lo tiene convergido.. tiene convergido
2 4843 9744 32 44 1522 99 712 0.02428256 pequeño continúe pequeño .. continúe
NO HA
NO HA NO HA convergi NO HA
6.2521 0.7817 19.36324 4.418091 1.3891 9.151808 1.042154 convergido.. convergido.. do.. convergido
3 3675 2255 786 17 5470 24 05 0.56182492 continúe continúe continúe .. continúe
Ya
convergi
- - - - Ya convergió.. NO HA ó.. Ya lo NO HA
10.675 17.329 18.81087 1.054295 1.5856 0.954889 2.029364 Ya lo tiene convergido.. tiene convergido
4 55154 22419 004 42 5000 93 82 3.19056280 pequeño continúe pequeño .. continúe
NO HA
- - NO HA Ya convergió.. convergi NO HA
22.246 13.781 55.24728 8.253296 1.4798 2.257411 1.340484 convergido.. Ya lo tiene do.. convergido
5 33830 66144 797 44 7904 83 95 0.87225766 continúe pequeño continúe .. continúe
Ya
convergi Ya
- - - - - Ya convergió.. NO HA ó.. Ya lo convergió..
41.522 46.748 89.02726 5.519854 1.5357 1.294807 1.620565 Ya lo tiene convergido.. tiene Ya lo tiene
6 49865 03815 923 87 6589 27 91 -2.49520171 pequeño continúe pequeño pequeño

Resultados:
No se llega a convergencia; método no es el adecuado.

pág. 27
Gauss-seidel relajado y/o S.O.R:

“La relajación representa una ligera modificación al método de Gauss-Seidel y ésta

permite mejorar la convergencia. Después de que se calcula cada nuevo valor de x,
ese valor se modifica mediante un promedio ponderado de los resultados de las
iteraciones anterior y actual.” (Chapra, 2007, p.315)

𝑥𝑖𝑛𝑢𝑒𝑣𝑜 = 𝜆𝑥𝑖𝑛𝑢𝑒𝑣𝑜 + (1 − 𝜆)𝑥𝑖𝑎𝑛𝑡𝑒𝑟𝑖𝑜𝑟

Tomemos lo ya realizado de Gauss-seidel, y ahora apliquemos la ponderación.

Las ecuaciones utilizadas para para un sistema SEL de 4x4, es de:

𝑏1 − 𝑎12 𝑥2 − 𝑎13 𝑥3 − 𝑎14 𝑥4

𝑥1 =
𝑎11

𝑏2 − 𝑎21 𝑥1 − 𝑎23 𝑥3 − 𝑎24 𝑥4

𝑥2 =
𝑎22

𝑏3 − 𝑎31 𝑥1 − 𝑎32 𝑥2 − 𝑎34 𝑥4

𝑥3 =
𝑎33

𝑏4 − 𝑎41 𝑥1 − 𝑎42 𝑥2 − 𝑎43 𝑥3

𝑥4 =
𝑎44

Teniendo en cuenta que no tenemos valores iniciales de x, iniciaremos con Cero en

todos los valores ocultos:

𝑥1 = 0
𝑥2 = 0
𝑥3 = 0
𝑥4 = 0

Aplicando factor ponderado de 0.99

𝜆 = 0.99

Iniciamos la búsqueda de los valores ocultos con las formulas; usaremos el número
de iteraciones que sean necesarios, hasta que converjan.

pág. 28
Realizando ejercicio en orden de:
6𝑥1 − 2𝑥2 + 2𝑥3 + 4𝑥4 = 16
12𝑥1 − 8𝑥2 + 6𝑥3 + 10𝑥4 = 26
3𝑥1 − 13𝑥2 + 9𝑥3 + 3𝑥4 = −19
−6𝑥1 + 4𝑥2 + 𝑥3 − 18𝑥4 = −34

Estableciendo en matriz ampliada:

6 −2 2 4 16
12 −8 6 10 26
| || |
3 −13 9 3 −19
−6 4 1 −18 −34

𝑥1 = 0
𝑥2 = 0
𝑥3 = 0
𝑥4 = 0

16 − (−2𝑥2 ) − (2𝑥3 ) − (4𝑥4 )

𝑥1 =
6

26 − (12𝑥1 ) − (6𝑥3 ) − (10𝑥4 )

𝑥2 =
−8

−19 − (3𝑥1 ) − (−13𝑥2 ) − (3𝑥4 )

𝑥3 =
9

−34 − (−6𝑥1 ) − (4𝑥2 ) − (1𝑥3 )

𝑥4 =
−18

Iteración 1:

16 − (−2 ∗ 0) − (2 ∗ 0) − (4 ∗ 0) 16 + 0 − 0 − 0
𝑥1 = = = 2.66666666
6 6

Aplicamos relajación:

𝑥1 = (0.99 ∗ 2.66666666) + (1 − 0.99) ∗ 0

𝑥1 = 2.64

pág. 29
Ojo.. ya tenemos un valor de x1 relajado, por lo tanto debemos usarlo para la
siguiente ecuación… así vamos haciendo en todo.. vamos buscando resultados de x
y los vamos usando inmediatamente.

26 − (12 ∗ 2.64) − (6 ∗ 0) − (10 ∗ 0)

𝑥2 =
−8

26 − 31.68 − 0 − 0
𝑥2 = = 0.71
−8

Aplicamos relajación:

𝑥2 = (0.99 ∗ 0.71) + (1 − 0.99) ∗ 0

𝑥2 = 0.7029

−19 − (3 ∗ 2.64) − (−13 ∗ 0.7029) − (3 ∗ 0)

𝑥3 =
9

−19 − 7.92 + 9.1377 − 0

𝑥3 = = −1.97581111111
9

Aplicamos relajación:

𝑥3 = (0.99 ∗ 1.97581111111) + (1 − 0.99) ∗ 0

𝑥3 = −1.95605299999

−34 − (−6 ∗ 2.64) − (4 ∗ 0.7029) − (1 ∗ −1.95605299999)

𝑥4 =
−18

−34 + 15.84 − 2.8116 + 1.95605299999

𝑥4 = = 1.05641927778
−18

Aplicamos relajación:

𝑥4 = (0.99 ∗ 1.05641927778) + (1 − 0.99) ∗ 0

𝑥4 = 1.045855085

pág. 30
Hasta este punto, tenemos los valores de x:

𝑥1 = 2.64
𝑥2 = 0.7029
𝑥3 = −1.95605299999
𝑥4 = 1.045855085

Seguimos con la siguiente iteración y repetimos.

Iteración 2:

𝑥1 = 2.64
𝑥2 = 0.7029
𝑥3 = −1.95605299999
𝑥4 = 1.045855085

16 − (−2 ∗ 0.7029) − (2 ∗ −1.95605299999) − (4 ∗ 1.045855085)

𝑥1 =
6

16 + 1.4058 + 3.91210599998 − 4.18342034

𝑥1 = = 2.85574760998
6

Aplicamos relajación:

𝑥1 = (0.99 ∗ 2.85574760998) + (1 − 0.99) ∗ 2.64

𝑥1 = 2.8535901339

Recuerden que los nuevos resultados de x se usan inmediatamente.

26 − (12 ∗ 2.8535901339) − (6 ∗ −1.95605299999) − (10 ∗ 1.045855085)

𝑥2 =
−8

26 − 34.2430816068 + 11.7363179999 − 10.45855085

𝑥2 = = 0.87066430711
−8

Aplicamos relajación:

𝑥2 = (0.99 ∗ 0.87066430711) + (1 − 0.99) ∗ 0.7029

𝑥2 = 0.868986664029

pág. 31
 −19 − (3 ∗ 2.8535901339) − (−13 ∗ 0.868986664029) − (3 ∗ 1.045855085)
𝑥3 =
9

−19 − 8.5607704197 + 11.2968266322 − 3.137565255

𝑥3 = = −2.15572322694
9

Aplicamos relajación:

𝑥3 = (0.99 ∗ −2.15572322694) + (1 − 0.99) ∗ −1.95605299999

𝑥3 = −2.15372652267553

−34 − (−6 ∗ 2.8535901339) − (4 ∗ 0.868986664029) − (1 ∗ −2.15372652267553)

𝑥4 =
−18

−34 + 17.1215408034 − 3.47594665608 + 2.15372652267553

𝑥4 = = 1.01114885166
−18

Aplicamos relajación:

𝑥4 = (0.99 ∗ 1.01114885166) + (1 − 0.99) ∗ 1.045855085

𝑥4 = 1.01149591400223

Hasta este punto, tenemos los valores de x:

𝑥1 = 2.8535901339
𝑥2 = 0.868986664029
𝑥3 = −2.15372652267553
𝑥4 = 1.01149591400223

Cada vez que se termine una iteración, debe verificarse el error porcentual para
validar si los 4 valores han convergido, esto en función de ‖𝐴‖ ≤ 10−8 ; si no han
convergido, se debe continuar con las iteraciones.

pág. 32
Tabla de resultados Gauss – Seidel relajado y/o SOR:

|(𝑥)𝑎𝑐𝑡𝑢𝑎𝑙 − (𝑥)𝑎𝑛𝑡𝑒𝑟|
𝐸𝑟𝑟𝑜𝑟% = 𝐸% = ∗ 100%
(𝑥)𝑎𝑐𝑡𝑢𝑎𝑙

𝑃𝑜𝑟𝑐𝑒𝑛𝑡𝑎𝑗𝑒 𝑑𝑒 𝑒𝑟𝑟𝑜𝑟 𝑒𝑛 𝑓𝑢𝑛𝑐𝑖𝑜𝑛 𝑐𝑜𝑛 𝑠𝑢 𝑙𝑖𝑚𝑖𝑡𝑒 𝑎 𝑐𝑒𝑟𝑜.

𝐸% > 𝑇𝑜𝑙𝑒𝑟𝑎𝑛𝑐𝑖𝑎 = 𝐶𝑜𝑛𝑣𝑒𝑟𝑔𝑒; 𝑑𝑒 𝑙𝑜 𝑐𝑜𝑛𝑡𝑟𝑎𝑟𝑖𝑜 𝑁𝑂 𝑐𝑜𝑛𝑣𝑒𝑟𝑔𝑒

𝑇𝑜𝑙𝑒𝑟𝑎𝑛𝑐𝑖𝑎 𝑠𝑒𝑔ú𝑛 𝐶é𝑑𝑢𝑙𝑎 = 10−8 = 0.00000001

iter Error(
ac x1 x2 x3 x4 Error(%) Error(%) %) Error(%) Evaluacion Evaluacion Evaluacion Evaluacion
0.000000 0.000000 0.000000 0.0000
0 00 00 00 0000
- -
2.640000 0.702900 1.956053 1.0458 1.000000 1.000000 1.0000 1.000000 NO HA Ya convergió.. NO HA
NO HA convergido.. convergido.. Ya lo tiene convergido..
1 00 00 00 5509 00 00 0000 00 continúe continúe pequeño continúe
- -
2.853590 0.868986 2.153726 1.0114 0.074849 0.191126 0.0917 0.033968 NO HA Ya convergió.. NO HA
NO HA convergido.. convergido.. Ya lo tiene convergido..
2 13 66 52 9591 62 83 8209 67 continúe continúe pequeño continúe
- -
2.998443 0.896463 2.152874 0.9694 0.048309 0.030650 0.0003 0.043379 NO HA Ya convergió.. NO HA
NO HA convergido.. convergido.. Ya lo tiene convergido..
3 95 38 78 4229 66 13 9563 20 continúe continúe pequeño continúe
- -
3.036434 0.901744 2.143973 0.9581 0.012511 0.005856 0.0041 0.011799 NO HA Ya convergió.. NO HA
NO HA convergido.. convergido.. Ya lo tiene convergido..
4 12 61 17 3645 44 68 5192 82 continúe continúe pequeño continúe
- -
3.043081 0.904286 2.138711 0.9566 0.002184 0.002811 0.0024 0.001523 NO HA Ya convergió.. NO HA
NO HA convergido.. convergido.. Ya lo tiene convergido..
5 15 74 51 7853 31 19 6020 94 continúe continúe pequeño continúe
- -
3.043212 0.906609 2.134899 0.9573 0.000043 0.002562 0.0017 0.000692 NO HA Ya convergió.. NO HA
NO HA convergido.. convergido.. Ya lo tiene convergido..
6 40 67 29 4136 13 22 8567 36 continúe continúe pequeño continúe
- -
3.042284 0.908906 2.131489 0.9583 0.000304 0.002526 0.0015 0.001049 NO HA Ya convergió.. NO HA
NO HA convergido.. convergido.. Ya lo tiene convergido..
7 78 21 74 4687 91 71 9961 21 continúe continúe pequeño continúe
- -
3.041244 0.911160 2.128220 0.9593 0.000342 0.002473 0.0015 0.001072 NO HA Ya convergió.. NO HA
NO HA convergido.. convergido.. Ya lo tiene convergido..
8 58 38 73 7590 03 95 3603 61 continúe continúe pequeño continúe
- -
3.040220 0.913362 2.125040 0.9603 0.000336 0.002410 0.0014 0.001049 NO HA Ya convergió.. NO HA
NO HA convergido.. convergido.. Ya lo tiene convergido..
9 12 26 86 8357 97 74 9638 24 continúe continúe pequeño continúe
- -
3.039222 0.915510 2.121940 0.9613 0.000328 0.002346 0.0014 0.001021 NO HA Ya convergió.. NO HA
NO HA convergido.. convergido.. Ya lo tiene convergido..
10 07 23 64 6607 39 20 6103 98 continúe continúe pequeño continúe
- -
3.038249 0.917605 2.118917 0.9623 0.000320 0.002282 0.0014 0.000995 NO HA Ya convergió.. NO HA
NO HA convergido.. convergido.. Ya lo tiene convergido..
11 40 05 29 2402 14 92 2684 46 continúe continúe pequeño continúe
- -
3.037301 0.919647 2.115968 0.9632 0.000312 0.002221 0.0013 0.000969 NO HA Ya convergió.. NO HA
NO HA convergido.. convergido.. Ya lo tiene convergido..
12 01 94 87 5816 25 38 9341 78 continúe continúe pequeño continúe

13 etc.. etc.. etc.. etc.. etc.. etc.. etc.. etc.. etc.. etc.. etc.. etc..

14 etc.. etc.. etc.. etc.. etc.. etc.. etc.. etc.. etc.. etc.. etc.. etc..

15 etc.. etc.. etc.. etc.. etc.. etc.. etc.. etc.. etc.. etc.. etc.. etc..

pág. 33
16 etc.. etc.. etc.. etc.. etc.. etc.. etc.. etc.. etc.. etc.. etc.. etc..
- -
47 3.000000 0.999999 2.000000 0.9999 0.000000 0.000000 0.0000 0.000000 NO HA Ya convergió..
Ya convergió.. Ya lo convergido.. Ya lo tiene Ya convergió.. Ya
7 32 32 99 9969 00 02 0001 01 tiene pequeño continúe pequeño lo tiene pequeño
- -
47 3.000000 0.999999 2.000000 0.9999 0.000000 0.000000 0.0000 0.000000 NO HA Ya convergió..
Ya convergió.. Ya lo convergido.. Ya lo tiene Ya convergió.. Ya
8 31 33 96 9970 00 02 0001 01 tiene pequeño continúe pequeño lo tiene pequeño
- -
47 3.000000 0.999999 2.000000 0.9999 0.000000 0.000000 0.0000 0.000000 NO HA Ya convergió..
Ya convergió.. Ya lo convergido.. Ya lo tiene Ya convergió.. Ya
9 30 35 94 9970 00 02 0001 01 tiene pequeño continúe pequeño lo tiene pequeño
- -
48 3.000000 0.999999 2.000000 0.9999 0.000000 0.000000 0.0000 0.000000 NO HA Ya convergió..
Ya convergió.. Ya lo convergido.. Ya lo tiene Ya convergió.. Ya
0 29 37 92 9971 00 02 0001 01 tiene pequeño continúe pequeño lo tiene pequeño
- -
48 3.000000 0.999999 2.000000 0.9999 0.000000 0.000000 0.0000 0.000000 NO HA Ya convergió..
Ya convergió.. Ya lo convergido.. Ya lo tiene Ya convergió.. Ya
1 29 38 89 9972 00 02 0001 01 tiene pequeño continúe pequeño lo tiene pequeño
- -
48 3.000000 0.999999 2.000000 0.9999 0.000000 0.000000 0.0000 0.000000 NO HA Ya convergió..
Ya convergió.. Ya lo convergido.. Ya lo tiene Ya convergió.. Ya
2 28 40 87 9972 00 02 0001 01 tiene pequeño continúe pequeño lo tiene pequeño
- -
48 3.000000 0.999999 2.000000 0.9999 0.000000 0.000000 0.0000 0.000000 NO HA Ya convergió..
Ya convergió.. Ya lo convergido.. Ya lo tiene Ya convergió.. Ya
3 27 41 85 9973 00 01 0001 01 tiene pequeño continúe pequeño lo tiene pequeño
- -
48 3.000000 0.999999 2.000000 0.9999 0.000000 0.000000 0.0000 0.000000 NO HA Ya convergió..
Ya convergió.. Ya lo convergido.. Ya lo tiene Ya convergió.. Ya
4 27 43 83 9974 00 01 0001 01 tiene pequeño continúe pequeño lo tiene pequeño
- -
48 3.000000 0.999999 2.000000 0.9999 0.000000 0.000000 0.0000 0.000000 NO HA Ya convergió..
Ya convergió.. Ya lo convergido.. Ya lo tiene Ya convergió.. Ya
5 26 44 81 9974 00 01 0001 01 tiene pequeño continúe pequeño lo tiene pequeño
- -
48 3.000000 0.999999 2.000000 0.9999 0.000000 0.000000 0.0000 0.000000 NO HA Ya convergió..
Ya convergió.. Ya lo convergido.. Ya lo tiene Ya convergió.. Ya
6 25 45 79 9975 00 01 0001 01 tiene pequeño continúe pequeño lo tiene pequeño
- -
48 3.000000 0.999999 2.000000 0.9999 0.000000 0.000000 0.0000 0.000000 NO HA Ya convergió..
Ya convergió.. Ya lo convergido.. Ya lo tiene Ya convergió.. Ya
7 25 47 77 9976 00 01 0001 01 tiene pequeño continúe pequeño lo tiene pequeño
- -
48 3.000000 0.999999 2.000000 0.9999 0.000000 0.000000 0.0000 0.000000 NO HA Ya convergió..
Ya convergió.. Ya lo convergido.. Ya lo tiene Ya convergió.. Ya
8 24 48 75 9976 00 01 0001 01 tiene pequeño continúe pequeño lo tiene pequeño
- -
48 3.000000 0.999999 2.000000 0.9999 0.000000 0.000000 0.0000 0.000000 NO HA Ya convergió..
Ya convergió.. Ya lo convergido.. Ya lo tiene Ya convergió.. Ya
9 23 49 73 9977 00 01 0001 01 tiene pequeño continúe pequeño lo tiene pequeño
- -
49 3.000000 0.999999 2.000000 0.9999 0.000000 0.000000 0.0000 0.000000 NO HA Ya convergió..
Ya convergió.. Ya lo convergido.. Ya lo tiene Ya convergió.. Ya
0 23 51 71 9977 00 01 0001 01 tiene pequeño continúe pequeño lo tiene pequeño
- -
49 3.000000 0.999999 2.000000 0.9999 0.000000 0.000000 0.0000 0.000000 NO HA Ya convergió..
Ya convergió.. Ya lo convergido.. Ya lo tiene Ya convergió.. Ya
1 22 52 69 9978 00 01 0001 01 tiene pequeño continúe pequeño lo tiene pequeño
- -
49 3.000000 0.999999 2.000000 0.9999 0.000000 0.000000 0.0000 0.000000 NO HA Ya convergió..
Ya convergió.. Ya lo convergido.. Ya lo tiene Ya convergió.. Ya
2 22 53 68 9979 00 01 0001 01 tiene pequeño continúe pequeño lo tiene pequeño
- -
49 3.000000 0.999999 2.000000 0.9999 0.000000 0.000000 0.0000 0.000000 NO HA Ya convergió..
Ya convergió.. Ya lo convergido.. Ya lo tiene Ya convergió.. Ya
3 21 54 66 9979 00 01 0001 01 tiene pequeño continúe pequeño lo tiene pequeño
- -
49 3.000000 0.999999 2.000000 0.9999 0.000000 0.000000 0.0000 0.000000 NO HA Ya convergió..
Ya convergió.. Ya lo convergido.. Ya lo tiene Ya convergió.. Ya
4 21 55 64 9980 00 01 0001 01 tiene pequeño continúe pequeño lo tiene pequeño
- -
49 3.000000 0.999999 2.000000 0.9999 0.000000 0.000000 0.0000 0.000000 NO HA Ya convergió..
Ya convergió.. Ya lo convergido.. Ya lo tiene Ya convergió.. Ya
5 20 56 63 9980 00 01 0001 01 tiene pequeño continúe pequeño lo tiene pequeño
- -
49 3.000000 0.999999 2.000000 0.9999 0.000000 0.000000 0.0000 0.000000 NO HA Ya convergió..
Ya convergió.. Ya lo convergido.. Ya lo tiene Ya convergió.. Ya
6 20 58 61 9981 00 01 0001 00 tiene pequeño continúe pequeño lo tiene pequeño
- -
49 3.000000 0.999999 2.000000 0.9999 0.000000 0.000000 0.0000 0.000000 NO HA Ya convergió..
Ya convergió.. Ya lo convergido.. Ya lo tiene Ya convergió.. Ya
7 19 59 60 9981 00 01 0001 00 tiene pequeño continúe pequeño lo tiene pequeño

pág. 34
 - -
49 3.000000 0.999999 2.000000 0.9999 0.000000 0.000000 0.0000 0.000000 NO HA Ya convergió..
Ya convergió.. Ya lo convergido.. Ya lo tiene Ya convergió.. Ya
8 19 60 58 9982 00 01 0001 00 tiene pequeño continúe pequeño lo tiene pequeño
- -
49 3.000000 0.999999 2.000000 0.9999 0.000000 0.000000 0.0000 0.000000 NO HA Ya convergió..
Ya convergió.. Ya lo convergido.. Ya lo tiene Ya convergió.. Ya
9 18 61 57 9982 00 01 0001 00 tiene pequeño continúe pequeño lo tiene pequeño
- -
50 3.000000 0.999999 2.000000 0.9999 0.000000 0.000000 0.0000 0.000000 Ya convergió.. Ya convergió..
Ya convergió.. Ya lo Ya lo tiene Ya lo tiene Ya convergió.. Ya
0 18 62 55 9982 00 01 0001 00 tiene pequeño pequeño pequeño lo tiene pequeño

Resultados:
𝒙𝟏 = 𝟑. 𝟎𝟎𝟎𝟎𝟎𝟎𝟏𝟕𝟖𝟏𝟒𝟒𝟎𝟏 ≅ 𝟑
𝒙𝟐 = 𝟎. 𝟗𝟗𝟗𝟗𝟗𝟗𝟔𝟏𝟔𝟐𝟓𝟎𝟔𝟒𝟔 ≅ 𝟏
𝒙𝟑 = −𝟐. 𝟎𝟎𝟎𝟎𝟎𝟎𝟓𝟓𝟑𝟖𝟒𝟗𝟗𝟒 ≅ −𝟐
𝒙𝟒 = 𝟎. 𝟗𝟗𝟗𝟗𝟗𝟗𝟖𝟐𝟒𝟓𝟐𝟔𝟓𝟐𝟓 ≅ 𝟏

pág. 35
Realice una gráfica de la forma como va convergiendo la solución (Número
de iteraciones vs norma del error). Realice una breve explicación,
sustentándose en la teoría revisada, acerca de los resultados. ¿Cuál
considera que es el mejor método para este SEL en particular? ¿por qué?

En la grafica observamos que el método jacobi no converge; el método gauss seidel

(475 iteraciones) y gauss seidel relajado (499 iteraciones a 0.99 de factor de
ponderación) si lo hacen.

Si bien es cierto que el método relajado tiene más iteraciones y lo conforma más
partes de cálculo, es el más preciso ya que se trata de usar técnicas estadísticas
para ir ajustando en la marcha el valor encontrado y asi ir consiguiendo mayor
precisión.

pág. 36
Aporte 2. Interpolación Numérica

Ejercicio 5. Para estudiantes cuyo documento de identidad termine en

dígito impar:

Considere la siguiente tabla de una función f(x) que describe la cantidad de

personas que entran/salen de un elevador en cada piso x en una subida
entre el piso 2 y el piso 2:

x -2 -1 0 1 2
f(x) a 2 1 0 4

Usando el método de interpolación de Newton, determine el número de

personas (a) que entran/salen del piso -2, de tal forma que el polinomio
interpolante de f(x) en los puntos de la tabla dada sea de grado 3.
Justifique su respuesta según lo revisado en la teoría.

Para resolver el problema, vale la pena primero citar de forma abstracta la siguiente
imagen para su resolución a tercer grado.

(Chapra, 2007, p. 509)

pág. 37
Ahora bien, para la interpolación de newton, necesitamos recordar lo que es una
interpolación cuadrática, en donde sus teoremas están dados por las siguientes
ecuaciones:

𝑏0 = 𝑓(𝑥0 )

𝑓(𝑥1 ) − 𝑓(𝑥0 )
𝑏1 =
𝑥1 − 𝑥0

𝑓(𝑥2 ) − 𝑓(𝑥1 ) 𝑓(𝑥1 ) − 𝑓(𝑥0 )

𝑥2 − 𝑥1 − 𝑥1 − 𝑥0
𝑏2 =
𝑥2 − 𝑥0

𝑓(𝑥3 ) − 𝑓(𝑥2 ) 𝑓(𝑥2 ) − 𝑓(𝑥1 ) 𝑓(𝑥2 ) − 𝑓(𝑥1 ) 𝑓(𝑥1 ) − 𝑓(𝑥0 )

− −
𝑥3 − 𝑥2 𝑥2 − 𝑥1 𝑥2 − 𝑥1 𝑥1 − 𝑥0
𝑥3 − 𝑥1 − 𝑥2 − 𝑥0
𝑏3 =
𝑥3 − 𝑥0

Etc….

Ahora bien, la ecuación cuadrática nos habla que se tienen tres puntos como datos,
éstos pueden ajustarse en un polinomio de segundo grado (también conocido como
polinomio cuadrático o parábola) para hallar un valor. (Chapra, 2007, p. 506)

Para hallar f2(x)

𝑓2 (𝑥) = 𝑏0 + 𝑏1 (𝑥 − 𝑥0 ) + 𝑏2 (𝑥 − 𝑥0 )(𝑥 − 𝑥1 )

Para hallar f3(x)

𝑓3 (𝑥) = 𝑏0 + 𝑏1 (𝑥 − 𝑥0 ) + 𝑏2 (𝑥 − 𝑥0 )(𝑥 − 𝑥1 ) + 𝑏3 (𝑥 − 𝑥0 )(𝑥 − 𝑥1 )(𝑥 − 𝑥2 )

Etc….

Ahora bien, para ensamblar a una interpolación de newton, tenemos:

𝑏0 = 𝑓(𝑥0 )
𝑏1 = 𝑓[𝑥1 , 𝑥0 ]
𝑏2 = 𝑓[𝑥2 , 𝑥1 , 𝑥0 ]
𝑏3 = 𝑓[𝑥3 , 𝑥2 , 𝑥1 , 𝑥0 ]

Etc…

Si observamos la representación gráfica, lo que debemos hacer es seguir las

diferencias divididas finitas.

pág. 38
Dicho todo esto, entramos en materia:

El ejercicio es el siguiente y lo que hay que hacer es hallar a.

x -2 -1 0 1 2
f(x) a 2 1 0 4

Preparamos la tabla de iteraciones con los valores completos para proceder a tercer
Grado.

Iteración x f(x)
0 -1 2
1 0 1
2 1 0
3 2 4

Hallamos el 1er grado de la secuencia que tenemos:

𝑓(𝑥1 ) − 𝑓(𝑥0 )
𝑏1 =
𝑥1 − 𝑥0

Donde
𝑏1 = 𝑓[𝑥1 , 𝑥0 ]

1−2 −1
𝑓[𝑥1 , 𝑥0 ] = = = −1
0 − (−1) 1

0 − 1 −1
𝑓[𝑥2 , 𝑥1 ] = = = −1
1−0 1

4−0 4
𝑓[𝑥3 , 𝑥2 ] = = =4
2−1 1

pág. 39
Hallamos el 2do grado de la secuencia que tenemos:

𝑓(𝑥2 ) − 𝑓(𝑥1 ) 𝑓(𝑥1 ) − 𝑓(𝑥0 )

𝑥2 − 𝑥1 − 𝑥1 − 𝑥0
𝑏2 =
𝑥2 − 𝑥0

Donde:
𝑏2 = 𝑓[𝑥2 , 𝑥1 , 𝑥0 ]

Esto puede re-expresarse de la siguiente forma:

𝑓(𝑥2 ) − 𝑓(𝑥1 ) 𝑓(𝑥1 ) − 𝑓(𝑥0 )

𝑥2 − 𝑥1 − 𝑥1 − 𝑥0
]
𝑓[𝑥2 , 𝑥1 , 𝑥0 =
𝑥2 − 𝑥0

𝑓[𝑥2 , 𝑥1 ] − 𝑓[𝑥1 , 𝑥0 ]
𝑓[𝑥2 , 𝑥1 , 𝑥0 ] =
𝑥2 − 𝑥0

Nota: Si ya lo notaron, el número de ecuaciones resultantes está relacionado

con el número de datos que tenemos… para explicar la formula, utilizo solo los
primeros términos, sin embargo, cuando paso a la siguiente ecuación, aplico
la misma fórmula pero con los datos del rango correspondiente que deseo
hallar.

Aplicaremos para el segundo grado:

𝑓[𝑥2 , 𝑥1 ] − 𝑓[𝑥1 , 𝑥0 ]
𝑓[𝑥2 , 𝑥1 , 𝑥0 ] =
𝑥2 − 𝑥0

𝑓[𝑥3 , 𝑥2 ] − 𝑓[𝑥2 , 𝑥1 ]
𝑓[𝑥3 , 𝑥2 , 𝑥1 ] =
𝑥3 − 𝑥1

Procedamos:

−1 − (−1) −1 + 1 0
𝑓[𝑥2 , 𝑥1 , 𝑥0 ] = = = =0
1 − (−1) 1+1 2

4 − (−1) 5
𝑓[𝑥3 , 𝑥2 , 𝑥1 ] = = = 2.5
2−0 2

pág. 40
Hallamos el 3er grado de la secuencia que tenemos:

𝑓(𝑥3 ) − 𝑓(𝑥2 ) 𝑓(𝑥2 ) − 𝑓(𝑥1 ) 𝑓(𝑥2 ) − 𝑓(𝑥1 ) 𝑓(𝑥1 ) − 𝑓(𝑥0 )

𝑥3 − 𝑥2 − 𝑥2 − 𝑥1 𝑥2 − 𝑥1 − 𝑥1 − 𝑥0
𝑥3 − 𝑥1 − 𝑥2 − 𝑥0
𝑏3 =
𝑥3 − 𝑥0

𝑏3 = 𝑓[𝑥3 , 𝑥2 , 𝑥1 , 𝑥0 ]

Esto puede re-expresarse de la siguiente forma:

𝑓(𝑥3 ) − 𝑓(𝑥2 ) 𝑓(𝑥2 ) − 𝑓(𝑥1 ) 𝑓(𝑥2 ) − 𝑓(𝑥1 ) 𝑓(𝑥1 ) − 𝑓(𝑥0 )

− −
𝑥3 − 𝑥2 𝑥2 − 𝑥1 𝑥2 − 𝑥1 𝑥1 − 𝑥0
𝑥3 − 𝑥1 − 𝑥2 − 𝑥0
𝑓[𝑥3 , 𝑥2 , 𝑥1 , 𝑥0 ] =
𝑥3 − 𝑥0

𝑓[𝑥3 , 𝑥2 , 𝑥1 ] − 𝑓[𝑥2 , 𝑥1 , 𝑥0 ]
𝑓[𝑥3 , 𝑥2 , 𝑥1 , 𝑥0 ] =
𝑥3 − 𝑥0

2.5 − 0 2.5
𝑓[𝑥3 , 𝑥2 , 𝑥1 , 𝑥0 ] = = 0.833333333333
2 − (−1) 3

Resumen hasta aquí:

Iteración x f(x) grado1 grado2 grado3

0 -1 2
1 0 1 -1
2 1 0 -1 0
3 2 4 4 2.5 0.83333333

𝑓[𝑥1 , 𝑥0 ] = −1
𝑓[𝑥2 , 𝑥1 ] = −1
𝑓[𝑥3 , 𝑥2 ] = 4
𝑓[𝑥2 , 𝑥1 , 𝑥0 ] = 0
𝑓[𝑥3 , 𝑥2 , 𝑥1 ] = 2.5
𝑓[𝑥3 , 𝑥2 , 𝑥1 , 𝑥0 ] = 0.833333333333

pág. 41
Teniendo estos datos, podemos proceder con la fórmula para hallar f(x) a tercer
grado, la cual es:

𝑓3 (𝑥) = 𝑏0 + 𝑏1 (𝑥 − 𝑥0 ) + 𝑏2 (𝑥 − 𝑥0 )(𝑥 − 𝑥1 ) + 𝑏3 (𝑥 − 𝑥0 )(𝑥 − 𝑥1 )(𝑥 − 𝑥2 )

𝑓3 (𝑥) = 2 + (−1)(𝑥 − (−1)) + 0(𝑥 − (−1))(𝑥 − 0) + 0.83333333(𝑥 − (−1))(𝑥 − 0)(𝑥 − 1)

𝑓3 (𝑥) = 2 + (−1)(𝑥 + 1) + 0(𝑥 + 1)(𝑥 − 0) + 0.83333333(𝑥 + 1)(𝑥 − 0)(𝑥 − 1)

¡AHORA SI…. PROCEDAMOS A HALLAR EL VALOR DE a!

El valor oculto de (a) corresponde cuando x = -2… reemplacemos.

𝑓3 (−2) = 2 + (−1)((−2) + 1) + 0((−2) + 1)((−2) − 0) + 0.83333333((−2) + 1)((−2) − 0)((−2)

− 1)

𝑓3 (−2) = 2 + (−1)(−1) + 0(−1)(−2) + 0.83333333(−1)(−2)(−3)

𝑓3 (−2) = 2 + 1 + 0 + 0.83333333(−6)

𝑓3 (−2) = 3 − 4.99999999998

𝑓3 (−2) = −1.99999999998 ≅ −2

Resultado:
𝒄𝒖𝒂𝒏𝒅𝒐 𝒙 = −𝟐, 𝒇(𝒙) = 𝒂 = −𝟏. 𝟗𝟗𝟗𝟗𝟗𝟗𝟗𝟗𝟗𝟗𝟖

x -2 -1 0 1 2
f(x) -1.99998 2 1 0 4

pág. 42
IMÁGENES DE REALIZACION EN EXCEL

pág. 43
pág. 44
 CONCLUSIONES

Los métodos iterativos constituyen una forma MACRO de resolver sistemas de

ecuaciones con varias incógnitas. Tanto los métodos directos como los iterativos
conducen a la misma vía, sin embargo, el manejo de grandes volúmenes implica la
aplicación de técnicas más mecánicas y aquí es donde entran las matrices
iterativas; es preciso resaltar que en ocasiones algún método iterativo puede no
funcionar, es decir, es posible que se pueda por gauss-seidel y no jacobi, o puede
suceder lo contrario; también puede suceder que se pueda con ambos o
simplemente, que no se pueda con ninguno (según material de estudio existen
casos en donde no se pueda.)

La interpolación de newton es una extensión de la interpolación cuadrática; el

concepto es sencillo; basados en una serie de puntos conocidos, se arma y obtiene
una serie de constantes y esas constantes se ensamblan a una formula y esa
fórmula constituye la regla para ese polinomio en donde simplemente se le dice el
valor en posición x y así sabemos f(x).

pág. 45
 REFERENCIAS

Álvarez, C. (09 de julio de 2019). Métodos Numéricos (100401) – 803 - Segunda

Webconferencia: Sistema de Ecuaciones Lineales e Interpolación Numérica.
[Archivo de Video]. Recuperado de
http://conferencia2.unad.edu.co/pu3g652djk5w/

Zuñiga, C. (2010). Módulo Algebra Lineal. Bogotá, Colombia: UNAD. Recuperado

de: http://hdl.handle.net/10596/7193

Chapra, E., Canale, R. (Ed). (2007). Métodos numéricos para ingenieros Quinta
Edición. Ciudad de México D.F., México: McGraw-Hill. ISBN 970-10-6114-4.
Recuperado de http://artemisa.unicauca.edu.co/~cardila/Chapra.pdf

Nieves, A., Domínguez, F. (2014). Métodos Numéricos aplicados a la ingeniera 4ª.

Edición. Ciudad de México D.F., México: Editorial Patria. ISBN 978-607-438-
926-5. Recuperado de https://ebookcentral-proquest-
com.bibliotecavirtual.unad.edu.co/lib/unadsp/reader.action?docID=3227640&
ppg=48

Garcia, A. [ALVARO GARCIA SEVERICHE]. (23 de noviembre de 2015). metodo

Gauss seidel bien explicado. [Archivo de Video]. Recuperado de
https://www.youtube.com/watch?v=GeW0WqiWZbA

Espejo, D. [Diana Espejero]. (06 de noviembre de 2013). Método Iterativo Para La

Solución De Sistemas De Ecuaciones Lineales (Jacobi). [Archivo de Video].
Recuperado de https://www.youtube.com/watch?v=7Dqf9vird1w

pág. 46