1

a)
(x0 , y0 ) = (1.31, 3.24) y (x1 , y1 ) = (1.93, 4.76) y aritmética de redondeo a tres
dı́gitos.
Primera Formula
x0 y1 − x1 y0
x1 =
y1 − y0
1.31 × 4.76 − 1.93 × 3.24
x1 =
4.76 − 3.24
6.24 − 6.25
x1 =
1.52
x1 = − 0.00658∗
Segunda Formula

(x1 − x0 )y0
x2 =x0 −
y1 − y0
(1.93 − 1.31) × 3.24
x2 =1.31 −
4.76 − 3.24
0.62 × 3.24
x2 =1.31 −
1.52
2.0088∗
x2 =1.31 −
1.52
x2 =1.31 − 1.32
x2 = − 0.01

Solución Real

Primera formula
(x1 − x0 )y0
x =x0 −
y1 − y0
(1.93 − 1.31) × 3.24
x =1.31 −
4.76 − 3.24
x = − 0.01157894737

Segunda formula
x0 y1 − x1 y0
x=
y1 − y0
1.31 × 4.76 − 1.93 × 3.24
x=
4.76 − 3.24
x = − 0.01157894737
Errores
|x − x1 | | − 0.01157894737 − (−0.00658))| 0.0049989
Er1 = = = = 0.4317231818
|x| | − 0.01157894737| 0.01157894737

1
 |x − x2 | | − 0.01157894737 − (−0.01))| 0.0015789
Er2 = = = = 0.1363595454
|x| | − 0.01157894737| 0.01157894737

Por lo tanto, podemos decir que la segunda formula presenta un menor error
con respecto a la solución real.

b)
Eliminación Gaussiana con sustitución por detrás.

2x1 + 4x2 − x3 = −5
x1 + x2 − 3x3 = −9 (1)
4x1 + x2 + 2x3 = 9

Sea:  
2 4 −1 −5
Ã(1) = 1 1 −3 −9
4 1 2 9
y
a2 1 1
m21 = = = 0.5
a1 1 2
a3 1 4
m31 = = =2
a1 1 2
Hacemos (E2 − m21 E1 ) → E2 y (E3 − m31 E1 ) → E3 , tal que:
 
2 4 −1 −5
Ã(2) = 0 1 − (0.5)(4) −3 − (0.5)(−1) −9 − (0.5)(−5)
0 1 − (2)(4) 2 − (2)(−1) 9 − (0.5)(−5)
 
2 4 −1 −5
Ã(2) = 0 −1 −2.5 −6.5
0 −7 4 11.5

Sea:
a32 −7
m32 = = =7
a11 −1
Hacemos (E3 − m32 E2 ) → E3 , tal que:
 
2 4 −1 −5
Ã(3) = 0 −1 −2.5 −6.5 
0 0 4 − (7)(−2.5) 11.5 − (7)(−6.5)
 
2 4 −1 −5
Ã(3) = 0 −1 −2.5 −6.5
0 0 21.5 57

2
Luego el sistema resultante es:

2x1 +4x2 −1x3 = −5

−x2 −2.5x3 = −6.5
21.5x3 = 57

Por sustitución hacia atrás, tenemos que:

57
x3 = = 2.6511
21.5
x2 =6.5 − 2.5(x3 ) = 6.5 − 2.5(2.6511) = −0.12775
x1 = − 5 − 4(x2 ) + (x3 ) = −5 − 4(−0.12775) + (2.6511) = −1.8379

Por lo tanto, la solución del sistema (1) usando eliminación Gaussiana es:
   
x1 2.6511
x2  = −0.12775
x3 −1.8379

2
Sea: f (x) = x − 2−x = 0, calcular la raı́z que se encuentra en el intervalo [0, 1]
por el método del punto fijo tomando p0 = 0.5 y g(x) = 2−x , para detenerse
use el método de paro: |pn − pn−1 | con una tolerancia de 10−4 .
Por el método del punto fijo tenemos que:

pn = g(pn−1 )

p1 =g(p0 ) = g(0.5) = 2−0.5 = 0.70710

p2 =g(p1 ) = g(0.70710) = 2−0.70710 = 0.61255
p3 =g(p2 ) = g(0.61255) = 2−0.61255 = 0.65403
p4 =g(p3 ) = g(0.65403) = 2−0.65403 = 0.63550
p5 =g(p4 ) = g(0.63550) = 2−0.63550 = 0.64371
p6 =g(p5 ) = g(0.6437) = 2−0.64371 = 0.64006
p7 =g(p6 ) = g(0.64006) = 2−0.64006 = 0.64168
p8 =g(p7 ) = g(0.64168) = 2−0.64168 = 0.64096
p9 =g(p8 ) = g(0.64096) = 2−0.64096 = 0.64128
p10 =g(p9 ) = g(0.64128) = 2−0.64128 = 0.64114
p11 =g(p1 0) = g(0.64114) = 2−0.64114 = 0.64120

3
 Table 1: Resultados de las iteraciones
Iteración i pi−1 pi |pi − pi−1 |
1 0.5 0.70710 0.20710
2 0.70710 0.61255 0.09455
3 0.61255 0.65403 0.04148
4 0.65403 0.63550 0.01853
5 0.63550 0.64371 0.00821
6 0.64371 0.64006 0.00365
7 0.64006 0.64168 0.00162
8 0.64168 0.64096 0.00072
9 0.64096 0.64128 0.00032
10 0.64128 0.64114 0.00014
11 0.64114 0.64120 0.00006 < 10−4

Por tanto, la raı́z aproximada con una tolerancia de 10−4 es 0.64120.

Ahora para demostrar que la función g satisface las condiciones del teorema
que garantizan que la sucesión converge al único punto fijo.
Si p0 es cualquier número en [0,1], entonces la sucesión definida por pn =
g(pn−1 ), n ≥ 1 converge al único punto fijo p en [0,1], es decir pn → p, ∀n > 0.
Si g ′ ∈ C[0, 1]:

lim g ′ (ξn ) = g ′ (p), donde pn < ξn < p

n→∞

Calculamos g ′ :
ln(2)
g(x) = 2−x −→ g ′ (x) = − = −ln(2) · 2x , x ∈ [0, 1]
2−x
Luego, nuestro punto fijo p es 0.64120, entonces:
ln(2)
g ′ (x) = −
2−x
ln(2)
g ′ (0.64120) = −
2−0.64120
g ′ (0.64120) = −1.08105
Como g ′ (p) ̸= 0 entonces la iteración del punto fijo es lineal. Por lo tanto, se
demuestra que las iteraciones convergen linealmente al punto fijo p = 0.64120.

3
Una matriz es estrictamente dominante diagonalmente por filas si y solo si:
n
X
|aii | > |aij |, para cada i = 1, 2, . . . , n.
j=1,j̸=i

4
Donde n es la dimensión de la matriz.
Del sistema:
2x1 − x2 +10x3 = −11
3x2 −x3 +8x4 = −11
10x1 − x2 +2x3 =6
−x1 + 11x2 −x3 +3x4 = 25
Sea la matriz:  
2 −1 10 0
0 3 −1 8
A= 
 10 −1 2 0
−1 11 −1 3
Entonces para que (A) sea una matriz estrictamente dominante diagonalmente
buscaremos los elementos de la diagonal que cumplan con la definición.
Para encontrar el primer elemento hacemos:

|a11 | =|a12 | + |a13 | + |a14 | →|2| > | − 1| + |10| + |0| = 11, 2 ≯ 11

|a21 | =|a22 | + |a23 | + |a24 | →|0| > |3| + | − 1| + |8| = 12, 0 ≯ 12
|a31 | =|a32 | + |a33 | + |a34 | →|10| > | − 1| + |2| + |0| = 3, 10 > 3
|a41 | =|a42 | + |a43 | + |a44 | →| − 1| > |11| + | − 1| + |3| = 12, 1 ≯ 12

Ahora con el primer elemento encontrado reordenamos (F1 ←→ F3 ):

 
10 −1 2 0
0 3 −1 8
A=  2 −1 10 0


−1 11 −1 3

Para encontrar el segundo elemento hacemos:

|a22 | =|a21 | + |a23 | + |a24 | →|3| > |0| + | − 1| + |8| = 9, 3 ≯ 9

|a32 | =|a31 | + |a33 | + |a34 | →| − 1| > |2| + |10| + |0| = 12, 1 ≯ 12
|a42 | =|a41 | + |a43 | + |a44 | →|11| > | − 1| + | − 1| + |3| = 5, 11 > 5

Ahora con el segundo elemento encontrado reordenamos (F2 ←→ F4 ):

 
10 −1 2 0
−1 11 −1 3
A=  2 −1

10 0
0 3 −1 8

Para encontrar el tercer y cuarto elemento verificamos que:

|a33 | =|a31 | + |a32 | + |a34 | →|10| > |2| + | − 1| + |0| = 3, 10 > 3

|a44 | =|a41 | + |a42 | + |a43 | →|8| > |0| + |3| + | − 1| = 4, 8 > 4

5
Ahora con los elementos encontrados, nuestra matriz estrictamente dominante
diagonalmente será:  
10 −1 2 0
−1 11 −1 3
A= 
2 −1 10 0
0 3 −1 8
Entonces, el sistema resultante es:

10x1 − x2 +2x3 =6
−x1 + 11x2 −x3 +3x4 = 25
2x1 − x2 +10x3 = −11
3x2 −x3 +8x4 = −11

Por el método de Jacobi, tenemos:

(k−1) (k−1)
(k) 6 + x2 − 2x3 (k−1) (k−1)
x1 = = 0.6 + 0.1 · x2 − 0.2 · x3
10
(k−1) (k−1) (k−1)
(k) 25 + x1 + x3 − 3x4 (k−1) (k−1) (k−1)
x2 = = 2.273 + 0.091 · x1 + 0.091 · x3 − 0.27 · x4
11
(k−1) (k−1)
(k) −11 − 2x1 + x2 (k−1) (k−1)
x3 = = −1.1 − 0.2 · x1 + 0.1 · x2
10
(k−1) (k−1)
(k) −11 − 3x2 + x3 (k−1) (k−1)
x4 = = −1.375 − 0.375 · x2 + 0.125 · x3
8
(2)
Entonces, matricialmente tenemos:
 (k)     (k−1)  
x1 x1

0 0.1 0.2 0 0.6
 (k)    (k−1)  
x2  0.091 0 0.091 −0.27 x 2.273 
 (k)  =  · 2(k−1)  + 
 
x3  −0.2 0.1 0 0  x3   −1.1 
x
(k) 0 −0.375 0.125 0 x
(k−1) −1.375
4 4

Donde:  
0 0.1 0.2 0
0.091 0 0.091 −0.27
Tj = 
 
−0.2 0.1 0 0 
0 −0.375 0.125 0
Entonces:
||Tj || = max{|0| + |0.1| + |0.2| + |0|,
|0.091| + |0| + |0.091| + | − 0.27|,
| − 0.2| + |0.1| + |0| + |0|,
|0| + | − 0.375| + |0.125| + |0|}
||Tj || = max{0.3, 0.452, 0.3, 0.5}
||Tj || = 0.5

6
Luego, calculo de la primera iteración, con x(0) = (0, 0, 0, 0)t . Entonces para
k = 1 en (2),tenemos:
(1) (0) (0)
x1 =0.6 + 0.1 · x2 − 0.2 · x3 = 0.6
(1) (0) (0) (0)
x2 =2.273 + 0.091 · x1 + 0.091 · x3 − 0.27 · x4 = 2.273
(1) (0) (0)
x3 = − 1.1 − 0.2 · x1 + 0.1 · x2 = −1.1
(1) (0) (0)
x4 = − 1.375 − 0.375 · x2 + 0.125 · x3 = −1.375

−→ x(1) = (0.6, 2.273, −1.1, −1.375)t

Luego, calculo de la segunda iteración, para k = 2 en (2),tenemos:
(2) (1) (1)
x1 =0.6 + 0.1 · x2 − 0.2 · x3
(2)
x1 =0.6 + 0.1 · 2.273 − 0.2(−1.1) = 1.0473
(2) (1) (1) (1)
x2 =2.273 + 0.091 · x1 + 0.091 · x3 − 0.27 · x4
(2)
x2 =2.273 + 0.091 · 0.6 + 0.091(−1.1) − 0.27(−1.375) = 2.59875
(2) (1) (1)
x3 = − 1.1 − 0.2 · x1 + 0.1 · x2
(2)
x3 = − 1.1 − 0.2 · 0.6 + 0.1 · 2.273 = −0.9927
(2) (1) (1)
x4 = − 1.375 − 0.375 · x2 + 0.125 · x3
(2)
x4 = − 1.375 − 0.375 · 2.273 + 0.125(−1.1) = −2.364875

−→ x(2) = (1.0473, 2.59875, −0.9927, −2.364875)t

Ahora, hallamos la cota del error, tal que:

||Tj ||k
|x − x(0) | ≤ · ||x(1) − x(0) ||
1 − ||Tj ||

Donde:
||x(1) − x(0) || = max{|0.6|, |2.273|, | − 1.1|, | − 1.375|}
||x(1) − x(0) || = 2.273
Entonces:
||Tj ||k (0.5)k
|x − x(0) | ≤ · ||x(1) − x(0) || = · 2.273
1 − ||Tj || 1 − 0.5

||Tj ||k (0.5)k

|x − x(0) | ≤ · ||x(1) − x(0) || = · 2.273
1 − ||Tj || 0.5
||Tj ||k 2.273
|x − x(0) | ≤ · ||x(1) − x(0) || = · (0.5)k
1 − ||Tj || 0.5

7
 ||Tj ||k
|x − x(0) | ≤ · ||x(1) − x(0) || = 4.546 · (0.5)k
1 − ||Tj ||
Por tanto, para k = 1, la cota del error es:

||Tj ||k
|x − x(0) | ≤ · ||x(1) − x(0) || = 4.546 · (0.5)1
1 − ||Tj ||

||Tj ||k
|x − x(0) | ≤ · ||x(1) − x(0) || = 2.273
1 − ||Tj ||
Ası́ mismo, para k = 2, la cota del error es:

||Tj ||k
|x − x(0) | ≤ · ||x(1) − x(0) || = 4.546 · (0.5)2
1 − ||Tj ||

||Tj ||k
|x − x(0) | ≤ · ||x(1) − x(0) || = 1.1365
1 − ||Tj ||