4.

1 La serie infinita
2 3 n
x x x
e x =1+ x + + +…+
2 3! n!
se utiliza para aproximar ex.
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 xi = 0 y h = x.
x2
f ( x )=f ( 0 ) + f ' ( 0 ) + f ' ' ( 0 ) + …
2
' '' 0
f ( 0 )=f ( 0 ) =f ( 0 )=e =1
2
x
f ( x )=1+ x + + …
2

b) Use la serie de Taylor para estimar f(x) = e–x en xi+1 = 1 para xi = 0.25. Emplee
versiones de cero, primero, segundo y tercer orden, y calcule |εt| para cada caso.

−x i −x i −x i h2 −x h3
f ( x i+1 ) =e −e h+ e −e i
+…
2 2

Para xi = 0.2, xi+1 = 1 y h = 0.8. verdad = e–1 = 0.367879.

orden cero:

f ( 1 )=e−0.2=0.818731

εt =|0.367879−0.818731
0.367879 |∗100 %=122.55 %
Primer orden:
f (1) =0.818731 -0.818731(0.8) =0.163746

εt =|0.367879−0.163746731
0.367879 |∗100 %=55.49 %
Segundo orden:
0.82
f ( 1 )=0.818731−0.818731 ( 0.8 ) +0.818731 A =0.42574
2

εt =|0.367879−0.42574
0.367879 |∗100 %=15.73 %
Tercer orden:
2 6
0.8 0.8
f ( 1 )=0.818731−0.818731 ( 0.8 ) +0.818731 A −0.818731 A =0.355875
2 6
 |0.367879−0.355875
εt =
0.367879 |∗100 %=3.26 %
4.2 La expansión en serie de Maclaurin para cos x es

x2 x4 x6 x8
e x =1− + + + …
2! 4! 6! 8!
Iniciando con el primer término cos x = 1, agregue los términos uno a uno para
estimar cos (p/4). Después de que agregue cada uno de los términos, calcule los
errores relativos porcentuales exactos 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 dos cifras
significativas.

ε s=0.5∗102−2 %=0.5 %

Valor verdadero: cos (pi/3) = 0.5

Orden cero:
π
cos()
3
=1

|0.5−1
ε s=
0.5 |
∗100 % 100 %

Primer orden:

π 2
π 3 ()
cos()
3
=1−
2
=0.451689

ε s=9.66 %

|0.451689−1
ε s=
0.451689 |
∗100 %=121.4 %

Segundo orden:
4
π
π 3 ()
cos()
3
=0.451689+
24
=0.501796

ε t =0.359 %
 |0.501796−0.451689
εt =
0.501796 |∗100 %=9.986 %

Tercer orden:
6
π
π 3 ()
cos()
3
=0.501796−
720
=0.499965

ε t =0.00709 %

ε a= |0.499965−0.501796
0.499965
89
|∗100 %=0.366 %
Dado que el error aproximado es inferior al 0,5%, el cálculo puede finalizarse.

4.3 Repita los cálculos del problema 4.2, pero ahora usando la
expansión de la serie de Maclaurin para sen x,

x3 x5 x7
sen x =x− + − + …
3 ! 5! 7 !
para evaluar el sen (p/4).

Use el criterio de parada:

2−2
ε s=0.5∗10 %=0.5 %

Valor verdadero: sin ( π /3) = 0.866025…

Orden cero:
π π
sen ()
3
= =1.047198
3

|0.866025−1.047198
εt =
0.866025 |∗100 %=20.92 %
Primer orden:
 3
π
π 3 ()
sin()
3
=1.047198−
6
=0.855801

ε t =1.18 %
0.855801−1.047198
ε a= | 0.855801 |
∗100 %=22.36 %

Segundo orden:

π 5
π 3 ()
sen ()
3
=0.855801+
120
=0.866295

ε t =0.031 %

εt =|0.866295−0.855801
0.866295 |∗100 %=1.211 %

Tercer orden:
7
π
π 3 ()
sin()
3
=0.866295−
5040
=0.866021

ε t =0.000477 %

ε a= |0.866021−0.866295
0.866021 |∗100 %=0.0316 %
Dado que el error aproximado es inferior al 0,5%, el cálculo puede finalizarse.

4.4 Emplee la expansión de la serie de Taylor de cero hasta tercer orden para
predecir f (2) sí.
3 2 2
senx=25 x −65 x +7 x−88=π r

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

para cada aproximación.

f ( 3 ) =554
Orden cero:

554−(−62)
εt =| 554 |
∗100 %=111.191 %

Primer orden:
f ( 3 ) =−62+f ' ( 1 ) (3−1 )=−62+ 70 ( 2 )=78

ε t =85.921 %

Segundo orden:

''
f ( 1)
f ( 3 ) =78+ ( 3−1 )3=78+ 138 4=354
2 2

ε t =36.101 %

Tercer orden:

3
f ( 1)
f ( 3 ) =78+ ( 3−1 )3 =354+ 150 8=554
6 6

ε t =0 %

Por lo tanto, el resultado de tercer orden es perfecto porque la función original es

un tercer orden
polinomio.

4.5 Use la expansión de la serie de Taylor de cero al cuarto orden para estimar f
(3) si f(x) = ln x utilizando x = 1 como punto base.
Calcule el error relativo porcentual et para cada aproximación.
Analice los resultados.
Valor verdadero:

f ( 2.5 ) =ln ( 2.5 )=0.916291.. .

orden cero:

f ( 2.5 ) =f ( 1 )=0

εt =|0.916291−0
0.916291 |
∗100 %=100 %
Primer orden:
f ( 2.5 ) =f ( 1 ) ( 2.5−1 )=0+1 (1.5 )=1.5

εt =|0.916291−1.5
0.9162911.5 |
∗100 %=63.704 %

Segundo orden:
2.5 1
¿
¿
5.1
¿
¿
(¿2¿)=0.375
−1
(¿ 2¿)=1.5+ ¿
2
f '' (1 )
f ( 2.5 ) =1.5− ¿
2

εt =|0.916291−0.375
0.916291 |∗100 %=59.074 %
Tercer orden:
f (3) (3) 2 3
f ( 2.5 ) =0.375+ ( 2.51) =0.375+ 1.5 =1.5
6 6

εt =|0.916291−1.5
0.916291 |
∗100 %=63.704 %

Cuarto orden:
(1 )
f (1) −6
f ( 2.5 ) =1.5+ (2.5 1)(4) =1.5+ 1.5 4=0.234375
24 24

εt =|0.916291−0.234375
0.916291 |∗100 %=74.421 %
Por lo tanto, el proceso parece ser divergente, lo que sugiere que se requeriría un
paso más pequeño para convergencia.

4.6 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 4.4. Evalúe la derivada en x = 2 usando
un tamaño del incremento 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.
 ' 2
f ( x )=75 x −12 x +7

f ' ( 2 )=75( 2)2−12 ( 2 ) +7=283

valores de función:

xi–1 = 1.8 f(xi–1) = 50.96

xi = 2 f(xi) = 102
xi+1 = 2.2 f(xi+1) = 164.56

adelante:

' 164.56−102
f (2)= =312.8
0.2
283−312.8
εt = | 283 |
∗100 %=10.53 %

hacia atrás:

' 102−50.96
f (2)= =255.2
0.2

εt = |283283255.2|∗100 %=9.823 %
centrado:

' 164.56−50.96
f (2)= =284
2( 0.2)

εt = |283−284
283 |
∗100 %=0.353%

Tanto el avance como el retroceso tienen errores que pueden aproximarse

mediante (recuperar la ecuación 4.15),
4.7 Con la aproximación en diferencias centrales de O(h2) estime la segunda
derivada de la función examinada en el problema 4.4. Realice la evaluación para x
= 2 usando un tamaño de incremento 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.

f ' ' ( x ) =150 x−12

f ' ' (2 )=150 ( 2 )−12=288

h = 0.25:

f ( 2.25 )−2 f ( 2 ) + f ( 1.75 ) 182.1406−2 (102 ) +39.85938

f ' ' (2 )= 2
= 2
=288
0.25 0.25

h = 0.125:

f ( 2.125 )−2 f ( 2 ) + f ( 1.875 ) 139.6738−2 ( 102 ) +68.82617

f ' ' (2 )= = =288
0.1252 0.1252

4.8 Recuerde que la velocidad de caída del paracaidista puede calcularse con
[ecuación (1.10)]

gm
v ( t )= ( 1−e(c/ m)t )
c

Use un análisis de error de primer orden para estimar el error de v para t = 6, si g

= 9.8 y m = 50, pero c = 12.5 ± 2.
− ( c/ m ) t
∂ v cgt −gm ( 1−e−(c /m )t )
= =−1.38666
∂c c2

∆ v (c )= |∂∂ vc|∆ c=1.38666 ( 1.5)=2.079989

9.8( 50)
v ( 12.5 )= ( 1−e−12.5 (6 )/ 50) =30.4533
12.5

v-30.4533 =2.079989

Por lo tanto, los límites calculados con el análisis de primer orden van desde
28.3733 a 32.5333. Esta
El resultado se puede verificar calculando los valores exactos como
9.8(50)
v ( c−∆ c )= ( 1−e−(11/50)/ 6 )=32.6458
11

9.8 (50)
v ( c +∆ c ) = ( 1−e−(14/ 50)/ 6 )=28.4769
11
Por lo tanto, el rango de 2.0844 es cercano a la estimación de primer orden.

4.9 Repita el problema 4.8 con g = 9.8, t = 6, c = 12.5 ± 1.5 y m = 50 ± 2.

 9.8( 50)
v ( 12.5 )= ( 1−e−12.5 (6 )/ 50) =30.4533
12.5

∆ v (c )=|∂∂ vc|∆ c+|∂∂ vc|∆ m

− ( c/ m ) t
∂ v cgt −gm ( 1−e−(c /m )t )
= =−1.38666
∂c c2

∂ v ¿ − (c/ m )t g
= e + ( 1−e−( c/ m) t ) =0.871467
∂c m c

v =30.4533 = 3.822923

4.10 La ley de Stefan-Boltzmann se utiliza para estimar la velocidad de cambio de

la energía H para una superficie, esto es, H = AeσT 4
donde H está en watts, A = área de la superficie (m2), e = la emisividad que
caracteriza la propiedad de emisión de la superficie (adimensional), σ = una
constante universal llamada constante de Stefan-Boltzmann (= 5.67 × 10–8 W m–
2 K–4) y T = temperatura absoluta (K). Determine el error de H para una placa de
acero con A = 0.15 m2, e = 0.90 y T = 650 ± 20. Compare los resultados con el
error exacto. Repita los cálculos, pero con T = 650 ± 40. Interprete los resultados.

Para ∆ (T )=20,

∆ H (T )= |∂∂TH |∆ T
|∂∂TH |=4 Aeσ T =4 ( 0.15 ) 0.9 (5.67∗10
3 −8
) 6503 =8.408

∆ H ( T )=8.408 ( 20 )=168.169

Error exacto:

H ( 670 ) −H ( 630 )
∆ H Verdad = =168.3286
2
Por lo tanto, la aproximación de primer orden está cerca del resultado exacto.

Para ∆ T =40

∆ H ( T )=8.408 ( 40 ) =336.3387
Error exacto:

H ( 690 ) −H ( 610 ) 1735.055−1059.83

∆ H Verdad = = =337.6124
2 2

Nuevamente, la aproximación de primer orden está cerca del resultado exacto.

Los resultados son buenos porque H (T) es casi lineal en los rangos que estamos
examinando. Esto se ilustra con el siguiente
trama.

4.11 Repita el problema 4.10, pero para una esfera de cobre con radio = 0.15 ±
0.01 m, e = 0.90 ± 0.05 y T = 550 ± 20. Para una esfera, A = 4r2. Por lo tanto,

H=4 π r 2 eσ T 4
En los valores medios de los parámetros,

H (0.15,0.9,550) = 4 π (0.15) +0.90(5.67 -10) (550) =1320.288

∆ H= |∂∂Hr |∆ r +|∂∂He |∆ e+|∂∂tH |∆ T

∂H 2 4
=8 πe r σ T =17,603.84
∂r
∂H
=4 πe r 2 σ T 4=1466.987
∂e
∂H
=16 πe r 2 σ T 3=9.6021
∂T
H 17603.84(0.01) 1466.987(0.05) 9.6021(20) 441.4297
Para verificar este resultado, podemos calcular

h ( 0.14,0 .85,530 )=4 ( 0.14 ) 0.85 ( 5.67∗10 ) (530 )=936.6372

H0.16,0.95,570) =4 (0.16) 0.95(5.67 *10) (570) =1829.178
1829.178−936.6372
∆ H Verdad = =446.2703
2
4.12 Evalúe e interprete los números de condición para

xf ' ( x)
CN=
f (x)

(a)CN =
1.00001
[
2 √ 1.000011
=
]
1.00001(158.1139)
=157.617
√1.000011+1 1.003162
El resultado está mal acondicionado porque la derivada es grande cerca de x = 1.
10 (−e−10) 10 (−4.54∗10−5 )
(a)CN = −10
= −5
=−10
e 4.54∗10
El resultado está mal acondicionado porque x es grande.

300

( C ) CN =
300
[ √3002 +1
+1
] =
300 (−5.555556∗10−6 )
=−0.99999444
√ 3002 +1−300 0.0016667

El resultado está bien acondicionado.

−x e−x −e−x +1
x
x2 0.001(0.499667)
(a)CN = −x
= =−0.0005
e −1 −0.9995
x ( )
El resultado está bien acondicionado.

( 1cos ) cosx +sinx (sin)

x
( 1+cosx )2 3.141907(20,264,237)
(e)CN = = =−10,001
senx −6366.2
1+ cosx ( )
El resultado está mal acondicionado porque, como en la siguiente gráfica, la
función tiene una singularidad en
x=
4.13 Empleando las ideas de la sección 4.2, muestre las relaciones de la tabla 4.3.

f ( u , v ) =u+v

∆ H= |∂∂uf |∆ u+|∂∂ fv|∆ v

|∂∂uf |=1
|∂∂ fv|=1
f (u~, v~) u~ v~
Multiplicación:
f ( u , v ) =u . v

|∂∂uf |=v
|∂∂ fv|=u
f (u~, v~) u~ v~
División:
f ( u , v ) =u/v

|∂∂uf |= 1v
|∂∂ fv|= vu 2

1 1
f (u~, v~)  ||
v
u~  ||
v2
v~
 |v| ∆ u +|v| ∆ v
f (u~, v~) =
|v 2|
4.14 Muestre que la ecuación (4.4) es exacta para todos los valores de x, si f(x) =
ax2 + bx + c.

f ( x )=ax2 +bx +c

f ' ( x )=2 ax 2 +b
f ' ' ( x )=2 a
Sustituir estas relaciones en Eq. (4.4),

2a 2
ax 2i+1 +bx 2i+1 +c=ax2i +1+ bx2i +1+ c+ ( 2a x i +b ) ( x i+ 1+ x i ) + ( x i+1 −2 a x i+ 1 x i+ x2i )
2!

Recoger términos
2a 2
ax 2i+1 +bx 2i+1 +c=ax2i +1+ bx2i +1+ c+2 a x i ( x i +1+ x i ) + a ( x 2i +1−2 a xi +1 x i + x 2i) + ( x i+1−2 a x i+1 xi + x 2i )
2!

ax 2i+1 +bx 2i+1 +c=ax2i +1+ bx2i +1+ c+2 a x i +b ( x i+1 + x i )−2 a x i+1 x i + x 2i+ b xi +1 x i + x 2i

ax 2i+1 +bx 2i+1 +c=( ax 2i+1 +2 a xi +bx 2i+1 ) + a xi + ( 2 a x i−2 a x i )( b xi +1+ b xi ) + x 2i+1 −2 a x i+ 1 x i+ x2i
2 2 2 2
ax i+1 +bx i+1 +c=ax i +1+ bx i +1+ c

4.15 La fórmula de Manning para un canal rectangular se escribe

como
5
3
1 ( BH )
Q= 2
S 1/ 2
n
( B+2 H )3

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 se debería intentar medir para una mejor precisión?
 ∆ Q= |∂Q∂n |∆ u+|∂∂QS |∆ S
5 /3
∂Q −1 ( BH )
= 2 S 0.5=−50.47
∂ n n ( B+ 2 H )2 /3
5 /3
∂Q −1 ( BH ) 1
= 2 2 /3 0.5
=2536 .9
∂ S n ( B+ 2 H ) S

Q 50.74]0.0032536.9]0.000030.1520.076 0.228

4.16 Si |x| < 1, se sabe que

ε s=0.5∗10 2−2
% =0.5 %

Valor verdadero:
1/ (1-0.1) =1.111111…….
orden cero:
1
=1
1−x

ε a= |1.11111−1
1.11111 |
∗100 %=10 %

Primer orden:
1
=1+ 0.1=1.1
1−x
ε t =1 %

ε a= |1.1−1
1.1 |
∗100 %=9.0909 %

segundo orden:
1
=1+ 0.1+ 0.01=1.11
1−x
ε t =0.1 %

ε a= |1.11−1.1
1.11 |
∗100 %=0.9009009 %
tercer orden:
1
=1+ 0.1+ 0.01=1.111
1−x
ε t =0. 0 1 %

ε a=|1.11−1.1
1.11 |
∗100 %=0.090009 %

4.17 Un misil sale de la Tierra con una velocidad inicial v0 formando con la vertical
un ángulo φ0 como se muestra en la figura

P4.17. La altitud máxima deseada es a R donde R es el radio de la Tierra. Usando

las leyes de la mecánica se demuestra que

2
α ve
sen ∅ 0=( 1+α ) 1−
√
1+α v o ( )
donde ve es la velocidad de escape del misil. Se quiere disparar el misil y alcanzar
la velocidad máxima proyectada con una exactitud de ±1%. Determine el rango de
valores de f0 si ve/v0 = 2 y
a = 0.2.
dsen ∅0
∆ ( sen ∅0 ) = | dα | ∆α

10.25−0.25(4)
dsen ∅0 β 0.25 ( 4 )
dα
=
2 √ ( 1+0.25 ) ( ¿ ) √
+ 1−
1+0.25
=−3.1305

donde = (ve/v0)2 = 4 y = 0.25

 10.25−0.25(4)
dsen ∅0 β 0.25 ( 4 )
dα
=
√
2 √ ( 1+0.25 ) ( ¿ )
+ 1−
1+0.25
=−3.1305

∆ ( sen ∅0 ) =3.1305 ∆ α

para ∆ α =0.25( 0.02)=0.005,

∆ ( sen ∅0 ) =3.1305 ( 0.005 )=0.015652

0.25
1− 4=¿ 0.559017
10.25
sen ∅0 =( 1+0.25 ) √ ¿
max sin ∅0=0.559017+0.015652=0.574669

max sin ∅0=0.559017−0.015652=0.543365

( 0.574669 )∗180
max ∅0=arcsin =35.076 °
π
( 0.574669 )∗180
min ∅0=arcsin =32.913 °
π
4.18 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 segundas 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.