lOMoARcPSD|17818758

lO

M SISTEMA DE ECUACIONES LINEALES Y NO LINEALES

ACTIVIDAD EVALUATIVA EJE 2

CRISTIAN ANDRES ANTONIO ESTEPA

DANILO DE JESÚS ARIZA AGÁMEZ

FUNDACION UNIVERSITARIA DEL AREA ANDINA

INGENIERIA DE SISTEMAS
ANALISIS NUMERICO
2023
 lOMoARcPSD|17818758

Taller Sistemas de Ecuaciones Lineales y No Lineales

1. Resolver el siguiente sistema de ecuaciones lineales por el método gráfico y

determinantes.

1 22 1
x+2 y= x − y=2
2 33 3

2. Aplicar la regla de Cramer para demostrar que el x=1, y=2 y z=3 son las soluciones

para el sistema dado por:

2 y−z=1−x + y =12 x− y + z=3

3. Mediante la eliminación de incógnitas, resolver:

4 x + 4 y=−10 8 x−5 y=5

 lOMoARcPSD|17818758

4. Considere el sistema propuesto

A = [ 1−35−3−27 8 4 9 ]

Encuentra la solución del sistema asociado aplicando la factorización LU.

5. Utilice el método de Jacobi para encontrar, a partir de cuatro iteraciones, la solución

a la matriz asociada a un sistema de ecuaciones lineales dada por:

A= [ 0.2 0.1−0.5 0.7 ] ,con b=[ 0.75 0.35 ] y x 0= [ 1 1 ]

6. Aplicando el método de Gauss-Seidel, resolver:

2 x1 +2 x2 −3 x3=2−4 x2−4 x3=6 x1 +3 x2 −10 x3=3

 lOMoARcPSD|17818758

Desarrollo

1. Resolver el siguiente sistema de ecuaciones lineales por el método gráfico y

determinantes.

1 22 1
x+2 y= x − y=2
2 33 3

Se ordena para mayor claridad

Método del determinante

1 2

[ ]
x+2 y=
2 3 1 2 2
2 x− 1 2 3
y=2 2 −1
3 2
3 3 3

[ ]
2 1

[ ]
2 3 2
1 2 ∣D∣= 2
2 −1 2 −1
2
3 3 3 3
 lOMoARcPSD|17818758

∣D∣
1−1
2 )( ) (
−(2
)
2
3

(=
3
)
¿− 1 −
4
63
¿− 3
2

[ ]
2
2

[ ] ( )(
2 3 3
1 2 ∣Dx∣=
2 −1
2 −1
3 3 2 3

)
2 −1
∣Dx∣= −(2) (2)
3 3
2
¿− −4
9
−38
¿ 9

[ ] [ ]
1 2 2 1 2
2 3 ∣Dy∣= 2 3
2 −1 2
2
3 3 3 2

∣Dy∣= ( ) ( )( )
2
1
(2)−
2 2
3 3
2 4
¿ 2 −9
5
¿9
 lOMoARcPSD|17818758

−38D 5
9 76 D 9 −10
X= = = =2.8148
x Y= y
= = =−0.37037
D −3 27 D −3 27
2 2

Método gráfico:

Tanto el método de los determinantes como el método gráfico nos da el mismo punto de

corte en la gráfica
 lOMoARcPSD|17818758

2. Aplicar la regla de Cramer para demostrar que el x=1, y=2 y z=3 son las soluciones

para el sistema dado por:

2 y−z=1−x + y =12 x− y + z=3

Se ordena para mayor claridad

Método de Cramer

0 x+ 2 y −z=1
−x+ y+ 0
z=1 2 x −
y+ z=3
0 2 −1

[
∣D∣= −1 1
2 −1
0
1 ]
∣D∣=0 ( 1∗1−0 (−1) ) +2 ( 0∗2−(−1)∗1 ) +(−1)((−1)∗(−1)−1∗2)
∣D∣=0 (1−0)+2 (0+1) +(−1)(1−2)
∣D∣=0+2+1
∣D∣=3

1 2 −1

∣DX∣= 1 1
3 −1
[ 0
1 ]
∣D X∣=1 ( 1∗1−0 (−1) ) +2 ( 0∗3−1∗1)+(−1)(1∗(−1)−1∗3)
∣D X∣=1(1−0)+ 2(0−1)+(−1)(−1−3)
∣D X∣=1−2+ 4
∣DX∣=3
0 1 1

∣Dy∣= −1
[ ]
2
1 0
3 1

∣D y∣=0 (1∗1−0∗3 )+1 ( 0∗2−(−1)∗1 ) +1(−1∗3−1∗2)

∣Dy∣=0+1 (0+1)+1(−3−2)
∣Dy∣=1−5
 lOMoARcPSD|17818758

∣Dy∣=−4

0 2 1

∣Dz∣=
[−1 1 1
2 −1 3 ]
∣Dz∣=0 ( 1∗3−1∗(−1) ) +2 ( 1∗2−(−1)∗3 ) +1((−1)∗(−1)−1∗2)
∣Dz∣=0+ 2(2+3 )+1 (1−2)
∣Dz∣=10−1
∣D z∣=9

∣D∣=9∣Dx∣=3∣Dy∣=−4∣Dz∣=9

∣ D x ∣ ∣ D y ∣ ∣ D z∣
x= ∣ ∣ y=∣ ∣ z=∣ ∣
D D D
3 −4 9
x= y = z=
9 9 9

x=0,333 y=−0,444 z=1

 lOMoARcPSD|17818758

3. Mediante la eliminación de incógnitas, resolver:

4 x + 4 y=−10 8 x−5 y=5

Método de eliminación

4x+ 4
y=−1

Eliminamos una incógnita (en este caso la

x)

4x+ 4 *(-2)
y=−10
8 x−5 y=5

-8x - 8y = 20
8x – 5y = 5
-13y = 25

−13 y=25

−25
y= 13

y=−1,9230

Reemplazamos la y en cualquiera de las ecuaciones (lo hacemos en ambas para

comprobar)
4 x + 4 y=−10 8 x−5 y=5

4x+ 4 ( )
−25 8 x−5 ( 13 )=5
−25

=−10
13 125
8 x+ =5
100 13
4x− =−10
13 125
8 x=5−
100 13
4 x =−10+
13
−60
8 x= 13
−30
4 x = 13
 lOMoARcPSD|17818758

−30 −60
− 13 13
x= 4 x=
8
−15 −15
x= 26 x= 26

x=−0,57692307 x=−0,57692307

Comprobamos por método gráfico:

 lOMoARcPSD|17818758

4. Considere el sistema propuesto

A = [ 1−35−3−27 8 4 9 ]

Encuentra la solución del sistema asociado aplicando la factorización LU.

(
1 −3 5
−3 −2 7
8 4 9

Descomponemos para obtener las matrices correspondientes para L y U

A = L * U
1 −3 5 1 0 0 U 11 U 12 U 13

( −83 −42 79
) ( L231 01
0 ) ( U0 2
2
U 23 )
L132
1 5 1 5
−3 −3

( − −2 7
4 9 ) L2 –(-3L1)
(
0 −11 22
8 4 9 )
3
18 5 1

)
−3 −3 5

( 0
81
−11 22
4
−3
95 ) L3-8L1
(0
01
−11 22
28 −31
−3
5

( ) ( )
−28
L− L
0
0
−11 22
28 −31 3 ( 11 2 ) 0
0
−11 22
0 25

Termino convertido Factor usado en la Factor a ubicar

a cero conversión (signo contrario) Posición a ubicar
-3 3 -3 L21
8 -8 8 L31
28 28/11 -28/11 L32

A = L * U
1 −3 5 1 0 0 1 −3 5

( −83 −42 79
) ( −83
−2
81/
11
01
)
 lOMoARcPSD|17818758

( 0
−011
252
)
 lOMoARcPSD|17818758

6. Aplicando el método de Gauss-Seidel, resolver:

2 x1 +2 x2 −3 x3=2−4 x2−4 x3=6 x1 +3 x2 −10 x3=3

2=2 x1+2 x2−3 x3

6=−4 x2−4 x3
23=x1 +3 x2−10 x3

Primera Iteración
Para x1
b −a ∗x k −a ∗x k
(k+1) 1 12 2 13 3
x1 =
a11

2−( 2) (0)−(−3)(0)
x 1=
2

x1=1
Para x2
b −a ∗x k −a ∗x k
(k+1) 2 21 1 23 3
x2 =
a22

6 − ( 0 ) ( 1 ) − (− 4 ) ( 0 )
x 2= −4
x2=−1,5
Para x3
b −a ∗x k−a ∗x k
(k+1) 3 31 1 32 2
x3 =
a33

23−( 1 ) ( 1 )−( 3) ( −1,5)

x 3= −10
x3=−2.65
 lOMoARcPSD|17818758

Segunda Iteración
Para x1
b −a ∗x k −a ∗x k
(k+1) 1 12 2 13 3
x1 =
a11

2 −( 2 ) (− 1.5 )−(−3 ) (−2.65 )

x 1=
2

x1=−1.475
Para x2
b −a ∗x k −a ∗x k
(k+1) 2 21 1 23 3
x2 =
a22

6 −( 0 ) (−1.475 )−(−4 ) (−2.65 )

x 2= −4
x2=1.15
Para x3
b −a ∗x k−a ∗x k
(k+1) 3 31 1 32 2
x3 =
a33

23−( 1 ) (−1.475 )−( 3 )( 1.15 )

x 3= −10
x3=−2.1025

Como el ejercicio no nos indica hasta que % llevarlo o cuantas iteraciones hacer, pasamos
a reemplazar en la formula inicial

2=2 x1+2 x2−3 x3 6=−4 x2−4 x3

2=2 (−1.475 )+2 (1.15)−3 (−2.1025) 6=0(−1.475)−4 (1.15)−4 (−2.1025 )
2=5.6575 6=3.81
23=x1 +3 x2−10 x3
23=(−1.475 )+ 3 (1.15 )−10 (−2.1025)
23=23
 lOMoARcPSD|17818758

Tercera Iteración
Para x1
b −a ∗x k −a ∗x k
(k+1) 1 12 2 13 3
x1 =
a11

2 −( 2 ) (1.15 )−( −3) ( −2.1025)

x 1=
2

x1=−3.30375
Para x2
b −a ∗x k −a ∗x k
(k+1) 2 21 1 23 3
x2 =
a22

6 −( 0 )(−3.30375 )−(−4 )
x= (−2.1025)
2
−4
x2=0.6025
Para x3
b −a ∗x k−a ∗x k
(k+1) 3 31 1 32 2
x3 =
a33

23−( 1 ) (−3.30375 )−(3 )( 0.6025 )

x =3 −10
x3=−2.449625

Reemplazamos

2=2 x1+2 x2−3 x3

2=2 (−3.30375 )+2 (0.6025)−3 (−2.449625 )
2=1.946375
6=−4 x2−4 x3
6=0(−3.30375)−4 (0.6025)−4 (−2.449625 )
6=7.3885
23=x1 +3 x2−10 x3
23=(−3.30375)+ 3 (0.6025)−10 (−2.449625)
23=23
 lOMoARcPSD|17818758

Cuarta Iteración
Para x1
b −a ∗x k −a ∗x k
(k+1) 1 12 2 13 3
x1 =
a11

2 −( 2 ) ( 0.6025 )−(−3 )
x= (−2.449625)
1
2

x1=−3.2769375
Para x2
b −a ∗x k −a ∗x k
(k+1) 2 21 1 23 3
x2 =
a22

6 −( 0 ) (−3.2769375 ) −(−4 ) (−2.449625 )

x 2= −4
x2=0.949625
Para x3
b −a ∗x k−a ∗x k
(k+1) 3 31 1 32 2
x3 =
a33

23−(1) (−3.2769375)−(3)(0.949625)
x 3= −10
x3=−2.34280625

Reemplazamos

2=2 x1+2 x2−3 x3

2=2 (−3.2769375)+2 (0.949625)−3 (−2.34280625)
2=2.37379375
6=−4 x2−4 x3
6=0(−3.2769375)−4 (0.949625)−4 (−2.34280625 )
6=5.572725
23=x1 +3 x2−10 x3
23=(−3.2769375)+ 3 (0.949625)−10 (−2.34280625)
23=23
lOMoARcPSD|17818758