Universidad del Magdalena

Taller 2 de Análisis Numérico

Solución de ecuaciones no lineales

Docente. Álvaro Espinoza Pérez

Guerra Zequeda Álvaro

Mercado Fernández Carlos

21/11/2021
Colombia
SOLUCIÓN DE TALLER.
PUNTO 1.
f(x) = (x − 2)2− ln(x)
MÉTODO DE BISECCION
Intervalo [1,2]
Exactitud 10-4
Se calcula el número de interacciones que sean necesarias para la aproximación
𝑏−𝑎
𝑙𝑜𝑔( )
n≥ 𝑡
log⁡(2)

b=2
a=1
t= 10-4
2−1
𝑙𝑜𝑔( −4 )
n≥ 10
log⁡(2)

log⁡(10000)
n⁡≥ log(2)

n ≥ 13,28771238 ≈ 14
n= 14
El número de interacciones para llegar a la aproximación es de n=14
a1= 1
f (1) = 1
b1= 2
f (2) = (2-2)2 – Ln (2)
f (2) = -0,6931471806
Se cumple que
⁡⁡f(a1). f(b1) < 0
calcula P1 en [1,2]
𝑎1 +⁡𝑏1
P1= 2
1+2
P1= 2
P1= 1,5
f (1,5) = (1,5 – 2)2 – Ln (1,5)
f (1,5) = -0,1554651081 < 0
Como f(P1) tiene el mismo signo que f(b1) entonces a2= a1 y b2=P1
a2= 1; b2= 1,5

Segunda iteración
Se calcula P2 en [1,1.5]
𝑎2 +⁡𝑏2
P2= 2
1+1,5
P2= 2

P2= 1,25
f (1,25) = (1,25 – 2)2 – Ln (1,25)
f (1,25) = 0,3393564487
Como f(P2) tiene el mismo signo que f(a2) entonces a3= P2 y b3=b2
a3= 1,25; b3= 1,5
Continuamos con las iteraciones en la tabla.
n an(+) bn(-) pn F(pn) Ea(%)
1 1 2 1,5 -0,1554651081
2 1 1,5 1,25 0,3393564487 20%
3 1,25 1,5 1,375 0,0721712689 9,0909%
4 1,375 1,5 1,4375 -0,0464992437 4,3478%
5 1,375 1,4375 1,40625 0,0116124755 2,2222%
6 1,40625 1,4375 1,421875 -0,0177479075 1,0989%
7 1,40625 1,421875 1,4140625 -0,0031440134 0,5525%
8 1,40625 1,4140625 1,41015625 0,0042151356 0,2770%
9 1,41015625 1,4140625 1,412109375 0,0005307898 0,1383%
10 1,412109375 1,4140625 1,413085938 -0,0013078043 0,0691%
11 1,412109375 1,413085938 1,412597656 -0,0003888054 0,0346%
12 1,412109375 1,412597656 1,412353516 0,0000709177 0,0173%
13 1,412353516 1,412597656 1,412475586 -0,0001589625 0,0086%
14 1,412353516 1,412475586 1,412414551 -0,0000440271 0,0043%
 MÉTODO DE NEWTON-RAPHSON
Intervalo [1,2]
X0=1,5
f(x)= (x-2)2 – Ln(x)
Derivamos
1
f´(x)= 2(x-2) - 𝑥

Aplicamos la ecuación
𝑓(𝑥 )
Xn+1= xn - 𝑓´(𝑥𝑛
𝑛)

Primera iteración
Calculamos X2
𝑓(𝑋1 )
X2= X1 - 𝑓´(𝑥1)

𝑓(1,406720935)⁡
X2= 1,406720935 -⁡𝑓´(1,406720935)
0,01071863069
X2= 1,406720935 - −1,897431179

X2= 1,412369957S
Tabla de interacción del método de newton-Raphson

n Xn F(Xn) F´(Xn) Ea(%)

1 1,5 -0,1554651081 -1,666666667
2 1,406720935 0,01071863069 -1,897431179 6,6310%
3 1,412369957 0,00003995300 -1,883289872 0,4000%
4 1,4123911717 0,00000000056 -1,883236809 0,0015%
5 1,4123911720 0,00000000000 -1,883236808 0,0000%

MÉTODO DE LA POSICIÓN FALSA

Intervalo [1,2]
F(x)= (x-2)2 – Ln (x)
X0 =1; X1=2
 f (0) = f (1) = 1
f(X1) = f (2) = - 0, 6931471806
se cumple que f(x0). f(x1) < 0
Tabla de interacción; método de la regule falsa
n Xn f(Xn) Xn+1 f(Xn+1) Xr f(Xr) Ea(%)

0 1 1 2 -0,6931471806 1,590616109 -
0,2965262611
1 1 1 1,590616109 -0,2965262611 1,455537328 - 9,2803%
0,0789355295
2 1 1 1,455537328 -0,0789355295 1,422209961 - 2,3434%
0,0183706441
3 1 1 1,422209961 -0,0183706441 1,41459361 - 0,5384%
0,0041416467
4 1 1 1,41459361 -0,0041416467 1,412883592 - 0,1210%
0,0009270400
5 1 1 1,412883592 -0,0009270400 1,412501187 - 0,0271%
0,0002071687
6 1 1 1,412501187 -0,0002071687 1,412415747 - 0,0060%
0,0000462800
7 1 1 1,412415747 -0,0000462800 1,412396661 - 0,0014%
0,0000103378
8 1 1 1,412396661 -0,0000103378 1,412392398 - 0,0003%
0,0000023092
9 1 1 1,412392398 -0,0000023092 1,412391446 - 0,0001%
0,0000005158
10 1 1 1,412391446 -0,0000005158 1,412391233 - 0,0000%
0,0000001152

segunda iteración
𝑐𝑎𝑙𝑐𝑢𝑙𝑎𝑚𝑜𝑠⁡⁡𝑥𝑟 ⁡⁡; de la siguiente forma:
𝒇(𝒙𝟎 )(𝒙𝟏 − 𝒙𝟎 )
𝒙𝒓 = ⁡ 𝒙𝟎 − ⁡
𝒇(𝒙𝟏 ) − 𝒇(𝒙𝟎 )
1(2 − 1)
𝑥𝑟 = 1 −
−0,6931471806 − 1
𝑥𝑟 = 1,590616109
Verificamos;
𝑓(𝑥𝑟 ) = 𝑓(1,590616109) = −0,2965262611
Como 𝑓(𝑥𝑟 ) tiene el mismo signo que 𝑓(𝑥1 )entonces; 𝑥1 = 𝑥0 ⁡𝑦⁡𝑥2 = 𝑥𝑟 ⁡⁡ quedaría de la
siguiente manera = 𝑥1 = 1 y 𝑥2 = 1,590616109
𝐴ℎ𝑜𝑟𝑎⁡𝑠𝑒⁡𝑐𝑎𝑙𝑐𝑢𝑙𝑎⁡⁡𝑥𝑟 ⁡⁡; de la siguiente forma:
𝒇(𝒙𝟏 )(𝒙𝟐 − 𝒙𝟏 )
𝒙𝒓 = ⁡ 𝒙𝒓 − ⁡
𝒇(𝒙𝟐 ) − 𝒇(𝒙𝟏 )
1(1,590616109 − 1)
𝑥𝑟 = 1 −
−0,2965262611 − 1
𝑥𝑟 = 1,455537328
Verificamos:
𝑓(𝑥𝑟 ) = 𝑓(1,455537328) = −0,0789356293
Método de la secante.
𝑖𝑛𝑡𝑒𝑟𝑣𝑎𝑙𝑜⁡[1,2]⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡𝑓(𝑥) = (𝑥 − 2)2 − 𝐿𝑛 (𝑥)
𝑥0 = 1⁡⁡⁡⁡⁡⁡⁡⁡;⁡⁡⁡⁡⁡⁡𝑥1 = 2
𝑓(𝑥0 ) = 𝑓(1) = 1
𝑓(𝑥1 ) = 𝑓(2) = −0,6931471806
Primera iteración.
Calculamos 𝑥2 ; de la siguiente forma:
𝒇(𝒙𝟏 )(𝒙𝟏 − 𝒙𝟎 )
𝒙𝟐 = ⁡ 𝒙𝟏 − ⁡
𝒇(𝒙𝟏 ) − 𝒇(𝒙𝟎 )
(−0,6931471806))(2 − 1)
𝑥2 = ⁡2 − ⁡
−0,6931471806 − 1
𝑥2 = 1,590616109
Verificamos:
𝑓(𝑥2 ) = 𝑓(1,590616109) = −0,2965262611
Segunda iteración
Calculamos 𝑥3⁡ ;⁡de la siguiente forma:
𝒇(𝒙𝟐 )(𝒙𝟐 − 𝒙𝟏 )
𝒙𝟑 = ⁡ 𝒙𝟐 −
𝒇(𝒙𝟐 ) − 𝒇(𝒙𝟏 )
(−0,2965262611)(1,590616109 − 2)
𝑥3 = 1,590616109 −
(−0,2965262611) + 0,6931471806
𝑥3 = 1,2845478496
Verificamos:
 𝑓(𝑥3 ) = 𝑓(1,2845478496) = 0,2614649911
Continuamos con las interacciones en la siguiente tabla del método de la secante:
i Xi-1 F(Xi-1) Xi F(Xi) Xi+1 F(Xi+1) Ea(%)
1 1 1 2 -0,6931471806 1,590616109 -0,2965262611
2 2 -0,6931471806 1,590616109 -0,2965262611 1,2845478496 0,2614649911 23,8269%
3 1,590616109 -0,2965262611 1,28454785 0,2614649911 1,4279661101 -0,0290283600 10,0435%
4 1,28454785 0,2614649911 1,42796611 -0,0290283600 1,4136346396 -0,0023398103 1,0138%
5 1,42796611 -0,0290283600 1,41363464 -0,0023398103 1,4123781860 0,0000244559 0,0890%
6 1,41363464 -0,0023398103 1,412378186 0,0000244559 1,4123911828 -0,0000000202 0,0009%
7 1,412378186 0,0000244559 1,412391183 -0,0000000202 1,4123911720 0,0000000000 0,0000%

Verificación en OCTAVE:
Gráfica de la función:

PUNTO 2.
𝟐 −𝟏
𝐞√𝐱 = 𝟐√𝐱 + 𝟏, 𝟐𝟒

𝒙𝟎 = 𝟏, 𝟓
𝟐 −𝟏
Punto fijo 𝒇(𝒙) = 𝐞√𝐱 − 𝟐√𝐱 − 𝟏, 𝟐𝟒

Despejamos x
𝟐
𝒈(𝒙) = √[𝑳𝒏⁡ (𝟐√𝒙 + 𝟏, 𝟐𝟒)] + 𝟏

𝟐
𝟐 −𝟏
𝒆√𝒙 − 𝟏, 𝟐𝟒
𝒈(𝒙) = ( )
𝟐

𝟐
Se toma 𝒈(𝒙) = √[𝑳𝒏⁡ (𝟐√𝒙 + 𝟏, 𝟐𝟒)] + 𝟏

Si 𝑥0 = 1,5⁡⁡⁡⁡; ⁡⁡𝒆𝒏𝒕𝒐𝒏𝒄𝒆𝒔⁡ → 𝑥𝑛+1 = 𝑔(𝑥𝑛 )

Primera iteración.
Hallamos x1⁡ de⁡la⁡siguiente⁡forma:⁡
𝟐
𝑥1 → 𝑔(𝑥0 ) = √[𝑳𝒏⁡ (𝟐√𝒙𝟎 + 𝟏, 𝟐𝟒)] + 𝟏

𝟐
𝑥1 = √[𝑳𝒏⁡ (𝟐√𝟏, 𝟓 + 𝟏, 𝟐𝟒)] + 𝟏

𝑥1 = 1,644475404
Segunda iteración.
Hallamos x2⁡ de⁡la⁡siguiente⁡forma:⁡
 𝟐
𝑥2 → 𝑔(𝑥1 ) = √[𝑳𝒏⁡ (𝟐√𝒙𝟏 + 𝟏, 𝟐𝟒)] + 𝟏

𝟐
𝑥2 = √[𝑳𝒏⁡ (𝟐√𝟏, 𝟔𝟒𝟒𝟒𝟕𝟓𝟒𝟎𝟒 + 𝟏, 𝟐𝟒)] + 𝟏

𝑥2 = 1,668999467
Continuamos con las iteraciones en la tabla.
Tabla de iteración del método de punto fijo

n Xn F(Xn) Ea(%)
0 1,5 -0,6306551581
1 1,644475404 -0,1152522749 8,786%
2 1,668999467 -0,0190532272 1,469%
3 1,673001491 -0,0030959327 0,239%
4 1,673650393 -0,0005016363 0,039%
5 1,673755499 -0,0000812433 0,006%
6 1,673772521 -0,0000131569 0,001%
7 1,673775277 -0,0000021307 0,000%

OCTAVE
Graficamos en OCTAVE:

PUNTO 3.
2
𝑓(𝑥) = 𝐿𝑛 (𝑥) + −1
𝑥2 + 1
Donde 𝑥0 = 40
1 4𝑥
𝑠𝑒⁡𝑑𝑒𝑟𝑖𝑣𝑎⁡𝑓´(𝑥) = − 2
𝑥 (𝑥 + 1)2
𝑓(𝑥𝑛 )
𝑠𝑒⁡𝑎𝑝𝑙𝑖𝑐𝑎⁡𝑙𝑎⁡𝑒𝑐𝑢𝑎𝑐𝑖𝑜𝑛⁡⁡𝑥𝑛 + 1 = 𝑥𝑛 −
𝑓´(𝑥𝑛 )
Primera iteración
𝑓(𝑥0 )
𝒙𝟏 = 𝒙𝟎 −
𝑓´(𝑥0 )
𝑓(40)
𝑥1 = 40 − ⁡⁡⁡
𝑓´(40)
2,689504064
⁡x1 = 40 − ⁡⁡⁡⁡⁡
0,02493751805
x1 = −67,87449638
Segunda iteración
𝑓(𝑥1)
𝒙𝟐 = 𝒙𝟏 −
𝑓´(𝑥1 )
f(−67,87449633)
x2 = −67,87449638 − ⁡⁡⁡⁡⁡
f´(−67,87449633)
x2 ⁡𝑛𝑜⁡𝑒𝑠⁡𝑝𝑜𝑠𝑖𝑏𝑙𝑒⁡𝑐𝑎𝑙𝑐𝑢𝑙𝑎𝑟⁡𝑒𝑙⁡𝑣𝑎𝑙𝑜𝑟⁡𝑑𝑎𝑑𝑜⁡𝑞𝑢𝑒⁡⁡𝑝𝑜𝑟⁡𝑠𝑒𝑟⁡𝑢𝑛⁡𝑛𝑢𝑚𝑒𝑟𝑜⁡𝑛𝑒𝑔𝑎𝑡𝑖𝑣𝑜, 𝑒𝑙⁡𝑳𝒏 (𝒙)⁡
𝑛𝑜⁡𝑠𝑒⁡𝑝𝑢𝑒𝑑𝑒⁡ℎ𝑎𝑙𝑙𝑎𝑟
Para x0 = 3
Primera iteración
𝑓(𝑥0 )
𝒙𝟏 = 𝒙𝟎 −
𝑓´(𝑥0 )
𝑓(3)
𝑥1 = 3 −
𝑓´(3)
0,29861228867
𝒙𝟏 =
0,21333333333
𝑥1 = 1,600254897
Segunda iteración
𝑓(𝑥1 )
𝒙𝟐 = 𝒙𝟏 −
𝑓´(𝑥1 )
𝑓(1,600254897)
𝒙𝟐 = 1,600254897 −
𝑓´(1,600254897)
0,0319319795
𝒙𝟐 = 1,600254897 −
0,1200646346
𝑥2 = 1,3351312016
Tabla de interacción; método de newton-Raphson

n Xn F(Xn) F´(Xn) Ea(%)

0 3 0,2986122887 0,2133333333
1 1,600254897 0,0318319796 0,1200646346 87,4701%
2 1,335131202 0,0077881850 0,0592426527 19,8575%
3 1,203668739 0,0020945874 0,0279075182 10,9218%
4 1,128614146 0,0005869441 0,0128449628 6,6502%
5 1,082919652 0,0001680782 0,0058352333 4,2196%
6 1,054115629 0,0000487395 0,0026300545 2,7325%
7 1,035583866 0,0000142423 0,0011796077 1,7895%
8 1,023510067 0,0000041820 0,0005274091 1,1796%
9 1,01558083 0,0000012318 0,0002353284 0,7808%
10 1,010346634 0,0000003635 0,0001048631 0,5181%
11 1,006879877 0,0000001074 0,0000466864 0,3443%
12 1,004578685 0,0000000318 0,0000207734 0,2291%
13 1,00304896 0,0000000094 0,0000092397 0,1525%
14 1,002031089 0,0000000028 0,0000041086 0,1016%
15 1,001353372 0,0000000008 0,0000018267 0,0677%
16 1,000901942 0,0000000002 0,0000008120 0,0451%
17 1,000601159 0,0000000001 0,0000003610 0,0301%

OCTAVE

Graficamos…
PUNTO 4
A. En particular cuando 𝑃(𝑥) se hace muy pequeña en el intervalo comprendido entre un
punto de partida y su raíz, el punto dado por la iteración siguiente tiene un valor muy
grande, puesto que la tangente es casi paralela y en este caso el método diverge.

B. Método de Horner
𝑃(𝑥) = 𝑥 3 − 3,9𝑥 2 + 4,79𝑥 − 1,881
⁡𝑥0 = 1
1⁡⁡⁡|1 − 3,9⁡⁡⁡4,79⁡1,881
⁡⁡⁡⁡|⁡⁡⁡⁡⁡⁡⁡⁡1⁡⁡⁡⁡⁡⁡⁡1⁡⁡ − 1,9
⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡1⁡ − 2,9⁡⁡1,89 − 0,004
⁡⁡𝑝(1) = 1⁡ − 1,9⁡ − 0,01
Donde ⁡⁡𝑄(𝑥) = 𝑥 2 − 2,9𝑥 + 1,89
⁡⁡𝑝′ (1) = ⁡𝑄(1)
𝑝′ (1) puede hallarse al averiguar 𝑄(1), a si
⁡𝑥0 = 1⁡⁡⁡⁡|1 − 2,9⁡⁡⁡1,89
⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡|⁡⁡⁡⁡⁡⁡⁡⁡⁡1⁡⁡ − 1,9
⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡1⁡ − 1,9⁡ − 0,01
1⁡ − 1,9⁡ − 0,01 = 𝑄(1) = 𝑃′(1)
Y al actualizar Newton-Raphson tenemos
𝑝(𝑥0 )
𝑥1 = 𝑥0 −
𝑝´(𝑥0 )
0,009
𝑥1 = 1 −
−0,01
𝑥1 = 1,9
Segunda iteración
1,9|1 − 3,9⁡⁡⁡4,79⁡⁡⁡⁡1,881
⁡⁡⁡⁡⁡⁡⁡|⁡⁡⁡⁡⁡⁡⁡1,9⁡ − 3,8⁡⁡1,9881
⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡1⁡⁡ − 2⁡⁡⁡⁡⁡0,99⁡⁡⁡⁡⁡0⁡⁡⁡
⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡1,9 − 0,19⁡⁡⁡⁡
1 -0,1 0,8= 𝑄(1,9) = 𝑃′(1,9)
𝑠𝑖⁡⁡𝑃(⁡𝑥1 ) = 0⁡⁡⁡⁡𝑦⁡⁡⁡⁡⁡p′ (𝑥1 ) = 0,8
 𝑝(𝑥1 )
𝑥2 = 𝑥1 −
𝑝´(𝑥1 )
0
𝑥2 = 1,9 −
0,8
𝑥2 = 1,9
El problema no se puede solucionar dado que los dos métodos calculan la misma raíz
r = 1,9⁡por⁡lo⁡tanto⁡las⁡dos⁡nos⁡sirven⁡para⁡hallar⁡la⁡primera⁡y⁡segunda⁡raiz
C. 𝑃(𝑥) = 𝑥 3 − 3,9𝑥 2 + 4,79𝑥 − 1,881
𝑥0 = 1,2
Método Newton-Raphson
Derivamos 𝑃(𝑥)
𝑃(𝑥) = 3𝑥 2 − 7,8𝑥 + 4,79
Aplicamos la ecuación
𝑓(𝑥𝑛 )
𝒙𝒏+𝟏 = 𝒙𝒏 −
𝑓´(𝑥𝑛 )
Primera iteración
𝑓(𝑥0 )
𝑥1 = 𝑥0 −
𝑓´(𝑥0 )
𝑓(1,2)
𝑥1 = 1,2 −
𝑓´(1,2)
(−0,021)
𝑥1 = 1,2 −
(−0,25)
𝑥1 = 1,116
Segunda iteración
Calculamos 𝑥2
𝑓(𝑥1 )
𝑥2 = 𝑥1 −
𝑓´(𝑥1 )
𝑓(1,116)
𝑥2 = 1,116 −
𝑓´(1,116)
(−0.0027095040)
𝑥2 = 1,116 −
(−0,1784320000)
𝑥2 = 1,1008149211
Tabla de interacción del método de newton-Raphson

n Xn F(Xn) F´(Xn) Ea(%)

1 1,2 -0,021 -0,25
2 1,116 -0,0027095040 -0,178432 7,5269%
3 1,100814921 -0,0001307853 -0,160975913 1,3794%
4 1,100002469 -0,0000003950 -0,160002962 0,0739%
5 1,1 0,0000000000 -0,16 0,0002%
6 1,1 0,0000000000 -0,16 0,0000%

OCTAVE
PUNTO 5
𝑓(𝑥) = cos(𝑥) + 𝑒 2 − 𝑥 2 + 1
𝑥0 = −2
Método de Newton-Raphson modificado
Derivamos 𝑓(𝑥)
𝑓′(𝑥) = −𝑠𝑒𝑛𝑥 + 𝑒 𝑥 − 2𝑥
𝑓 ′′ (𝑥) = −𝑐𝑜𝑠𝑥 + 𝑒 𝑥 − 2
Aplicamos la ecuación
𝑓(𝑥𝑛 )𝑓′(𝑥𝑛 )
𝑥𝑛+1 = 𝑥𝑛 −
[𝑓′(𝑥𝑛 )]2 − 𝑓(𝑥𝑛 )𝑓′′(𝑥𝑛 )
Primera iteración
𝑓(𝑥0 )𝑓′(𝑥0 )
𝑥1 = 𝑥0 −
[𝑓′(𝑥0 )]2 − 𝑓(𝑥0 )𝑓′′(𝑥0 )
(−3,2808115533)(5,0446327101)
𝑥1 = −2 −
[5,0446327101]2 − (−3,2808115533)(−1,4485178802)
𝑥1 = −1,2003051173
Segunda iteración
𝑓(𝑥1 )𝑓′(𝑥1 )
𝑥2 = 𝑥1 −
[𝑓′(𝑥1 )]2 − 𝑓(𝑥1 )𝑓′′(𝑥1 )
𝑓(−1,2003051173)𝑓′(−1,2003051173)
𝑥2 = −1,2003051173 −
[𝑓′(−1,2003051173)]2 − 𝑓(−1,2003051173)𝑓′′(= −1,2003051173)
(0,2224433082)(3,6338621651)
𝑥2 = −1,2003051173 −
[3,6338621651]2 − (0,3224433062)(−2,0609710300)
𝑥2 = −1,2594652178
Continuamos con las iteraciones en la tabla.
Tabla de interacción del método de newton-Raphson Modificado

n Xn F(Xn) F´(Xn) F´´(Xn) Ea(%)

0 -2 -3,2808115533 5,0446327101 -1,4485178802
1 -1,200305117 0,2224433082 3,6338621651 -2,0609710300 66,6243%
2 -1,259465218 0,0038791509 3,7546628558 -2,0225202654 4,6972%
3 -1,260497799 0,0000010796 3,7567509154 -2,0218300621 0,0819%
 4 -1,260498086 0,0000000000 3,7567514964 -2,0218298699 0,0000%

OCTAVE

Graficamos…

PUNTO 6
𝑓(𝑥) = 𝑥 3 − 3𝑥 + 2;
Con 𝑥0 = −2,4
Método de Halley
Derivamos 𝑓(𝑥)
𝑓 ′(𝑥) = 3𝑥 3 − 3

𝑓 ′(𝑥) = 6𝑥
Aplicamos la ecuación
𝑓(𝑥𝑛 ) 𝑓(𝑥𝑛 )𝑓′′(𝑥1 ) −1
𝑥𝑛+1 = 𝑥𝑥 − (1 − ⁡⁡⁡⁡)
𝑓′(𝑥𝑛 ) 2(𝑓′′(𝑥𝑥 ))2
Primera iteración
𝑓(𝑥0 ) 𝑓(𝑥0 )𝑓′′(𝑥0 ) −1
𝑥1 = 𝑥0 − (1 − ⁡⁡⁡⁡)
𝑓′(𝑥0 ) 2(𝑓′(𝑥0 ))2
𝑓(−2,4) 𝑓(−2,4 )𝑓′′(−2,4) −1
𝑥1 = −2,4 − (1 − ⁡⁡⁡⁡)
𝑓(−2,4) 2(𝑓′(−2,4))2
𝑓(−4,624) (−4,624)(−14,4) −1
𝑥1 = −2,4 − (1 − ⁡⁡⁡⁡)
14,28 2(14,28)2
𝑥1 = −2,01301
Segunda iteración
𝑓(𝑥1 ) 𝑓(𝑥1 )𝑓′′(𝑥1 ) −1
𝑥2 = 𝑥1 − (1 − ⁡⁡⁡⁡)
𝑓′(𝑥1 ) 2(𝑓′(𝑥1 ))2
𝑓(−2,01301) 𝑓(−2,01301)𝑓 ′′ (−2,01301) −1
𝑥2 = −2,01301 − (1 − ⁡⁡⁡⁡)
𝑓 ′ (−2,01301) 2(𝑓′(−2,01301))2
(−0,11809) (−0,11809)(−12,07805) −1
𝑥2 = −2,01301 − (1 − ⁡⁡⁡⁡)
(9,15661) 2(9,15661)2
𝑥2 = −2,00000
Tabla de interacción del método de Halley

n Xn F(Xn) F´(Xn) F´´(Xn) Ea(%)

0 -2,4 -4,6240000000 14,2800000000 -14,4000000000
1 -2,01300813 -0,1180906405 9,1566051953 -12,0780487805 19,2246%
2 -2,000000721 -0,0000064901 9,0000086534 -12,0000043267 0,6504%
3 -2 0,0000000000 9,0000000000 -12,0000000000 0,0000%

Método de Newton-Raphson
𝑓(𝑥) = 𝑥 3 − 3𝑥 + 2;
Con 𝑥0 = −2,4
Derivamos 𝑓(𝑥)

𝑓 ′(𝑥) = 3𝑥 3 − 3
Aplicamos la ecuación
𝑓(𝑥𝑛 )
𝑥𝑛+1 = 𝑥𝑥 −
𝑓′(𝑥𝑛 )
Primeras iteraciones
𝑓(𝑥0 )
𝑥1 = 𝑥0 −
𝑓′(𝑥0 )
𝑓(−2,4)
𝑥1 = −2,4 −
𝑓 ′ (−2,4)
𝑓(−4,62,4)
𝑥1 = −2,4 −
𝑓 ′ (14,28)
𝑥1 = −2,07619
Segundas iteraciones
𝑓(𝑥1 )
𝑥2 = 𝑥1 −
𝑓′(𝑥1 )
𝑓(−2,07619 )
𝑥1 = −2,07619 −
𝑓 ′ (−2,07619)
(−0,72099)
𝑥1 = −2,07619 −
(9,93110)
𝑥1 = −2,00360

Tabla de interacción del método de newton-Raphson

n Xn F(Xn) F´(Xn) Ea(%)
0 -2,4 -4,62400000 14,28000000
1 -2,07619048 -0,72098650 9,93170068 15,5963%
2 -2,00359601 -0,03244173 9,04319092 3,6232%
3 -2,00000859 -0,00007731 9,00010308 0,1794%
4 -2,00000000 0,00000000 9,00000000 0,0004%
5 -2,00000000 0,00000000 9,00000000 0,0000%
En conclusión, se puede evidenciar que el método de Halley utiliza solo (3) iteraciones para
hallar la raíz de la función 𝑓(𝑥) = 𝑥 3 − 3𝑥 + 2; en el puto inicial 𝑥0 = −2,4, en cambio el
método de Newton-Raphson.
OCTAVE

Graficamos…

PUNTO 7
𝑓(𝑥) = 𝑥 2 − cos(𝑥) − 1

𝑔(𝑥) = √cos(𝑥) + 1
Método de punto fijo
𝑥0 = 2 entonces sabemos que
𝑥𝑛+1 = 𝑔(𝑥𝑛)

Primera iteración

𝑥1 = 𝑔(𝑥0 ) = √cos(𝑥0 ) + 1

𝑥1 = √cos(2) + 1
𝑥1 = 0,7641028487
Segundad interacción

𝑥2 = 𝑔(𝑥1 ) = √cos(𝑥1 ) + 1

𝑥2 = √cos(0,7641028481) + 1
𝑥2 = 1,312251263
Tabla de interacción del método de punto fijo

n Xn F(Xn) Ea(%)
0 2 3,4161468365
1 0,764102849 -1,1381502140 161,745%
2 1,312251263 0,4663291330 41,772%
3 1,120568715 -0,1794962310 17,106%
4 1,197986008 0,0709363379 6,462%
5 1,16800434 -0,0277542740 2,567%
6 1,179825585 0,0109023281 1,002%
7 1,175196189 -0,0042759608 0,394%
8 1,177014037 0,0016780832 0,154%
9 1,176300965 -0,0006583991 0,061%
10 1,176580792 0,0002583484 0,024%
11 1,176470999 -0,0001013693 0,009%
12 1,17651408 0,0000397753 0,004%
13 1,176497176 -0,0000156069 0,001%
14 1,176503809 0,0000061238 0,001%
15 1,176501206 -0,0000024029 0,000%

Método de Aiken
𝑥0 = 2 𝑥2 = 1,312251263
𝑥1 = 0,764102849 𝑥3 =1,120568715
Utilizaremos las intercesiones de punto fijo para Aiken
Primera iteración
(𝑥𝑛+1 − 𝑥𝑛 )2
𝑥1 ≈ 𝑥0 −
𝑥𝑛+2 − 2𝑥1 + 𝑥0
(0,764102849⁡⁡ − 2)2
𝑥1 ≈ 𝑥0 −
1,312251263 − 2(0,764102849) + 2
𝑥1 = 1,1438695
Segunda iteración
(𝑥2 − 𝑥1 )2
𝑥2 = 𝑥1 −
𝑥3 − 2𝑥2 + 𝑥1
(1,312251263⁡⁡ − 0,764102849)2
𝑥2 = 0,764102849⁡⁡ −
1,120568715 − 2(1,312251263) + 0,764102849
𝑥2 = 1,170231365
Tabla de interacción del Método de Aitken
n Punto Fijo Aitken F(Xn) Ea(%)
0 2
1 0,764102849 1,143832535 -0,1057562063 100,0000%
2 1,312251263 1,17023168 -0,0204961697 2,2559%
3 1,120568715 1,175713838 -0,0025812939 0,4663%
4 1,197986008 1,176374071 -0,0004189142 0,0561%
5 1,16800434 1,176482712 -0,0000629957 0,0092%
6 1,179825585 1,176498954 -0,0000097841 0,0014%
7 1,175196189 1,176501482 -0,0000015011 0,0002%
8 1,177014037 1,176501869 -0,0000002314 0,0000%

Método de steffense
Primera Iteración
𝑥0 = 2 𝑥2 = 1,312251263
𝑥1 = 0,764102849
(𝑥1 − 𝑥0 )2
𝑥3 = 𝑥0 −
𝑥2 − 2𝑥1 + 𝑥0
(0,764102849⁡⁡ − 2)2
𝑥3 = 2 −
1,312251263 − 2(0,764102849⁡⁡) + 2
𝑥3 = 1,143832536
Segunda Iteración
𝑥4 = 𝑔(𝑥3 ) = 1,189163183
𝑥5 = 𝑔(𝑥4 ) = 1,1715104
(𝑥4 − 𝑥3 )2
𝑥6 = 𝑥3 −
𝑥5 − 2𝑥4 + 𝑥3
(1,189163183 − 1,143832536)2
𝑥6 = 1,143832536 −
1,1715104 − 2(1,189163183) + 1,143832536
𝑥6 = 1,176458062
Tabla de interacción del Método de Steffensen
n Xn Xn+1 Xn+2 X
0 2 0,764102849 1,312251263 1,143832535
3 1,143832535 1,189163183 1,171510400 1,176458062
6 1,176458062 1,176519156 1,176495185 1,176501940
9 1,176501940 1,176501940 1,176501940 1,176501940

OCTAVE
Graficamos…

PUNTO 8.
Método de Newton
𝑓(𝑥, 𝑦) = 12𝑥 2 − 12𝑦 2 − 24𝑥 + 7
𝑔(𝑥, 𝑦) = 9𝑥𝑦 − 9𝑦 + 1
Se halla las derivadas parciales
𝑓𝑥 = 24𝑥 − 24 ; 𝑓𝑦 = −24𝑦
𝑔𝑥 = 9𝑦⁡ ; 𝑔𝑦 = 9𝑥 − 9
24𝑥 − 24 −24𝑦 ∆𝑥 −12𝑥 2 + 12𝑦 2 + 24𝑥 − 7
[ ][ ] = [ ]
9𝑦 9𝑥 − 9 ∆𝑦 −9𝑥𝑦 + 9𝑦 − 1
Si tomamos (𝑥0 , 𝑦0 ) = (0,0)
24𝑥0 − 24 −24𝑦0 ∆𝑥0 −12𝑥0 2 + 12𝑦0 2 + 24𝑥0 − 7
[ ][ ]=[ ]
9𝑦0 9𝑥0 − 9 ∆𝑦0 −9𝑥0 𝑦0 + 9𝑦0 − 1
−24 0 ∆𝑥0 −7
[ ][ ] = [ ] nos queda un sistema
0 −9 ∆𝑦0 −1
−24∆𝑥0 = −7
−9∆𝑦0 = −1
7 1
Con lo que ∆𝑥0 = 24 y ∆𝑦0 = 9

7 7
𝑥1⁡=⁡ 𝑥0⁡ + ∆𝑥0 = 0 + =
24 24
 1 1
𝑦1 = 𝑦0 + ∆𝑦0 = ⁡0 + =
9 9
7 1
El nuevo punto de iteración es (𝑥1 , 𝑦1 ) = (24 ; 9)

24𝑥1 − 24 −24𝑦1 ∆𝑥1 −12𝑥1 2 + 12𝑦1 2 + 24𝑥1 − 7

⁡[ ][ ]=[ ]
9𝑦1 9𝑥1 − 9 ∆𝑦1 −9(𝑥1 𝑦1 + 9𝑦1 − 1
Tenemos entonces
8 377
−17 − ∆𝑥 −
[ 3 ] [ 1 ] = [ 432]
51 ∆𝑦1 7
1 − −
8 24
Queda el siguiente sistema
3 377
−17∆𝑥1 − 8 ∆𝑦1 = − 432 (1)
51 7
∆𝑥1 − ∆𝑦1 = − 24 (2)
8

Multiplicamos 2 por 17
3 377
−17∆𝑥1 − 8 ∆𝑦1 = − 432⁡
51 7
∆𝑥1 − ∆𝑦1 = − ⁡⁡⁡⁡⁡⁡
8 24
2665 2519
− 24 ∆𝑦1 = − 432

2519
∆𝑦1 = − = 0,0525198666
432
7 51 2519
∆𝑥1 = − + ( )
24 8 47970
∆𝑥1 = 0,04304724828
Con que
∆𝑥1 = 0,04304724828
∆𝑦1 = 0,0525198666
7
𝑥2 = 𝑥1 + ∆𝑥1 = 24 + 0,04304724828

𝑥2 = 0,33476339149
1
𝑦2 = 𝑦1 + ∆𝑦1 = 8+0,0525198666

𝑦2 = 0,1636230978
 𝑘𝑛 0 1 2

𝑥𝑘 0 0,291666667 0,3347639149

𝑦𝑘 0 0,1111111111 0,1636230978

Si tomamos ⁡⁡(𝑥0 , 𝑦0 ) = (2, −1)

24𝑥0 − 24 −24𝑦0 ∆𝑥0 −12𝑥0 2 + 12𝑦0 2 + 24𝑥0 − 7
[ ][ ]=[ ]
9𝑦0 9𝑥0 − 9 ∆𝑦0 −9𝑥0 𝑦0 + 9𝑦0 − 1
Tenemos que
24 24 ∆𝑥0 5
[ ][ ]=[ ]
−9 9 ∆𝑦0 8
Nos queda el siguiente sistema
24∆𝑥0 + 24∆𝑦0 = 5 por 9
−9∆𝑥0 + 9∆𝑦0 = 8 por 24
Tenemos que
216∆𝑥0 + 216∆𝑦0 = 45
−216∆𝑥0 + 216∆𝑦0 = 192⁡
216∆𝑥0 +216∆𝑦0 =45
−216∆𝑥0 +216∆𝑦0 =192⁡
⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡432∆𝑦0 =237

237
∆𝑦0 =
432
79
∆𝑦0 =
144
8
∆𝑥0 = − + ∆𝑦0
9
8 79
∆𝑥0 = − +
9 144
−49
∆𝑥0 =
144
−49 79
Con lo que ⁡⁡∆𝑥0 = ⁡ 144 y ⁡⁡∆𝑦0 = ⁡ 144
49 239
𝑥1 = 𝑥0 + ∆𝑥0 = 2 − 144⁡⁡=⁡⁡⁡144
 79 65
𝑦1 = 𝑦0 + ∆𝑦0 = −1 + ⁡ −= ⁡
144 144
234 −65
El nuevo punto de iteración es (𝑥1 , 𝑦1 )= 144 , 144

24𝑥1 − 24 −24𝑦1 ∆𝑥1 −12𝑥1 2 + 12𝑦1 2 + 24𝑥1 − 7

[ ][ ]=[ ]
9𝑦1 9𝑥1 − 9 ∆𝑦1 −9𝑥1 𝑦1 + 9𝑦1 − 1
95 65
20
[ 6 6 ] [∆𝑥1 ] = [ ]
65 95 ∆𝑦1 9
− 1,680121523
16 16
Nos queda el siguiente sistema
95 65 20 6
∆𝑥1 + ∆𝑦1 = ⁡⁡⁡⁡⁡⁡𝑝𝑜𝑟⁡
6 6 9 95
65 95 16
− ∆𝑥1 + ∆𝑦1 = 1,680121523⁡⁡⁡⁡𝑝𝑜𝑟⁡ ⁡
16 16 55
Tenemos
13 8
∆𝑥1 + ⁡⁡∆𝑦1 = ⁡
19 57
19
−∆𝑥1 + ⁡⁡∆𝑦1 = ⁡0,4135683761
13
13 8
∆𝑥1 + 19 ⁡⁡∆𝑦1 ⁡⁡ = ⁡⁡⁡⁡
57
19
−∆𝑥1 + 13 ⁡⁡∆𝑦1 = ⁡0,4135683761
⁡
530
⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡247 ∆𝑦1 = 0,5539192533

∆𝑦1 = 0,2581472746
8 13
∆𝑥1 = − (0,2581472746
57 19
∆𝑥1 = −0,03627620556
Con que ∆𝑥1 = −0,03627620556 y ∆𝑦1 = 0,2581472746
239
𝑥2 = 𝑥1 + ∆𝑥1 = − 0,03627620546⁡
144
𝑥2 = 1,623446017
−65
𝑦2 = 𝑦1 + ∆𝑦1 = ⁡ + 0,2581472746
144
𝑦2 = −0,1932416143

𝑘𝑛 0 1 2

𝑥𝑘 1 1,659722222 1,623446017

𝑦𝑘 2 -0,4513888889 −0,1932416143

PUNTO 9
1 𝑝𝑄𝑥
𝐹=
4𝜋𝑒0 2 2 32
(𝑥 𝑎 )
𝐷𝑎𝑡𝑜𝑠
𝐹 = 1𝑁
𝑞 = 2𝑥10−5
𝑎 = 0,8𝑚

𝑒0 = 8,85𝑥10−12 /𝑁𝑚2
𝑅𝑒𝑚𝑝𝑙𝑎𝑧𝑎𝑚𝑜𝑠
1 (2𝑥10−5 )(2𝑥10−5 )𝑥
𝐹= ( )
4𝜋(8,85𝑥10−12 ) √(𝑥 2 + (0,8)2 )3
𝑥
0 = 3,596721878 ( )−1
√(𝑥 2 + 0,64)3
𝑥
𝑓(𝑥) = 3,596721878 ( )−1
√(𝑥 2 + 0,64)3
Analizaremos la función en el intervalo [0,1], ya que gráficamente se observa una raíz,
usaremos el método de bisección para hallar la raíz deseada.
Método de bisección
𝑥
𝑎 = 0⁡⁡⁡⁡⁡⁡𝑓(0) = 3,596721878 ( )−1
√(02 + 0,64)3
𝑓(0) = −1
𝑥
𝑏 = 1⁡⁡⁡⁡⁡⁡𝑓(1) = 3,596721878 ( )−1
√(12 + 0,64)3
𝑓(1) = 0,7125414212
Calculamos el número de iteraciones con una exactitud de 10−5
𝑏 − 𝑎)
log⁡(
𝑛≥ 𝑒 )
log⁡(2)
1 − 0)
log⁡( )
𝑛≥ 105
log⁡(2)
log⁡(100000)
𝑛≥
log⁡(2)
𝑛 ≥ 16,60964047
𝑛 = 17
Primera iteración
𝑎1 + 𝑏1 0 + 1
𝑝1 = = = 0,5
2 2
𝑝1 = 0,5
0,5
𝑓(0,5) = 3,596721878 ( )−1
√((0,5)2 + 0,64)3
𝑓(0,5) = 1,141863801
Como 𝑓(𝑝1 )⁡𝑡𝑖𝑒𝑛𝑒⁡𝑒𝑙⁡𝑚𝑖𝑠𝑚𝑜⁡𝑠𝑖𝑔𝑛𝑜⁡𝑞𝑢𝑒⁡𝑓(𝑏1 )⁡𝑒𝑛𝑡𝑜𝑛𝑐𝑒𝑠⁡⁡
𝑎1 = 𝑎1 y 𝑏1 = 𝑝1 esta seria así
𝑎1=0 ⁡⁡⁡⁡⁡⁡𝑦⁡⁡⁡𝑏2 = 0,5
Segunda iteración
𝑎2 + 𝑏2 0 + 0,5
𝑝2 = = = 0,25
2 2
𝑝2 = 0,25
𝑓(0,25) = 0,527134795
Como 𝑓(𝑝2 )⁡⁡⁡𝑡𝑖𝑒𝑛𝑒⁡𝑒𝑙⁡𝑚𝑖𝑠𝑚𝑜⁡𝑠𝑖𝑔𝑛𝑜⁡𝑞𝑢𝑒⁡𝑓(𝑏2 )⁡𝑒𝑛𝑡𝑜𝑛𝑐𝑒𝑠⁡
𝑎3 = 𝑎2 y 𝑏3 = 𝑏2 . Está así 𝑎3 = 0 y 𝑏3 = 0,25
Continuamos las iteraciones en la tabla

n an(-) bn(+) pn F(pn) Ea(%)

0 0 1 0,5 1,1418638009
1 0 0,5 0,25 0,5271347950 100%
2 0 0,25 0,125 -0,1530971169 100%
3 0,125 0,25 0,1875 0,2156307863 33%
4 0,125 0,1875 0,15625 0,0376924547 20%
5 0,125 0,15625 0,140625 -0,0562103218 11,1111%
6 0,140625 0,15625 0,1484375 -0,0088706189 5,2632%
7 0,1484375 0,15625 0,15234375 0,0145098236 2,5641%
8 0,1484375 0,15234375 0,150390625 0,0028441043 1,2987%
9 0,1484375 0,150390625 0,149414063 -0,0030071602 0,6536%
10 0,149414063 0,150390625 0,149902344 -0,0000800001 0,3257%
11 0,149902344 0,150390625 0,150146484 0,0013824345 0,1626%
12 0,149902344 0,150146484 0,150024414 0,0006513128 0,0814%
13 0,149902344 0,150024414 0,149963379 0,0002856802 0,0407%
14 0,149902344 0,149963379 0,149932861 0,0001028460 0,0204%
15 0,149902344 0,149932861 0,149917603 0,0000114245 0,0102%
16 0,149902344 0,149917603 0,149909973 -0,0000342874 0,0051%
17 0,149909973 0,149917603 0,149913788 -0,0000114314 0,0025%
18 0,149913788 0,149917603 0,149915695 -0,0000000034 0,0013%
19 0,149915695 0,149917603 0,149916649 0,0000057105 0,0006%
20 0,149915695 0,149916649 0,149916172 0,0000028535 0,0003%
21 0,149915695 0,149916172 0,149915934 0,0000014250 0,0002%
22 0,149915695 0,149915934 0,149915814 0,0000007108 0,0001%
23 0,149915695 0,149915814 0,149915755 0,0000003537 0,0000%
B) El desplazamiento de una estructura está definido por la ecuación siguiente para una oscilación
amortiguada:

𝑦 = 9𝑒 −𝑘𝑡 cos⁡(𝑤𝑡)
donde 𝑘 = 0,7⁡⁡⁡𝑦⁡⁡⁡𝑤 = 4
1) Utilice el método gráfico para realizar una estimación inicial del tiempo que se requiere
para que el desplazamiento disminuya a 3,5.

Solución:
𝑦 = 9𝑒 −𝑘𝑡 cos⁡(𝑤𝑡)
Datos
𝑘 = 0,7⁡⁡⁡
𝑤=4
𝑦 = 3,5
Remplazamos

3,5 = 9𝑒 −0,7𝑡 cos⁡(4𝑡)

0 = 9𝑒 −0,7𝑡 cos(4𝑡) − 3,5
𝑓(𝑡) = 9𝑒 −0,7𝑡 cos(4𝑡) − 3,5
2) analizamos la función en el intervalo [0,1], ya que gráficamente se observa una raíz,
usaremos el método de Newton-Raphson con 𝑥0 = 0,5⁡𝑝𝑎𝑟𝑎⁡ℎ𝑎𝑙𝑙𝑎𝑟⁡𝑙𝑎⁡𝑟𝑎𝑖𝑧⁡𝑑𝑒𝑠𝑒𝑎𝑑𝑎.

Derivemos 𝑓(𝑡)
𝑓(𝑡) = 6,3𝑒 −0,7𝑡 cos(4𝑡) + 36𝑒 −0,7𝑡 𝑠𝑒𝑛(4𝑡)
𝑓(𝑡𝑛 )
𝑡𝑛+1 = 𝑡𝑛 −
𝑓′(𝑡𝑛 )
Primero iteración
𝑓(𝑡0 )
𝑡1 = 𝑡0 −
𝑓′(𝑡0 )
𝑓(0,5)
𝑡1 = 0,5 −
𝑓′(0,5)
(−6,1392834736)
𝑡1 = 0,5 −
(−24,9152568326)
𝑡1 = 0,2535934101
Segunda iteración
𝑓(𝑡1 )
𝑡2 = 𝑡1 −
𝑓′(𝑡1 )
𝑓(0,2535934101)
𝑡2 = 0,2535934101 −
𝑓′(0,2535934101)
(0,4802230444)
𝑡2 = 0,2535934101 −
(−22,8110704774)
𝑡2 = 0,27464560

Tabla de interacción del método de newton-Raphson

n Xn F(Xn) F´(Xn) Ea(%)
0 0,5 -6,13928347 -24,91525683
1 0,25359341 0,48022304 -22,81107048 97,1660%
2 0,27464560 -0,12226648 -24,08849443 7,6652%
3 0,26956988 0,02379358 -23,79966722 1,8829%
4 0,27056963 -0,00494090 -23,85752289 0,3695%
5 0,27036253 0,00101299 -23,84557679 0,0766%
6 0,27040501 -0,00020823 -23,84802887 0,0157%
7 0,27039627 0,00004278 -23,84752494 0,0032%
8 0,27039807 -0,00000879 -23,84762848 0,0007%
9 0,27039770 0,00000181 -23,84760720 0,0001%
10 0,27039778 -0,00000037 -23,84761157 0,0000%

2ℎ ℎ
C) t = √ 𝑔 + 𝑐

2ℎ ℎ
1. t = √ 𝑔 + 𝑐

ℎ 2ℎ
t - ⁡𝑐 = √ 𝑔

ℎ 2 2ℎ
(t⁡ − ⁡⁡⁡ 𝑐 ⁡⁡) = √ 𝑔

𝑡ℎ ℎ2 2ℎ
t2 - ⁡2 + 𝑐2 =
𝑐 𝑔

𝑡ℎ ℎ2 2ℎ
t2 - ⁡2 + 𝑐2 - = 0 X c2
𝑐 𝑔

2ℎ𝑐 2
c2 t2 – 2cth + h2 - =0
𝑔

𝑐2
h2 – 2 (tc⁡ − ⁡⁡⁡ 𝑔 ⁡⁡) ℎ + ⁡ 𝑐 2 𝑡 2 = 0

2.
1
h = 2 g t1 2
h = ct2
t1 + t2 = 10
t2 = 10 – t1
Se iguala…
1
g t12 = ct2
2
1
g t12 = c (10 – t1)
2
1
g t12 = 10c – ct1
2
1
g t12 + ct1 – 10c = 0
2
1
f(t) = 2 g t1 + ct1 – 10c
1
f(t) = 2 (9,8) t12 + 343 t1 - 3430

f(t) = 4,9 t12 + 343 t1 – 3430

Para X0 = 5
Derivamos f(t)
f’(t) = 9,8t + 343
Aplicamos la ecuación:
𝑓(Xn)
Xn+1 = Xn - 𝑓′ (Xn)

Primera iteración
𝑓(𝑋o)
X1 = X0 + 𝑓′ (Xo)

𝑓(5)
X1 = 5 + 𝑓′ (5)

(−1532,5)
X1 = 5 - 392

X1 = 9,0625
Segunda iteración
𝑓(X1⁡)
X2 = X1 - 𝑓′ (X1⁡)

𝑓(9,0625)
X2 = 9,0625 - 𝑓′(9,0625)
 80,8691
X2 = 9,0625 - 431,8125

X2 = 8,875221531
Continuamos con las interacciones en la tabla

n Xn F(Xn) F´(Xn) Ea(%)

0 5 -1592,50000000 392,00000000

1 9,06250000 80,86914063 431,81250000 44,8276%

2 8,87522163 0,17185862 429,97717199 2,1101%

3 8,87482194 0,00000078 429,97325500 0,0045%

4 8,87482194 0,00000000 429,97325498 0,0000%

D) 𝜋⁡𝑟 2 ℎ = 1000
Ahora la superficie del material usado:
S =⁡2𝜋 (r + 0,25)2 + (2𝜋r + 0,25) h
Despejamos h
1000
h= 𝜋𝑟 2

reemplazamos
1000
S = 2𝜋(r + 0,25)2 + (2𝜋r + 0,25) ( 𝜋𝑟 2 )
2000𝜋𝑟⁡+⁡250
S = 2𝜋(r + 0,25)2 + ( )
𝜋𝑟 2

Derivamos la función en r y lo igualamos a 0

(2000𝜋𝑟)(𝜋𝑟 2 )−(⁡2000𝜋𝑟⁡+⁡250)(2𝜋𝑟)
S’ = 4𝜋(r + 0,25) + (𝜋𝑟 2 )2

2000𝜋 2 𝑟 2 −⁡4000𝜋2 𝑟 2 −⁡500𝜋𝑟

S’ = 4𝜋(r + 0,25) + 𝜋2𝑟 4

2000𝜋 2 𝑟 2 −⁡4000𝜋2 𝑟 2 −⁡500𝜋𝑟

S’ = 4𝜋(r + 0,25) + 𝜋2𝑟 4

2000 500
S’ = 4𝜋r + 1 - ( ) − (𝜋𝑟 3) = 0
𝑟2

2000 500
f(X) = 4𝜋𝑋 - ( ) − (𝜋𝑋 3) + 1
𝑋2

Analizamos la función en el intervalo [5, 6] ya que gráficamente se observa una raíz

usaremos el método de la secante para hallar dicha raíz deseada.
 Método de la secante
f(X0) = f (5) = -17,44138647
f(X1) = f (6) = 20,10583969
Primera iteración
𝑓(X1)(X1−⁡X0)
X2 = X1 - ( 𝑓(X1)−𝑓(X0) )

(20,10583969)(6−⁡5)
X2 = 6 - 20,10583969⁡+⁡17,44138647

X2 = 5,464443676
Segunda iteración
f(X1) = f (6) = 20,10583969
f(X2) = f (5,464443676) = 1,713913622
𝑓(X2)(X2−⁡X1)
X3= X2 - ( 𝑓(X2)−𝑓(X1) )

(1,713913622)(5,464443676−⁡6)
X3 = 5,464443676 - (1,713913622)−(20,10583969)

X3 = 5,414536055
Luego…
f(X3) = 0,1808904941
Continuamos en la siguiente tabla.
i Xi-1 F(Xi-1) Xi F(Xi) Xi+1 F(Xi+1) Ea(%)
1 5 -17,44138647 6 20,10583969 5,464518641 10,3046133768
2 6 20,10583969 5,464518641 1,7167334985 5,4145282404 9,9699173486 0,9233%
3 5,464518641 1,71673350 5,41452824 -0,1811899408 5,4193006972 10,0016504737 0,0881%
4 5,41452824 -0,18118994 5,419300697 0,0015299374 5,4192607368 10,0013845757 0,0007%
5 5,419300697 0,00152994 5,419260737 0,0000013558 5,4192607014 10,0013843399 0,0000%

OCTAVE
A.
Graficamos…
B.

Graficamos…
C.
D.
 ANEXOS
A continuación, se anexa los códigos o algoritmos que se usaron para el
desarrollo del taller dependiendo a cada método usado para la solución de las
distintas ecuaciones y problemas.