Desafío 5: La velocidad ascendente, v , de un cohete puede ser calculada

por la siguiente expresión:

m0
v =u∗ln ( m0−qt)−¿

En dónde u es la velocidad relativa a la que el combustible es expedido,

m0 es la masa inicial del cohete en el instante t=0 , q es la tasa de
consumo de combustible y g es la aceleración de la gravedad.
2 5
Considerando u=2200 m/ s , g=9.81m/s , m0=1.6 ×10 kg y
q=2680 kg /s , calcule, utilizando los métodos indicados en el Ejercicio 2,
el tiempo para el cual el cohete alcanza la velocidad de v =1000 m/ s ,
sabiendo que ese instante está entre 20 s y 30 s.

Ejercicio 2:

Determine, con una exactitud de 10−6 o un máximo de 10 iteraciones, la raíz

visualizada en el ejercicio anterior empleando cada uno de los siguientes
métodos:

 Bisección
 Regula Falsi https:
 Newton – Raphson
 Secante

Para cada método haga una descripción teórica (las ecuaciones deberán ir en
un editor de texto) y elabore una tabla con los resultados. Debe contener al
menos la siguiente información por columna: El número de iteración niter
(empezando en cero), valor aproximado de la raíz en cada iteración x i , el
x
valor de la función evaluada en la raíz aproximada en esa interacción, f (¿¿ i)
¿
, y el error relativo, Erel ( ) . Los cálculos deberán ser consignados en tablas y
no es necesario que hagan cada cálculo en un editor de ecuaciones, para eso
son los programas de cálculo.
x
Note que, si el método converge, en cada iteración i el valor de f (¿¿ i) se
¿
debe ir aproximando cada vez más a cero. Si esto no ocurre revise
cuidadosamente sus cálculos. Realice los cálculos empleando 6 cifras
decimales. Escoja los valores iniciales de forma adecuada y que tengan sentido
con el enunciado del ejercicio.

Bisección

Dada una función en la que en los puntos ubicados entra a y b, existe

una solución de la ecuación. Así mismo en ese punto la ecuación es
continua. Se debe encontrar un punto en que la función vale 0. La mitad
que genera los intervalos a y b, se crea el nuevo punto c.
 Bisección
1400

1300

1200 1179
1242
1118
1100 1059

1000 948
1003
896
900
845
795
800 748
700 656
613701
600 571
529
489
500 451
413
400 376
340
305
300 271
238
206
200 174
143
113
100 5584
27
0
-27
0
-40 -35 -30 -25 -20 -15 -10 -53 0
-5-78 5 10 15 20 25 30 35
-103
-128-100
-152
-175
-198
-200

-300

Segundo Velocida
s d
x f(x)
-30 -601
-29 -587
-28 -571
-27 -556
-26 -540
-25 -524
-24 -508
-23 -491
-22 -474
-21 -457
-20 -439
-19 -421
-18 -403
-17 -384
-16 -365
-15 -346
-14 -326
-13 -306
-12 -285
-11 -264
-10 -243
-9 -221
-8 -198
-7 -175
-6 -152
-5 -128
-4 -103
-3 -78
-2 -53
-1 -27
0 0
1 27
2 55
3 84
4 113
5 143
6 174
7 206
8 238
9 271
10 305
11 340
12 376
13 413
14 451
15 489
16 529
17 571
18 613
19 656
20 701
21 748
22 795
23 845
 24 896
25 948
26 1003
27 1059
28 1118
29 1179
30 1242

Bisección
Función: $G$3*LN(($G$4)/(($G$4)-($G$6*$C$4)))-($G$5*$C$4)
Xi: 20
Xu: 30
Tolerancia 1,00E-06
Iteraciones 10

u=Velocidad Relativa 2200

g=Gravedad 9.81
m0=masa inicial del cohete 160000
q
¿ Tasade consumo de combustible 2680

Error
Iteraciones a b Xr f(a) f(b) f(Xr) xa*xr Convergencia
(%)
30,0
0 20,00 0 25,00 701,3301243 1241,651389 948,413509 665.150,963875
30,0 1088,38227 9,09090
1 25,00 0 27,50 948,4135086 1241,651389 8 1.032.236,455 9 El método no converge
30,0 1163,22125 4,34782
2 27,50 0 28,75 1088,382278 1241,651389 7 1.266.029,401 6 El método no converge
30,0 1201,96901 2,12766
3 28,75 0 29,38 1163,221257 1241,651389 8 1.398.155,912 0 El método no converge
30,0 1221,69095 1,05263
4 29,38 0 29,69 1201,969018 1241,651389 1 1.468.434,672 2 El método no converge
30,0 1231,64104 0,52356
5 29,69 0 29,84 1221,690951 1241,651389 5 1.504.684,719 0 El método no converge
30,0 1236,63864 0,26109
6 29,84 0 29,92 1231,641045 1241,651389 6 1.523.094,915 7 El método no converge
7 29,92 30,0 29,96 1236,638646 1241,651389 1239,14312 1.532.372,271 0,13037 El método no converge
 0 0 8
30,0 1240,39678 0,06514
8 29,96 0 29,98 1239,14312 1241,651389 0 1.537.029,136 7 El método no converge
30,0 1241,02396 0,03256
9 29,98 0 29,99 1240,39678 1241,651389 6 1.539.362,131 3 El método no converge
30,0 1241,33764 0,01627
10 29,99 0 30,00 1241,023966 1241,651389 8 1.540.529,770 9 El método no converge
30,0 1241,49451 0,00813
11 30,00 0 30,00 1241,337648 1241,651389 1 1.541.113,876 9 El método no converge
30,0 1241,57294 0,00406
12 30,00 0 30,00 1241,494511 1241,651389 8 1.541.406,000 9 El método no converge
30,0 1241,61216 0,00203
13 30,00 0 30,00 1241,572948 1241,651389 8 1.541.552,080 5 El método no converge
30,0 1241,63177 0,00101
14 30,00 0 30,00 1241,612168 1241,651389 9 1.541.625,125 7 El método no converge
30,0 1241,64158 0,00050
15 30,00 0 30,00 1241,631779 1241,651389 4 1.541.661,648 9 El método no converge
30,0 1241,64648 0,00025
16 30,00 0 30,00 1241,641584 1241,651389 7 1.541.679,910 4 El método no converge
30,0 1241,64893 0,00012
17 30,00 0 30,00 1241,646487 1241,651389 8 1.541.689,041 7 El método no converge
30,0 1241,65016 0,00006
18 30,00 0 30,00 1241,648938 1241,651389 4 1.541.693,607 4 El método no converge
30,0 1241,65077 0,00003
19 30,00 0 30,00 1241,650164 1241,651389 6 1.541.695,890 2 El método no converge
30,0 1241,65108 0,00001
20 30,00 0 30,00 1241,650776 1241,651389 3 1.541.697,031 6 El método no converge
30,0 1241,65123 0,00000
21 30,00 0 30,00 1241,651083 1241,651389 6 1.541.697,602 8 El método no converge
30,0 1241,65131 0,00000
22 30,00 0 30,00 1241,651236 1241,651389 3 1.541.697,887 4 El método no converge
30,0 1241,65135 0,00000
23 30,00 0 30,00 1241,651313 1241,651389 1 1.541.698,030 2 El método no converge
30,0 1241,65137 0,00000
24 30,00 0 30,00 1241,651351 1241,651389 0 1.541.698,101 1 El método converge
30,0 1241,65138 0,00000
25 30,00 0 30,00 1241,65137 1241,651389 0 1.541.698,137 0 El método converge
30,0 1241,65138 0,00000
26 30,00 0 30,00 1241,65138 1241,651389 4 1.541.698,155 0 El método converge

Regla Falsa
Función: $G$3*LN(($G$4)/(($G$4)-($G$6*$C$4)))-($G$5*$C$4)
 Xi: 20
Xu: 30
Tolerancia 1,00E-06
Iteraciones 10

u=Velocidad Relativa 2200

g=Gravedad 9.81
m0=masa inicial del cohete 160000
q ¿ Tasade consumo de combustible 2680

Iteración Xm f(xm) xi xf f(xi) f(xf) f(xf)*f(xm) Error Convergencia

206,342053 1241,65138
1 7,020 5 20,00 30,00 701,330124 9 256204,897344
561,978984 206,34205 4001,19584 2248587,97541
2 16,795 3 7,020127 3 206,342053 2 2 58,20 No converge
2060,84923 16,79508 206,34205 4001,19584 8245861,39602
3 40,139 6 3 3 561,978984 2 5 58,16 No converge
108,107077 40,13878 206,34205 2060,84923 4001,19584
4 96,641 4 9 3 6 2 -432557,588796 58,47 No converge
85,7009619 96,64115 206,34205 4001,19584
5 99,527 1 2 3 108,107077 2 -342906,332444 2,90 No converge
70,3289555 99,52715 206,34205 4001,19584
6 97,189 3 1 3 -85,700962 2 -281399,924432 2,41 No converge
70,3289555 99,52715 206,34205 4001,19584
7 97,189 3 1 3 -85,700962 2 -281399,924432 0,00 Converge
70,3289555 99,52715 206,34205 4001,19584
8 97,189 3 1 3 -85,700962 2 -281399,924432 0,00 Converge
70,3289555 99,52715 206,34205 4001,19584
9 97,189 3 1 3 -85,700962 2 -281399,924432 0,00 Converge
70,3289555 99,52715 206,34205 4001,19584
10 97,189 3 1 3 -85,700962 2 -281399,924432 0,00 Converge

El error debe ser menor al 1%, para hallar la raíz del método y así se podrá
saber en cual iteración el método converge.
 Newton – Raphson
Función: $G$3*LN(($G$4)/(($G$4)-($G$6*$C$4)))-($G$5*$C$4)
Xi: 20
Xu: 30
Tolerancia 1,00E-06
Iteraciones 10

u=Velocidad Relativa 2200

g=Gravedad 9.81
m0=masa inicial del cohete 160000
q ¿ Tasade consumo de combustible 2680

Iteración Xi f(Xi) f´(Xi) Error Convergencia

20,00
701,3301243 45,60353383
1 0
2 4,621 131,907456 30,13165083 332,79 No converge
3 0,243 6,600983433 27,19087611 1798,26 No converge
4 0,001 0,018287826 27,04041745 35894,95 No converge
5 0,000 1,41164E-07 27,04 12954760,39 No converge
6 0,000 4,01189E-13 27,04 58658293,49 No converge
7 0,000 -4,30257E-13 27,04 249,91 No converge
8 0,000 3,90644E-13 27,04 159,52 No converge
9 0,000 4,38694E-14 27,04 323,06 No converge
10 0,000 -4,28713E-13 27,04 26,62 No converge

Iteración Xi f(Xi) f´(Xi) Error Convergencia

60,00
11067,69821 -7379,81
1 0
61,50
7102,301751 1233,23241 2,44
2 0 No converge
67,25
3887,237363 300,918162 8,56
3 9 No converge
80,17
1567,774495 117,256793 16,11
4 7 No converge
93,54
330,9021517 74,8108967 14,29
5 7 No converge
97,97
17,2954165 67,2979819 4,51
6 0 No converge
 98,22
0,049387891 66,9144923 0,26
7 7 No converge
98,22
4,03723E-07 66,9133983 0,00
8 8 No converge
98,22
0 66,9133983 0,00
9 8 Converge
98,22
0 66,9133983 0,00
10 8 Converge

Cuando el método inicia partiendo del rango de 20, y mostrando 10 iteraciones

queda demostrado que no converge, pero solo hasta el rango inicia en 60, en
la novena iteración converge. De resto en ninguna otra ocasión iniciando en el
rango de 20 hasta 60, no converge el método.

Secante
$G$3*LN(($G$4)/(($G$4)-($G$6*$C$4)))-
Función: ($G$5*$C$4)
Xi: 20
Xu: 30
Tolerancia 1,00E-06
Iteraciones 10

u=Velocidad Relativa 2200

g=Gravedad 9.81
m0=masa inicial del cohete 160000
q ¿ Tasade consumo de combustible 2680

Iteración Xi f(Xi) Error Convergencia

-1 20,000 701,3301243
0 21,000 747,64348 4,76 No converge
1 4,856849 139,0298106 332,38 No converge
2 1,169158 32,04148988 315,41 No converge
3 0,064747 1,752051379 1705,74 No converge
4 0,000864 0,023355909 7396,06 No converge
5 6,38E-07 1,72592E-05 135223,21 No converge
6 6,31E-12 1,70647E-10 10119096,29 No converge
7 -3,3E-15 3,202E-14 193348,22 No converge
8 -4,4E-15 4,36389E-14 26,63 No converge
9 0 0 #¡DIV/0! #¡DIV/0!
10 0 0 #¡DIV/0! #¡DIV/0!

Iteración Xi f(Xi) Error Convergencia

 -1 60,000 11067,69821
0 61,000 7823,501347 1,64 No converge
1 63,41154 5490,221664 3,80 No converge
2 69,0859 3392,943987 8,21 No converge
3 78,2658 1802,068738 11,73 No converge
4 88,66435 721,5567113 11,73 No converge
5 95,60842 180,6204143 7,26 No converge
6 97,92707 20,21051113 2,37 No converge
7 98,2192 0,595444184 0,30 No converge
8 98,22807 0,001988878 0,01 No converge
9 98,2281 1,9601E-07 0,00 No converge
10 98,2281 -9,09495E-13 0,00 Converge

Como el método anterior al iniciar el rango partiendo de 20, y mostrando 10

iteraciones queda demostrado que no converge, pero solo hasta el rango inicia
en 60, en la décima iteración converge. De resto en ninguna otra ocasión
iniciando en el rango de 20 hasta 60, no converge el método.