Solucionario-Tarea 3 - Series-Métodos Numéricos
Solucionario-Tarea 3 - Series-Métodos Numéricos
Solucionario-Tarea 3 - Series-Métodos Numéricos
DE TRUJILLO
FACULTAD DE INGENIERÍA
QUÍMICA
ESCUELA DE INGENIERÍA
AMBIENTAL
TAREA 3
EXPERIENCIA CURRICULAR:
MÉTODOS NUMÉRICOS
DOCENTE:
ALUMNO(S):
: • Pérez Castro Nayeli
• Pereda Zelada Daniel Arturo
• Polo Tolentino Anghely Mirella
• Rivera Salas Ana Laura
• Ramos Rosales Josh
• Reyes Dolores Exequiel
• Rodríguez Rodríguez Sofía Melisa
TRUJILLO-2021
RESOLUCIÓN DE PROBLEMAS
4.1. La serie infinita se utiliza para aproximar 𝒆𝒙
𝑥2 𝑥3 𝑥𝑛
𝑒𝑥 = 1 + 𝑥 + + + ⋯+
2! 3! 𝑛!
a) Muestre que la expansión en serie de Maclaurin es un caso especial de la expansión en la
serie de Taylor [ecuación (4.7)] con 𝒙𝒊 = 0 y h = x.
𝑒0 𝑒2 1
𝑓(𝑥),, = 𝑒 𝑥 n=0 𝑎0 = 0!
=1 n=2 𝑎2 = 2!
= 2!
𝑒1 𝑒3 1
𝑓(𝑥),,, = 𝑒 𝑥 n=1 𝑎1 = 1!
=1 n=3 𝑎3 = 3!
= 3!
𝑥2
• 𝑓(𝑥) = 𝑓(0) + 𝑓(0), 𝑥 + 𝑓(0),, +⋯
2!
• ,
𝑓(0) = 𝑓(0) = 𝑓(0) = 1; ,,
𝑥2 𝑥3 𝑥𝑛
• 𝑓(𝑥) = 1 + 𝑥 + + + ⋯+
2! 3! 𝑛!
b) Use la serie de Taylor para estimar f (x) = 𝒆−𝒙 en 𝒙𝒊+𝟏 = 1 para 𝒙𝒊 = 0.25. Emplee
versiones de cero, primero, segundo y tercer orden, y calcule |𝜺𝒕 | para cada caso.
ℎ2 ℎ3
𝑓( 𝑥𝑖+1 ) = 𝑒 −𝑥𝑖 − 𝑒 −𝑥𝑖 ℎ + 𝑒 −𝑥𝑖 + 𝑒 −𝑥𝑖 +⋯
2! 3!
𝜀 0,367879−0,778801 = −111,7%
𝑡= ∗100
0,367879
ORDEN UNO
𝑓(1) ≅ 0,778801 − 0,778801(0,75) ≅ 0,1947
𝜀 0,367879−0,1947
𝑡= ∗100 ≅47,1%
0,367879
SEGUNDO ORDEN
(0.75)2
𝑓(1) ≅ 0,778801 − 0,778801(0,75) + 0,778801 ≅ 0,413738
2
𝜀 0,367879−0,413738
𝑡= ∗100 ≅−12,46%
0,367879
TERCER ODEN
(0.75)3
𝑓(1) ≅ 0,413738 − 0,778801 = 0,358978
6
𝜀 0,367879−0,358978
𝑡= ∗100 ≅2,42%
0,367879
𝜀 0,707107−1
𝑡= ∗100 ≅−41,42%
0,707107
PRIMER ORDEN
𝜋 2
𝜋 (4 )
cos ( ) = 1 − = 0,691575
4 2
𝜀 0,707107−0,691575
𝑡= ∗100 ≅2,19%
0,707107
𝜀 0,691575−1
𝑎= ∗100 ≅−44,6%
0,691575
SEGUNDO ORDEN
𝜋 4
𝜋 ( )
cos ( ) = 0,691575 + 4 = 0,707429
4 24
𝜀 0,707107−0,707429
𝑡= ∗100 ≅−0.45%
0,707107
𝜀 0,707429−0,691575
𝑎= ∗100 ≅2,24%
0,707429
TERCER ORDEN
𝜋 4
𝜋 (4 )
cos ( ) = 0,707429 − = 0,707130
4 120
𝜀 0,707107−0,707130
𝑡= ∗100 ≅0,0005%
0,707107
𝜀 0,707130−0,707429
𝑎= ∗100 ≅ −0,046 %
0,707130
4.3. Repita los cálculos del problema 4.2, pero ahora usando la expansión de la
serie de Maclaurin para sen x, para evaluar el sen (𝝅/3).
𝑥3 𝑥5 𝑥7
𝑆𝑒𝑛𝑥 = 𝑥 − + − +⋯
3! 5! 7!
𝜀𝑠 = 0,5 𝑥 102−2 = 0,5%
ORDEN CERO
sen (π/3) =1
Valor verdadero=0.8660
𝜀 0,8660−1
𝑡= ∗100 ≅−0,154%
0,8660
PRIMER ORDEN
𝜋 3
𝜋 (3 )
Sen ( ) = 1 − = 0.8086
3 6
𝜀 0.8660−0.8086
𝑡= ∗100 ≅6.62%
0,8660
𝜀 0.8086−1
𝑎= ∗100 ≅−23.6%
0.8086
SEGUNDO ORDEN
𝜋 5
𝜋 (3 )
Sen ( ) = 0.8086 + = 0.8190
3 120
𝜀 0.8660−0.8190
𝑡= ∗100 ≅5.42%
0,8660
𝜀 0,8190−0.8086
𝑎= ∗100 ≅1.26%
0.8190
4.5. Emplee la expansión de la serie de Taylor de cero hasta tercer orden para
predecir f (3) si:
𝑓(𝑥) = 25𝑥 3 − 6𝑥 2 + 7𝑥 − 88
usando como punto base x = 1. Calcule el error relativo porcentual verdadero para cada
aproximación.
Donde: h =x-a
✓ ORDEN CERO: FUNCIÓN POLINOMIAL
Recuerda
𝑓(𝑥) = 25𝑥 3 − 6𝑥 2 + 7𝑥 − 88 = 102 h =x-a
h=3-1
✓ PRIMER ORDEN h=2
𝑓(2), = 75𝑥 2 − 12𝑥 + 7=283
✓ SEGUNDO ORDEN
𝑓(2),, = 150𝑥 − 12 = 288
✓ TERCER ODEN
𝑓(2),,, = 150
1
Reemplazando en la ecuación
𝑓(2), 𝑓(2), 𝑓(3),
F (3) =f (1) + (2) + (4) + (8)
1! 2! 3!
F (3) =102+566+576+200
F (3) =1444 ← 𝑉𝑎𝑙𝑜𝑟 𝑣𝑒𝑟𝑑𝑎𝑑𝑒𝑟𝑜
ERRORES REALES
• 1 er término: 𝑓(3) = 102
𝜀 1444−102
𝑟= ∗100 ≅92,93%
1444
• 2do término
F (3) =102+566
F (3) =385
𝜀 1444−668
𝑟= ∗100 ≅53,74%
1444
• 3er término
F (3) =102+566+576
F(3)=1244
𝜀 1444−1244
𝑟= ∗100 ≅13.85%
1444
• 4to término
F (3) =102+566+576+200
F(3)=1444
𝜀𝑟=1444−1444∗100 ≅0%
1444
4.6. 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.
𝑓(2,5) = 0,916
• Orden 0:
𝑓(2,5) ≅ 𝑓(1) ≅ 0
• Orden 1:
1
𝑓 ′ (1) = = 1
1
• Orden 2:
1
𝑓 ′′ (1) = − 12
= −1
• Orden 3:
2
𝑓 ′′′ (1) = 13
=2
• Orden 4:
6
𝑓 ′′′′ (1) = − 14
= −6
102 − 39,859
• 𝑓 ′ (2) = 0,25
= 248,56
182,14 − 39,859
• 𝑓 ′ (2) = 2(0,25)
= 284,5
𝑓 ′′ (𝑥𝔦 )ℎ
ℰ𝑡 ≈
2
𝑓 ′′ (2) = 150(2) − 12 = 288
288(0,25)
ℰ𝑡 ≈ = 36
2
𝑓 ′′′ (𝑥𝔦 )ℎ2
ℰ𝑡 ≈ −
2
150(0,25)2
ℰ𝑡 ≈ − = −1,56
6
ℰ𝑡 = 283 − 284,56
• Error exacto:
𝐻(670) − 𝐻(630)
𝛥𝐻 = = 168,1693
2
𝑇 = 650 ± 40
𝛥𝐻 = 8,4084(40) = 3366,3387
4.10. Repita el problema 4.9, pero para una esfera de cobre con 𝒓𝒂𝒅𝒊𝒐 = 𝟎. 𝟏𝟓 ±
𝟎. 𝟎𝟏 𝒎, 𝒆 = 𝟎. 𝟗𝟎 ± 𝟎. 𝟎𝟓 𝒚 𝑻 = 𝟓𝟓𝟎 ± 𝟐𝟎.
𝐴 = 4𝜋𝑟 2
𝑑𝐴
𝛥𝐴(𝑟) = ( ) 𝛥(𝑟)
𝑑𝑟
𝑑𝐴
= 8𝜋𝑟1 = 8 × 3.1415 × 0.15 = 3.7698
𝑑𝑟
ΔA = 3.7698(0.01) = 0.037698
ⅆH
= 4AeσT 3
ⅆT
ⅆH
= 4(0.037698)(0,90)(5,67 × 10−8 )(550)3
ⅆT
H(570) − H(530)
ΔH = = 25.6048
2
T = 550 ± 20
ΔH = 1.28024(20) = 25.6048
4.11. Recuerde que la velocidad de caída del paracaidista puede calcularse con
[Ec. (1.10)]
𝒈𝒎
𝒗(𝒕) = (𝟏 − 𝒆−𝒄×𝒕⁄𝒎 )
𝒄
Use un análisis de error de primer orden para estimar el error de 𝒗 𝒑𝒂𝒓𝒂 𝒕 =
𝟔, 𝒔𝒊 𝒈 = 𝟗. 𝟖𝟏 𝒚 𝒎 = 𝟓𝟎, 𝒑𝒆𝒓𝒐 𝒄 = 𝟏𝟐. 𝟓 ± 𝟏. 𝟓.
𝑑𝑣 𝑐 × 𝑔 × 𝑡 × 𝑒 −𝑐×𝑡⁄𝑚 − 𝑔 × 𝑚(1 − 𝑒 −𝑐×𝑡⁄𝑚 )
=
𝑑𝑐 𝑐2
𝑑𝑣 12.5 × 9.81 × 6 × 𝑒 −12.5×6⁄50 − 9.81 × 50(1 − 𝑒 −12.5×6⁄50 )
=
𝑑𝑐 12.52
𝑑𝑣 735.75 × 𝑒 −1.5 − 490.5(1 − 𝑒 −1.5 ) 164.19354 − 381.03763
= = = −1.38780
𝑑𝑐 156.25 156.25
𝑑𝑣
𝛥𝑣(𝑐) = ( ) 𝛥(𝑐)
𝑑𝑐
𝛥𝑣(𝑐) = 1.38780(2)
𝛥𝑣(𝑐) = 2.77560
490.5
𝑣(12.5) = 9.81(50)(1 − 𝑒 −12.5×6⁄50 ) = (1 − 𝑒 −1.5 ) = 30.48437
12.5
𝑣 = 30.48437 ± 2.77560
𝑑𝑣 𝑔×𝑡 𝑔
| |= × 𝑒 −𝑐⁄𝑚𝑡 + (1 − 𝑒 −𝑐⁄𝑚𝑡 ) = 0.871467
𝑑𝑚 𝑚 𝑐
𝑑𝑣
| | 𝛥𝑚 = 0.871467(1.5) = 1.3072
𝑑𝑚
𝑣 = 2.77560 + 1.3072 = 4.0805
𝑣 = 30.48437 ± 4.0805
𝑣 = 30.4533 ± 0.692
SOLUCIÓN:
a) 1
1.0001 [ ]
2√1.0001 − 1
𝑁𝐶 = = 50.00
√1.0001 − 1 + 1
✓ Es una mala condición porque f’(1) es grande cerca a x=1
0 1 2 x
b) 9(−𝑒 −9 )
𝑁𝐶 = = −9
𝑒 −9
✓ Es una mala condición porque ‘‘x’’ es grande.
c) 200
(200) −1
√2002 + 1
𝑁𝐶 =
√2002 + 1 − 200
200(−1.2𝑥10−5 )
𝑁𝐶 = = −0.96
0.0025
✓ Es una buena condición.
d)
𝑥(−𝑥𝑒 −𝑥 − 𝑒 −𝑥 + 1)
𝑁𝐶 = 𝑥2
−𝑥
𝑒 −1
( 𝑥 )
0.01(−0.1(0.99) − (0.99) + 1
0.012 0.9
𝑁𝐶 = = = −0.9
0.99 − 1 −1
0.01
✓ Es una buena condición.
𝑠𝑒𝑛𝑥
𝑓(𝑥) =
1 + 𝑐𝑜𝑠𝑥
(1 + 𝑐𝑜𝑠𝑥)(cos (𝑥)) + 𝑠𝑖𝑛𝑥(𝑠𝑖𝑛𝑥)
𝑓′(𝑥) =
(1 + 𝑐𝑜𝑠𝑥)2
1.001𝜋 𝑓′(1.001𝜋)
𝑁𝐶 =
𝑓(1.001𝜋)
3.144(202642)
𝑁𝐶 = = −1001
−636.6
✓ Es una mala condición porque
1 + cos(1.001𝜋) 𝑥0
4.14. Empleando las ideas de la sección 4.2, muestre las relaciones de la tabla 4.3.
SOLUCIÓN:
𝑓(𝑢, 𝑣) = 𝑢 + 𝑣
𝑑𝑓 𝑑𝑓
∆𝑓 = | | ∆𝑢 + | | ∆𝑣̃
𝑑𝑢 𝑑𝑣
𝑑𝑓 𝑑𝑓
| |=1 | |=1
𝑑𝑢 𝑑𝑣
∆𝑓(𝑢̃, 𝑣̃) = ∆𝑢̃ + ∆𝑣̃
• Multiplicamos • Dividimos
𝑓(𝑢, 𝑣) = 𝑢. 𝑣 𝑢
𝑓(𝑢, 𝑣) =
𝑣
𝑑𝑓 𝑑𝑓
| |=𝑣 | |=𝑢 𝑑𝑓 1 𝑑𝑓 𝑢
𝑑𝑢 𝑑𝑣 = =− 2
𝑑𝑢 𝑣 𝑑𝑣 𝑣
̃ ∆𝑢̃ + |𝑢̃|∆𝑣̃
∆𝑓(𝑢̃, 𝑣̃) = |𝑣|
1 𝑢 |𝑣|∆𝑢+|𝑢|∆𝑣
∆𝑓(𝑢, 𝑣) = |𝑣| ∆𝑢 + |𝑣2 | ∆𝑣 = |𝑣 2 |
4.15 Muestre que la ecuación (4.4) es exacta para todos los valores de x, si f(x) = ax2 + bx + c
SOLUCIÓN:
𝑓(𝑥) = 𝑎𝑥 2 + 𝑏𝑥 + 𝑐
𝑓 ′ (𝑥) = 2𝑎𝑥 + 𝑏
𝑓 ′′ (𝑥) = 2𝑎
2𝑎
𝑎𝑥𝑖+1 2 + 𝑏𝑥𝑖 + 𝑐 + 2𝑎𝑥𝑖 + 𝑏(𝑥𝑖+1 − 𝑥𝑖 ) + (𝑥 2 − (2𝑥𝑖+1 )𝑥𝑖+ 𝑥𝑖 2 )
𝑎 𝑖+1
𝑎(𝑥𝑖 2 + 2𝑥𝑖 𝑥𝑖+1 − 2𝑥𝑖 2 + 𝑥𝑖+1 2 − 2𝑥𝑖+1 (𝑥𝑖 ) + 𝑥𝑖 2 )
𝑎𝑥𝑖+1 2
donde 𝑄 = flujo (𝑚3 /s), 𝑛 = coeficiente de rugosidad, 𝐵 = ancho (m), 𝐻 = 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?
SOLUCIÓN
Datos
Q= Flujo (m3/s)
N= Coeficiente de rugosidad=0.03
B = Ancho (m)=20 m
H= Profundidad (m) =0.3m
S= Pendiente=0.0003
A partir los datos encontramos el flujo
5
1 (20 ∗ 0.3)3
Q= 2 √0.0003
0.03
(20 + 2 ∗ 0.3)3
5
1 (6)3
Q= 2 (0.017)
0.03
(20.6)3
19.81 ∗ 0.017
Q= = 1.50(m3/s)
0.225
El error estimado usando números inexacto
• ∆𝑛 =0.03-0.027=0.003
∆𝑛 =0.03-0.033=-0.003 ∆𝒏 =±𝟎. 𝟎𝟎𝟑
• ∆𝑠 =0.0003 − 0.00027 = 0.00003
∆𝒔 =±𝟎. 𝟎𝟎𝟎𝟎𝟑
∆𝑠 =0.0003 − 0.00033 = −0.00003
∞
𝑓 (𝑛) (0)
𝑓(𝑥) = ∑ (𝑥 − 0)𝑛
𝑛!
𝑛=0
1
Derivamos la función 1−𝑥
1
• 𝑓(𝑥), = (1−𝑥)2
2
• 𝑓(𝑥),, = (1−𝑥)3
6
• 𝑓(𝑥),, = (1−𝑥)4
• …
• …
𝑛!
• 𝑓(𝑥)𝑛 = (1−𝑥)𝑛+1
Sustituyamos x =0
✓ 𝑓(𝑥), = 1
✓ 𝑓(0),, = 2
✓ 𝑓(0),, = 6
✓ …
✓ 𝑓 𝑛 (x) = 𝑛!
1 𝑓(0),, 𝑥 2 𝑓(0),,, 𝑥 3
= 𝑓(0) + 𝑓(0), + + + ⋯,
1−𝑥 2! 3!
= 𝟏 + 𝒙 + 𝒙𝟐 + 𝒙𝟑 + 𝒙𝟒 …
4.19. Para calcular las coordenadas espaciales de un planeta tenemos que resolver
la función
𝒇(𝒙) = 𝒙 − 𝟏 − 𝟎. 𝟓𝒔𝒆𝒏𝒙
Sea 𝒂 = 𝒙 = 𝝅⁄𝟐 en el intervalo [𝟎, 𝝅]el punto base. Determine la expansión de la
serie de Taylor de orden superior que da un error máximo de 𝟎, 𝟎𝟎𝟏𝟓 en el intervalo
dado. El error es igual al valor absoluto de la diferencia entre la función dada y la
expansión de la serie de Taylor especificada. (Sugerencia: Resuelva gráficamente).
3 3𝜀 3𝜀 1/3
ℎó𝑝𝑡 = √ = ( ) = (3𝜀. 𝑀−1 )1/3
𝑀 𝑀
𝒅 𝒏 𝒉′ó𝒑𝒕 = (𝟑𝜺. 𝑀 −1 )𝟏/𝟑
𝑽 = 𝒏𝑽𝒏−𝟏 . 𝑽′
𝒅𝒙
1 2
ℎ′ó𝑝𝑡 = (3𝜀. 𝑀−1 )−3 . (3𝜀. 𝑀−1 )′
3
1 2
ℎ′ó𝑝𝑡 = (3𝜀. 𝑀−1 )−3 . (−3𝜀. 𝑀−2 )
3
1 1 −2
ℎ′ó𝑝𝑡 = . 2 . (−3𝜀. 𝑀 )
3 −1
(3𝜀. 𝑀 )3
1 1
ℎ′ó𝑝𝑡 = . 3 . (−3𝜀. 𝑀−2 )
3 √(3𝜀. 𝑀−1 )2
1 1 3𝜀
ℎ′ó𝑝𝑡 = . . (− 2 )
3 3 3𝜀 2 𝑀
√( )
𝑀
3𝜀
ℎ′ó𝑝𝑡 = −
3
3𝜀 2 2
3. √( ) .𝑀
𝑀
𝜀
ℎ′ó𝑝𝑡 = −
2
√(3𝜀 ) . 𝑀2
3