Universidad Católica Nordestana

(Ucne)

Asignatura:
Metodos Numéricos.

Maestro:
Juan Jose bello

Tema:
Tarea numero 3

Alumno:
Anneury Sosa Abreu. . . .2018-0906.
 1. expansión en serie de Mac-Laurin para cos x es:

Iniciando con el primer término cos x = 1, agregue los términos uno a uno para estimar
cos(π/3). Después de que agregue cada uno de los términos, calcule los errores relativos
porcentuales verdaderos y aproximados. Use una calculadora para determinar el valor
exacto. Agregue términos hasta que el valor absoluto del error aproximado se encuentre
dentro de cierto criterio de error, considerando tres cifras significativas.

𝑉𝑎𝑙𝑜𝑟 𝑒𝑥𝑎𝑐𝑡𝑜: 𝑐𝑜𝑠(𝜋/3) ≈ 1

(𝜋/3)2
𝑐𝑜𝑠(𝜋/3) ≈ 1 − = 0.452
2!
|0.452 − 0.5|
𝜀𝑎 = = 0.106 = 10.6%
0.452

(𝜋/3)2 4
𝑐𝑜𝑠(𝜋/3) ≈ 1 − + (𝜋/3) = 0.502
2! 4!
|0.502 − 0.412|
𝜀𝑎 = = 0.179 = 17.9%
0.502

(𝜋/3)2 (𝜋/3)4 6
𝑐𝑜𝑠(𝜋/3) ≈ 1 − + − (𝜋/3) = 0.499
2! 4! 6!
|0.499 − 0.502|
𝜀𝑎 = = 0.006 = 0.6%
0.499

(𝜋/3)2 (𝜋/3)4 (𝜋/3)6 8

𝑐𝑜𝑠(𝜋/3) ≈ 1 − + − + (𝜋/3) = 0.500
2! 4! 6! 8!
|0.500 − 0.499|
𝜀𝑎 = = 0.002 = 0.2%
0.5
2. Repita los cálculos del problema anterior, pero ahora usando la expansión de la serie de
Mac-Laurin para sen x, para evaluar el sen (π/3).

𝑠𝑒𝑛 (𝜋/3) ≈ 0.866

(𝜋/3)3
𝑠𝑒𝑛(𝜋/3) ≈ 𝜋/3 − = 0.856
3!
|0.856 − 0.866|
𝜀𝑎 = = 0.0116 = 1.16%
0.856

(𝜋/3)3 5
𝑠𝑒𝑛(𝜋/3) ≈ 𝜋/3 − + (𝜋/3) = 0.866
3! 5!
|0.866 − 0.856|
𝜀𝑎 = = 0.0115 = 1.15%
0.866

(𝜋/3)3 (𝜋/3)5 − (𝜋/3)

7

𝑠𝑒𝑛(𝜋/3) ≈ 𝜋/3 − + = 0.865

3! 5! 7!
|0.865 − 0.866|
𝜀𝑎 = = 0.001 = 0.1%
0.865
3. Emplee la expansión de la serie de Taylor de cero hasta tercer orden para predecir f(3) si

usando como punto base x = 1. Calcule el error relativo porcentual verdadero et

para cada aproximación

𝑓(3) = 25(3)3 − 6(3)2 + 7(3) − 88 = 554

ℎ = 𝑥𝑖+1 − 𝑥𝑖 = 3 − 1 = 2

𝑓(3) ≈ 𝑓(1) = 25(1)3 − 6(1)2 + 7(1) − 88 = −62

554 − (−62)
𝜀𝑡 = = 1.112 = 111.2%
554

𝑓′(𝑥) = 75𝑥2 − 12𝑥 + 7

𝑓′(1) = 75(1)2 − 12(1) + 7 = 70
𝑓(3) = 𝑓(1) + 𝑓′(1)ℎ = −62 + 70 ∗ 2 = 78
554 − 78
𝜀𝑡 = = 0.859 = 85.9%
554

𝑓′′(𝑥) = 150𝑥 − 12
𝑓′′(1) = 150(1) − 12 = 138
𝑓′′(1) 2 138 2
𝑓(3) = 𝑓(1) + 𝑓′(1)ℎ + ℎ = −62 + 70 ∗ 2 + ∗ 2 = 354
2! 2!
554 − 354
𝜀𝑡 = = 0.361 = 36.1%
554

𝑓(3)(𝑥) = 150
𝑓(3)(1) = 150
(3)
𝑓′′(1) 2 𝑓 (1) ℎ3 = −62 + 70 ∗ 2 + 138 150
𝑓(3) = 𝑓(1) + 𝑓′(1)ℎ + ℎ + ∗ 22 + ∗ 23
2! 3! 2! 3!
= 554
554 − 554
𝜀𝑡 = = 0 = 0%
554
4. Use la expansión de la serie de Taylor de cero al cuarto orden para estimar f(2.5)
si f(x) = ln x utilizando x = 1 como punto base. Calcule el error relativo porcentual
verdadero et para cada aproximación. Analice los resultados. Repita el ejemplo
usando Excel.

valor 0.91629073
exacto
x 1
h 1.5

iteración f(x) Error

relativo
0 0 100.00%
1 1.5 -63.70%
2 0.375 59.07%
3 1.5 -63.70%
4 0.234375 74.42%
5. Utilice aproximaciones en diferencias de O(h) hacia atrás y hacia adelante y una
aproximación de diferencia central de O(h2) para estimar la primera derivada de la
función mencionada en el problema
3. Evalúe la derivada en x = 2 usando un tamaño del incremento h = 0.2. Compare los
resultados con el valor exacto de las derivadas. Interprete los resultados considerando el
término residual de la expansión en la serie de Taylor. Repita el ejemplo usando Excel.

𝑓(𝑥) = 25𝑥3 − 6𝑥2 + 7𝑥 − 88

𝑓′(𝑥) = 75𝑥2 − 12𝑥 + 7

𝑓′(2) = 75(2)2 − 12(2) + 7 = 283

𝑥𝑖 − 1 = 1.8 𝑓(𝑥𝑖 − 1) = 25(1.8)3 − 6(1.8)2 + 7(1.8) − 88 = 50.96

𝑥𝑖 = 2 𝑓(𝑥𝑖 ) = 25(2)3 − 6(2)2 + 7(2) − 88 = 102

𝑥𝑖 + 1 = 2.2 𝑓(𝑥𝑖 + 1) = 25(2.2)3 − 6(2.2)2 + 7(2.2) − 88 = 164.56

Hacia
adelante:
𝑓(𝑥𝑖+1) − 𝑓(𝑥𝑖) 164.56 − 102
𝑓′(𝑥) ≈ = = 312.8
ℎ 0.2
𝐸𝑡 |283 − 312.8|
𝜀% = ∗ 100 = ∗ 100 = 10.53%
𝑉𝑎𝑙𝑜𝑟 𝑣𝑒𝑟𝑑𝑎𝑑𝑒𝑟𝑜 283

Hacia
atrás:
𝑓(𝑥𝑖) − 𝑓(𝑥𝑖−1) 102 − 50.96
𝑓′(𝑥) ≈ = = 255.2
ℎ 0.2
𝐸𝑡 |283 − 255.2|
𝜀% = ∗ 100 = ∗ 100 = 9.82%
𝑉𝑎𝑙𝑜𝑟 𝑣𝑒𝑟𝑑𝑎𝑑𝑒𝑟𝑜 283

Centrada
:
𝑓(𝑥𝑖+1) − 𝑓(𝑥𝑖−1) 164.56 − 50.96
𝑓′(𝑥) ≈ = = 284
2ℎ 2 ∗ 0.2
𝐸𝑡 |283 − 284|
𝜀% = ∗ 100 = ∗ 100 = 0.35%
𝑉𝑎𝑙𝑜𝑟 𝑣𝑒𝑟𝑑𝑎𝑑𝑒𝑟𝑜 283
En Excel:

valor verdadero 283

h 0.2

x-h 1.8
x 2
x+h 2.2

x-1 50.96
x 102
x+1 164.56

hacia adelante 312.8

error 10.53%

hacia atras 255.2

error 9.82%

centrada 284
error 0.35%
6. Con la aproximación en diferencias centrales de O(h2) estime la segunda derivada de
la función examinada en el problema 3. Realice la evaluación para x = 2 usando un
tamaño de incremento h =
0.25 y 0.125. Compare lo estimado con el valor exacto de la segunda derivada. Interprete
sus resultados considerando el término residual de la expansión en la serie de Taylor.

𝑓(𝑥) = 25𝑥3 − 6𝑥2 + 7𝑥 − 88

𝑓′′(𝑥) = 150𝑥 − 12

𝑓′′(2) = 150(2) − 12 = 288

h = 0.25

𝑥𝑖 − 1 = 1.75 𝑓(𝑥𝑖 − 1) = 25(1.75)3 − 6(1.75)2 + 7(1.75) − 88 = 39.859375

𝑥𝑖 = 2 𝑓(𝑥𝑖 ) = 25(2)3 − 6(2)2 + 7(2) − 88 = 102

𝑥𝑖 + 1 = 2.25 𝑓(𝑥𝑖 + 1) = 25(2.25)3 − 6(2.25)2 + 7(2.25) − 88 = 182.140625

𝑓(2.25) − 2 ∗ 𝑓(2) + 𝑓(1.75) 39.859375 − 204 + 182.140625

𝑓′′(2) ≈ = = 288
0.252 0.0625

h = 0.125

𝑥𝑖 − 1 = 1.875 𝑓(𝑥𝑖 − 1) = 25(1.875)3 − 6(1.875)2 + 7(1.875) − 88 = 68.82617187

𝑥𝑖 = 2 𝑓(𝑥𝑖 ) = 25(2)3 − 6(2)2 + 7(2) − 88 = 102

𝑥𝑖 + 1 = 2.125 𝑓(𝑥𝑖 + 1) = 25(2.125)3 − 6(2.125)2 + 7(2.125) − 88 = 139.6738281

𝑓(2.25) − 2 ∗ 𝑓(2) + 𝑓(1.75) 68.82617187 − 204 + 139.6738281

𝑓′′(2) ≈ = ≈ 288
0.1252 0.015625
7. Recuerde que la velocidad de caída del paracaidista puede calcularse con

Use un análisis de error de primer orden para estimar el error de v(t) para t = 6, si g =
9.81 m/s2 y m = 50 kg., pero c = 12.5 ± 1.5.

𝑔𝑚 −𝑐∗𝑡⁄𝑚
𝑣(𝑡) = (1 − 𝑒 )
𝑐
9.81 ∗ 50 −12.5∗6⁄
𝑣(6) = (1 − 𝑒 50) = 30.48437252
12.5

−𝑔𝑚 −𝑐∗𝑡⁄𝑚 𝑔𝑚 −𝑡 −𝑐∗𝑡⁄𝑚

𝑣′(𝑡) = (1 − 𝑒 )+ −𝑒
𝑐2 𝑐𝑚
−𝑔𝑚 −𝑐∗𝑡 𝑔𝑡 −𝑐∗𝑡
𝑣′(𝑡) = (1 − 𝑒 ⁄𝑚) + 𝑒 ⁄𝑚
𝑐2 𝑐

∆𝑣(𝑐̃) ≈ |𝑣′(𝑐̃) ∗ ∆𝑐̃|

−9.81 ∗ 50 −12.5∗6⁄ 9.81 ∗ 6 −12.5∗6⁄
∆𝑣(𝑐̃) ≈ 𝑣′(12.5) = |[ (1 − 𝑒 50) + 𝑒 50] ∗ 1.5 |
12.5 2 12.5
∆𝑣(𝑐̃) ≈ |−1.388074503 ∗ 1.5| = 2.082111755

𝑣 = 30.48437252 ± 2.082111755
8. Repita el problema anterior con g = 9.81, t = 6, c = 12.5 ± 1.5 y m = 50 ± 2.

𝜕𝑣 −𝑔𝑚 𝑔𝑡 −𝑐∗𝑡⁄𝑚
= 2 (1 − 𝑒−𝑐∗𝑡⁄𝑚) + 𝑒
𝜕𝑐 𝑐 𝑐

𝑔 𝑔𝑚 𝑐𝑡 −𝑐∗𝑡⁄
𝜕𝑣 = (1 − 𝑒−𝑐∗𝑡⁄𝑚) + −𝑒 𝑚
𝜕𝑚 𝑐 𝑐𝑚 2

𝜕𝑣 𝑔 𝑔𝑡
= (1 − 𝑒−𝑐∗𝑡⁄𝑚) − 𝑒−𝑐∗𝑡⁄𝑚
𝜕𝑚 𝑐 𝑚

𝜕𝑣 𝜕𝑣
∆𝑣(𝑐̃ , 𝑚
̃) ≈ | ∆𝑐̃| + | ∆𝑚̃|
𝜕𝑐 𝜕𝑚
𝑔 −𝑐̃∗𝑡⁄ 𝑔𝑡 −𝑐̃∗𝑡⁄
∆𝑣(𝑐̃, 𝑚
̃ ) ≈ 2.082111755 + |[ (1 − 𝑒 𝑚̃) − 𝑒 𝑚̃ ] ∗ ∆𝑚
̃|
𝑐 𝑚̃
9.81 −12.5∗6⁄ 9.81 ∗ 6 −12.5∗6⁄
∆𝑣(𝑐̃ , 𝑚
̃ ) ≈ 2.082111755+ |[ (1 − 𝑒 50 ) − 𝑒 50 ] ∗ 2|
12.5 50
∆𝑣(𝑐̃ , 𝑚
̃ ) ≈ 2.082111755 + |0.3470186258 ∗ 2|
∆𝑣(𝑐̃ , 𝑚
̃ ) ≈ 2.082111755 + 0.6940372516 = 2.776149007

𝑣 = 30.48437252 ± 2.776149007
9. Evalúe e interprete los números de condición para:

𝑓′(𝑥) = −𝑒−𝑥
𝑓(10) = 𝑒−10 = 4.54 ∗ 105
𝑓′(10) = −𝑒−10 = −4.54 ∗ 105
10(−4.54 ∗ 105)
𝑁ú𝑚𝑒𝑟𝑜 𝑑𝑒 𝑐𝑜𝑛𝑑𝑖𝑐𝑖ó𝑛 = −10
= 4.54 ∗ 105

𝑥𝑒𝑥 − 𝑒𝑥 + 1
𝑓′(𝑥) =
𝑥2
𝑒0.001 − 1
𝑓(10) = =1
0.001
0.001𝑒0.001 − 𝑒0.001 + 1
𝑓′(0.001) = = 0.5
0.0022
0.001(0.5)
𝑁ú𝑚𝑒𝑟𝑜 𝑑𝑒 𝑐𝑜𝑛𝑑𝑖𝑐𝑖ó𝑛 = = 5 ∗ 104
1
10. La fórmula de Manning para un canal rectangular se escribe como

donde Q = flujo (m3/s), n = coeficiente de rugosidad, B = ancho (m), H = profundidad (m) y

S = pendiente. Aplique la fórmula para un arroyo donde se conoce que el ancho = 20 m y la
profundidad =
0.3 m. Por desgracia conocemos el coeficiente de rugosidad y la pendiente con una
precisión de sólo ± 10%. Es decir, la rugosidad tiene un valor de 0.03 con un rango de
0.027 a 0.033, y la pendiente es 0.0003 con un rango de 0.00027 a 0.00033. Use un análisis
de error de primer orden para determinar la sensibilidad en la predicción del flujo para
cada uno de esos dos factores. ¿Cuál variable se debería intentar medir para una mejor
precisión?

5⁄
𝜕𝑄 = − 1 (𝐵𝐻) 3
𝜕𝑛 𝑛2 2 √𝑆
(𝐵 + 2𝐻) ⁄3
5
𝜕𝑄 1 (𝐵𝐻) ⁄3 1
= 2
𝜕𝑆 (𝐵 + 2𝐻)⁄𝑛3 2√𝑆

𝜕𝑄 𝜕𝑄
∆𝑄(𝑛̃ , 𝑆̃) = | ∆𝑛̃ | + | ∆𝑆̃ |
𝜕 𝜕𝑆
5⁄ 𝑛 5⁄
1 3 3
1 1
∆𝑄(𝑛̃ , 𝑆̃ ) = |[ 𝑛̃ 2
(𝐵𝐻) √𝑆̃ ] ∆𝑛̃ | + |[ 𝑛̃ (𝐵𝐻)
2⁄
3
2⁄
3 √ ̃] ∆𝑆̃|
(𝐵 + 2𝐻) (𝐵 + 2𝐻) 2𝑆
5⁄
1 3
̃ (20 ∗ 0.3)
∆𝑄(𝑛̃ , 𝑆) = |[ 20.03 2⁄ √0.0003 ] ∗ 0.0003|
3
(20 + 2 ∗ 0.3)
1 5⁄
+ |[ (20 ∗ 0.3) 3 1
0.03 2⁄ ] 0.00003|
(20 + 2 ∗ 0.3) 3 2√0.0003

∆𝑄 ≈ 0.152211211 + 0.076105605 = 0.228316816

11. Considere la función f(x) = x3 – 2x + 4 en el intervalo [–2, 2] con h = 0.25. Use las aproximaciones en diferencias finitas
hacia adelante, hacia atrás y centrada para la primera y segunda derivadas, e ilustre gráficamente qué aproximación es
más exacta. Grafique las tres aproximaciones a la primera derivada por diferencias finitas, junto con los valores exactos, y
haga lo mismo con la segunda derivada. Use Excel para ilustrar los resultados.

Primera derivada
f'(x)
intervalo 1ra derivada valores verdaderos f(x-1) f(x) f(x+1) hacia delante hacia atrás centrada
-2 10 -2.890625 0 2.140625 8.5625 11.5625 10.0625
-1.75 7.1875 0 2.140625 3.625 5.9375 8.5625 7.25
-1.5 4.75 2.140625 3.625 4.546875 3.6875 5.9375 4.8125
-1.25 2.6875 3.625 4.546875 5 1.8125 3.6875 2.75
-1 1 4.546875 5 5.078125 0.3125 1.8125 1.0625
-0.75 -0.3125 5 5.078125 4.875 -0.8125 0.3125 -0.25
-0.5 -1.25 5.078125 4.875 4.484375 -1.5625 -0.8125 -1.1875
-0.25 -1.8125 4.875 4.484375 4 -1.9375 -1.5625 -1.75
0 -2 4.484375 4 3.515625 -1.9375 -1.9375 -1.9375
0.25 -1.8125 4 3.515625 3.125 -1.5625 -1.9375 -1.75
0.5 -1.25 3.515625 3.125 2.921875 -0.8125 -1.5625 -1.1875
0.75 -0.3125 3.125 2.921875 3 0.3125 -0.8125 -0.25
1 1 2.921875 3 3.453125 1.8125 0.3125 1.0625
1.25 2.6875 3 3.453125 4.375 3.6875 1.8125 2.75
1.5 4.75 3.453125 4.375 5.859375 5.9375 3.6875 4.8125
1.75 7.1875 4.375 5.859375 8 8.5625 5.9375 7.25
2 10 5.859375 8 10.890625 11.5625 8.5625 10.0625
 Hacia delante
15

10

0
-3 -2 -1 0 1 2 3
-5

Valores verdaderos Hacia delante

Hacia atrás
15

10

0
-3 -2 -1 0 1 2 3
-5

Valores verdaderos Hacia atrás

Centrada
12
10
8
6
4
2
0
-3 -2 -1 -2 0 1 2 3
-4

Valores verdaderos Centrada

Segunda derivada

f'(x)
2da derivada hacia hacia
Intérvalo valores f(x-2) f(x-1) f(x) f(x+1) f(x+2) delante atrás centrada
verdaderos
-2 -12 -6.625 -2.890625 0 2.140625 3.625 -10.5 -13.5 -12
-1.75 -10.5 -2.890625 0 2.140625 3.625 4.546875 -9 -12 -10.5
-1.5 -9 0 2.140625 3.625 4.546875 5 -7.5 -10.5 -9
-1.25 -7.5 2.140625 3.625 4.546875 5 5.078125 -6 -9 -7.5
-1 -6 3.625 4.546875 5 5.078125 4.875 -4.5 -7.5 -6
-0.75 -4.5 4.546875 5 5.078125 4.875 4.484375 -3 -6 -4.5
-0.5 -3 5 5.078125 4.875 4.484375 4 -1.5 -4.5 -3
-0.25 -1.5 5.078125 4.875 4.484375 4 3.515625 0 -3 -1.5
0 0 4.875 4.484375 4 3.515625 3.125 1.5 -1.5 0
0.25 1.5 4.484375 4 3.515625 3.125 2.921875 3 0 1.5
0.5 3 4 3.515625 3.125 2.921875 3 4.5 1.5 3
0.75 4.5 3.515625 3.125 2.921875 3 3.453125 6 3 4.5
1 6 3.125 2.921875 3 3.453125 4.375 7.5 4.5 6
1.25 7.5 2.921875 3 3.453125 4.375 5.859375 9 6 7.5
1.5 9 3 3.453125 4.375 5.859375 8 10.5 7.5 9
1.75 10.5 3.453125 4.375 5.859375 8 10.890625 12 9 10.5
2 12 4.375 5.859375 8 10.890625 14.625 13.5 10.5 12
 Hacia delante
15
10
5
0
-3 -2 -1 0 1 2 3
-5
-10
-15

Valores verdaderos Hacia delantes

Hacia atrás
15
10
5
0
-3 -2 -1 0 1 2 3
-5
-10
-15

Valores verdaderos Hacia atrás

Centrada
15
10
5
0
-3 -2 -1 0 1 2 3
-5
-10
-15

Valores verdaderos Centrada