UNIVERSIDAD DEL VALLE SEDE TULU

ASIGNATURA: ANLISIS Y MTODOS NUMRICOS,

750086M, Sem 2017 01
TALLER No 7: Diferenciacin e integracin numrica. Solu-
cin numrica de EDO: Mtodo de Euler

Profesor: Efran Vsquez Milln

1. Diferenciacin Numrica
1. Use las frmulas de diferencia progresiva y de diferencia regresiva para determinar las aproxima-
ciones con que se complementarn las siguientes tablas

a. b.
x f (x) f 0 (x) x f (x) f 0 (x)
0 5 0 4794 0 5 0 4794
0 6 0 5646 0 6 0 5646
0 7 0 6442 0 7 0 6442

2. Los datos del ejercicio anterior se tomaron de las siguientes funciones. Calcule los errores reales del
ejercicio anterior y obtenga las cotas de error por medio de las frmulas de error.

a. f (x) = sin x b. f (x) = ex 2x2 + 3x 1

3. Use la frmula de tres puntos ms conveniente para determnar las aproximaciones con que se
complementarn las siguientes tablas.

a. c.
x f (x) f 0 (x) x f (x) f 0 (x)
1 1 9 025013 2 9 4 827866
1 2 11 02318 3 0 4 240058
1 3 13 46374 3 1 3 496909
1 4 16 44465 3 2 2 596792
b. d.
x f (x) f 0 (x) x f (x) f 0 (x)
8 1 16 94410 2 0 3 6887983
8 3 17 56492 2 1 3 6905701
8 5 18 19056 2 2 3 6688192
8 7 18 82091 2 3 3 6245909

4. Los datos del ejercicio anterior se tomaron de las siguientes funciones. Calcule los errores reales del
ejercicio anterior y obtenga las cotas de error por medio de las frmulas de error.

a. f (x) = e2x c. f (x) = x cos x x2 sin x

b. f (x) = x ln x d. f (x) = 2(ln x)2 + 3 sin x

5. Sea f (x) = cos x. Use la frmula

1h i h2
f 00 (x0 ) = 2
f (x0 h) 2f (x0 ) + f (x0 + h) f (4) ()
h 12

1
 para alguna , donde x0 h < < x0 + h y los valores de f (x) en x = 0 25, 0 5, 0 75 para aproximar
f 00 (0 5). Compare este resultado con el valor exacto y con la aproximacin conseguida en el ejercicio
17 del taller 4 (Interpolacin). Explique por qu este mtodo resulta tan exacto en este problema.
Obtenga el error.
6. En un circuito con un voltaje impreso E(t) y una inductancia L, la primera ley de Kirchhoff nos
da la siguiente relacin

di
E=L
+ Ri
dt
donde R es la resistencia del circuito e i es la corriente. Suponga que medimos la corriente con
varios valores de t y obtenemos

t 1 00 1 01 1 02 1 03 1 04
i 3 10 3 12 3 14 3 18 3 24
donde t se mide en segundos, i se da en amperes, la inductancia L es una constante de 0 98 henries y
la resistencia es de 0 142 ohms. Aproxime el voltaje E en los valores t = 1 00, 1 01, 1 02, 1 03, 1 04.

2. Integracin Numrica
7. Aplique la regla compuesta del trapecio con los valores indicados de n para aproximar las siguientes
integrales.

Z2 Z2
2
a. x ln xdx, n=4 c. dx, n=6
x2 + 4
1 0
Z2 Z5
3 x 1
b. x e dx, n=4 d. dx, n = 8
x2 + 4
2 3

8. Aplique la regla compuesta de Simpson para aproximar las integrales del ejercicio anterior.
Z2
2
9. Aproxime x2 ex dx, por medio de h = 0 25
0

a. Aplique la regla compuesta del traprecio. b. Aplique la regla compuesta de Simpson.

10. Suponga que f (0 25) = f (0 75) = . Obtenga si la regla compuesta del trapecio con n = 2 da el
Z1
valor de 2 para f (x)dx y con n = 4 da el valor 1 75.
0

11. Un automvil recorre una pista de carreras en 84 s. Su velocidad en cada intervalo de 6 s. se

determina mediante una pistola de radar y est dada, en pies/s, desde el principio del recorrido,
por los datos de la tabla siguiente:

T iempo 0 6 12 18 24 30 36 42 48 54 60 66 72 78 84
V elocidad 124 134 148 156 147 133 121 109 99 85 78 89 104 116 123

Qu longitud tiene la pista?

2
 12. Para simular las caracteraticas trmicas de los frenos de disco (vease la figura anexa), D. A. Secrist
y R. W. Hornbeck (SH) tuvieron que aproximar numricamente la "temperatura promedio exterior
promediada del rea", T , en el cojn del freno, basndose para ello en la ecuacin

Zr0
T (r)rp dr
re
T =
Rr0
rp dr
re

donde re representa el radio donde comienza el contacto

entre cojn y disco, r0 representa el radio exterior de dicho
contacto, p representa el ngulo subtendido por los cojines
del freno del sector y T (r) es la temperatura en cada punto
del cojn, la cual se obtuvo numricamente al analizar la
ecuacin del calor. Si re = 0 308 pies, r0 = 0 478 pies,
p = 0 7051 radianes y si las temperaturas dadas en la tabla
siguiente se calcularon en varios puntos del disco, obtenga
una aproximacin de T .

r (f t) 0 308 0 325 0 342 0 359 0 376 0 393 0 410 0 427 0 444 0 461 0 478
T (r) ( F ) 640 794 885 943 1034 1064 1114 1152 1204 1222 1239

3. Solucin Numrica de Ecuaciones Diferenciales Ordinarias:

Mtodo de Euler
13. Aplique el mtodo de Euler para aproximar las soluciones de los siguientes problemas de valor inicial

a. y0 = te3t 2y, 0 t 1, y(0) = 0, con h = 0 5, h = 0 2

b. y0 = 1 + (t y)2 , 2 t 3, y(2) = 1, con h = 0 5, h = 0 2
y
c. y0 = 1 + t , 1 t 2, y(1) = 2, con h = 0 25, h = 0 1
d. y0 = cos(2t) + sin(3t), 0 t 1, y(0) = 1, con h = 0 25, h = 0 1

14. A continuacin se dan las soluciones reales de los problemas de valor inicial del ejercicio anterior.
En cada paso, compare el error real con la cota de error

a. 1 3t
y(t) = 15 te3t 25 e + 251 2t
e
1
b. y(t) = t + 1t
c. y(t) = t ln t + 2t
d. y(t) = 12 sin(2t) 13 cos(3t) + 43

15. Dado el problema de valor inicial

2
y0 = y + t 2 et 1 t 2, y(1) = 0
t
con las soluciones exactas y(t) = t2 (et e):

3
a. Use el mtodo de Euler con h = 0 1 para aproximar la solucin y compararla con los valores
reales de y.
b. Use las respuestas obtenidas en la parte (a) y la interpolacin lineal para aproximar los sigu-
ientes valores de y y comprelos con los valores reales.

i. y(1 04) ii. y(1 55) iii. y(1 97)