REPÚBLICA BOLIVARIANA DE VENEZUELA.

MINISTERIO DEL PODER POPULAR PARA LA EDUCACIÓN.

I.U.P. “SANTIAGO MARIÑO”

VALENCIA, EDO-CARABOBO

METODO DE GAUSS SEIDEL

ALUMNO:

José Nabil Nassr De Sousa

CI. 21.238.549

Esc#44

Valencia, marzo 2021

Método de Gauss Seidel

Este es uno de los métodos más interesantes del análisis numérico y

particularmente útil ya que nos permite encontrar la solución de un sistema de “n”
ecuaciones con “n” incógnitas.

Para comenzar es preciso mencionar que es un método iterativo, es decir que

debe aplicarse recursivamente hasta encontrar una solución adecuada o con un error
considerablemente pequeño.

En cada iteración obtenemos una solución posible del sistema con un error
determinado, a medida que aplicamos nuevamente el método, la solución puede ser
más precisa, entonces se dice que el sistema converge, pero si al aplicar el método
reiteradas veces la solución tiene un error cada vez mayor se dice que el sistema no
converge y no se puede resolver el sistema de ecuaciones por este método.

Teniendo el siguiente sistema de ecuaciones:

Despejamos x1 de la ecuación 1, x2 de la ecuación 2, … , xn de la ecuación n,

quedando:

Ahora despejamos x2 de la ecuación 2 y reemplazamos a x1 por el valor

obtenido en la ecuación anterior. De esto nos queda:
 Una vez que tenemos x2, despejamos x3 de la ecuación 3 y así sucesivamente
con las n ecuaciones, cada vez asignando el valor de las x1, x2, … xn-1 obtenido en el
paso anterior. Cuando hemos despejado las xn, tenemos lo que se conoce como
primera solución o solución de la primera iteración:

Con los nuevos valores de x1, x2, …, xn aplicamos los mismos pasos anteriores,
pero con los nuevos valores de las xn, de esta manera conseguimos una segunda
solución:

Al tener esta segunda solución estamos en condiciones de calcular el error que

se calcula como sigue:
 Así, repetimos el método tantas veces hasta que el error sea muy pequeño o los
suficientemente aceptable.

Ahora solo queda mencionar que para que un sistema sea convergente se debe
cumplir que la matriz de coeficientes sea diagonalmente dominante, y para ello se debe
verificar la siguiente expresión:

Si no se cumple esa condición, se puede permutar las filas de la matriz, con el

fin de poder convertirla en una diagonalmente dominante.

Este método es prácticamente idéntico al de Jacob, la única diferencia estriba en

que el método de Gauss-Seidel se acerca más rápido a la solución cuando el método
converge, 4 debido a que una vez que se calcula la componente ( 1) k i x  la utiliza
inmediatamente en la misma iteración, esto es:
El criterio de convergencia de este método es el mismo que el de Jacobi.

Ejercicios propuestos:

a) X -7Y + 2Z = 19
2X + Y + 9Z = 13            
8X + Y – Z = 21

Reordenamos las ecuaciones

8X + Y – Z = 21
X -7Y + 2Z = 19
2X + Y + 9Z = 13
Despejamos las ecuaciones y queda
21− y + z
x=
8
19− x−2 z
y=
−7
 13−2 x− y
z=
9

Valores iniciales ( x , y , z )=( 0,0,0)

Iteración 1=

21−0+0
x 1= =2.625
8
19−2.625+ 2∗0
y 1= =−2.33929
−7
13−2∗( 2.625 )−(−2.33929)
z 1= =1.1210 3
9
Iteración 2=

21−(−2.33929)+1.12103
x 2= =3.05754
8
19−3.05754 +2∗1.12103
y 2= =−1.95720
−7
13−2∗( 3.05754 )−(−1.95720)
z 2= =0.98246
9
Iteración 3=

21−(−1.95720)+0.98246
x 3= =2.99246
8
19−2.99246+2∗0.98246
y 3= =−2.00609
−7
13−2∗( 2.99246 )−(−2.00609)
z 3= =1.12103
9
Iteración 4=

21−(−2.00609)+1.12103
x4 = =3.00106
8
 19−3.00106+2∗1.12103
y4= =−1.99918
−7
13−2∗( 3.00106 )−(−1.99918)
z4 = =0.99967
9

Iteración 5=

21−(−1.99918)+0.99967
x 5= =2.99986
8
19−2.99986+2∗0.99967
y 5= =−2.00011
−7
13−2∗( 2.99986 )−(−2.00011)
z 5= =1.0000 4
9
Iteración 6=

21−(−2.00011)+1.00004
x 6= =3.00002
8
19−3.00002+2∗1.00004
y 6= =−1.99998
−7
13−2∗( 3.00002 )−(−1.99998)
z 6= =0.99999
9
Iteración 7=

21−(−1.99998)+ 0.99999
x 7= =3.00000
8
19−3.00000+ 2∗0.99999
y 7= =−2.00000
−7
13−2∗( 3.00000 )−(−2.00000)
z 7= =1.00000
9
Iteración 8=

21−(−2.00000)+1.00000
x 8= =3.00000
8
19−3.00000+ 2∗1.00000
y 8= =−2.00000
−7
 13−2∗( 3.00000 )−(−2.00000)
z 8= =1.00000
9

Iteración 9=

21−(−2.00000)+1.00000
x 9= =3.00000
8
19−3.00000+ 2∗1.00000
y 9= =−2.00000
−7
13−2∗( 3.00000 )−(−2.00000)
z 9= =1.0000 0
9

iteració
n x y z Err x err y err z
0 0 0 0      
1 2.625 -2.33929 1.12103 100.00000 -100.00000 100.00000
2 3.05754 -1.95720 0.98246 14.14666 -19.52209 14.10482
3 2.99246 -2.00609 1.00235 2.17489 -2.43708 1.98483
4 3.00106 -1.99918 0.99967 0.28650 -0.34577 0.26796
5 2.99986 -2.00011 1.00004 0.03997 -0.04683 0.03705
6 3.00002 -1.99998 0.99999 0.00545 -0.00646 0.00507
7 3.00000 -2.00000 1.00000 0.00075 -0.00088 0.00070
8 3.00000 -2.00000 1.00000 0.00010 -0.00012 0.00010
9 3.00000 -2.00000 1.00000 0.00001 -0.00002 0.00001

x=3
y=−2
z=1

b) −2 X +10 Y +3 Z +2 W =35
−X +3 Y −2 Z +9 W =35
−8 X +Y + Z +5 W =17
3 X – 4 Y +7 Z+W =2 0
 Reordenamos las ecuaciones
9 W −X +3 Y −2 Z=3 5
5 W−8 X +Y + Z=17
2 W −2 X +10 Y +3 Z=35
W +3 X – 4 Y +7 Z=2 0

Despejamos las ecuaciones y queda

35+ x−3 y +2 z
w=
9
17−5 w− y −z
x=
−8
35−2 w+2 x−3 z
y=
10
20−w−3 x + 4 y
z=
7

Iteración 1=

35+0−3(0)+2(0)
w 1= =3.88889
9
17−5 ( 3.88889 )−0−0
x 1= =0.30556
−8
35−2 ( 3.88889 ) +2(0.30556)−3(0)
y 1= =2.78333
−7
20−3.88889−3 ( 0.30556 ) +4 (2.78333)
z 1= =3.76111
9

Iteración 2=

35+0.30556−3(2.78333)+ 2(3.76111)
w 2= =3.83086
9
17−5 ( 3.83086 ) −0.30556−3.76111
x 2= =1.08735
−8
35−2 ( 3.83086 ) +2(1.08735)−3 (3.76111)
y 2= =1.82296
−7
 20−3.83086−3 ( 1.08735 ) +4 (1.82296)
z 2= =2.88556
9
Iteración 3=
35+1.08735−3 (1.82296)+2( 2.88556)
w 3= =4.04329
9
17−5 ( 4.04329 )−1.82296−2.88556
x 3= =0.99062
−8
35−2 ( 4.04329 )+ 2(0.99062)−3(2.88556)
y 3= =2.02380
−7
20−4.04329−3 ( 0.99062 ) + 4(2.02380)
z 3= =3.01143
9

Iteración 4=

35+0.99062−3(2.02380)+2(3.01143)
w 4= =3.99357
9
17−5 ( 3.99357 )−2.02380−3.01143
x4 = =1.00038
−8
35−2 ( 3.99357 ) +2(1.00038)−3(3.01143)
y4= =1.99793
−7
20−3.99357−3 ( 1.00038 ) + 4(1.99793)
z4 = =2.99957
9
Iteración 5=
35+1.00038−3 (1.99793)+2(2.99957)
w 5= =4.00064
9
17−5 ( 4.00064 )−1.99793−2.99957
x 5= =1.00009
−8
35−2 ( 4.00064 ) +2(1.00009)−3(2.99957)
y 5= =2.00002
−7
20−4.00064−3 (1.00009 )+ 4 (2.00002)
z 5= =2.99988
9
Iteración 6=
 35+1.00009−3 (2.00002)+2(2.99988)
w 6= =3.99998
9
17−5 ( 3.99998 ) −2.00002−2.99988
x 6= =0.99997
−8
35−2 ( 3.99998 ) +2( 0.99997)−3(2.99988)
y 6= =2.00003
−7
20−3.99998−3 ( 0.99997 ) + 4(2.00003)
z 6= =3.00003
9
Iteración 7=
35+0.99997−3(2.00003)+2(3.00003)
w 7= =3.99999
9
17−5 ( 3.99999 ) −2.00003−3.00003
x 7= =1.00000
−8
35−2 ( 3.99999 ) +2(1.00000)−3 (3.00003)
y 7= =1.99999
−7
20−3.99999−3 ( 1.00000 ) +4 (1.99999)
z 7= =2.99999
9

Iteración 8=

35+1.00000−3 (1.99999)+2( 2.99999)

w 8= =4.00000
9
17−5 ( 4.00000 )−1.99999−2.99999
x 8= =1.00000
−8
35−2 ( 4.00000 ) +2(1.00000)−3(2.99999)
y 8= =2.00000
−7
20−4.00000−3 ( 1.00000 ) + 4(2.00000)
z 8= =3.00000
9
iteració
n w x y z Err x err y err z err w
0 0 0 0 0        
1 3.88889 0.30556 2.78333 3.76111 100.00000 100.00000 100.00000 100.00000
2 3.83086 1.08735 1.82296 2.88556 1.51466 71.89895 52.68184 30.34230
3 4.04329 0.99062 2.02380 3.01143 5.25372 9.76410 9.92364 4.17974
4 3.99357 1.00038 1.99793 2.99957 1.24503 0.97587 1.29455 0.39541
5 4.00064 1.00009 2.00002 2.99988 0.17673 0.02965 0.10424 0.01028
6 3.99998 0.99997 2.00003 3.00003 0.01648 0.01130 0.00084 0.00507
7 3.99999 1.00000 1.99999 2.99999 0.00039 0.00309 0.00213 0.00133
8 4.00000 1.00000 2.00000 3.00000 0.00022 0.00048 0.00046 0.00020

w=4
x=1
y=2
z=3