Universidad del Bo-Bo
Facultad de Ciencias
Departamento de Matematicas
Ejercicios Integracion Numerica
Ejercicio 1:
Aproxime el area bajo la curva f (x) = 1 + 2x + 3x2 entre [-1,0], usando el Metodo del
Trapecio con h = 0,1.
Recordemos el Metodo del Trapecio es de la forma:
n1
h X
I(f ) = [f (a) + 2 f (xj ) + f (b)] E
2 j=1
donde E es el error, de la forma:
(b a) 2 (2)
E= h f (), [a, b]
12
y
xj = a + j h, j = 0...n
Desarrollo:
Ahora lo primero que tenemos que obtener es n, lo cual es bastante simpre pues h es
conosido, Si h = 0,1:
ba 0 (1)
h= 0,1 =
n n
n = 10
luego obtenemos xj para realizar la iteracion:
xj = a + j h, j = 0...n
xj = 1 + 0,1j, j = 0, 1,2..,10
esto es:
j xj f (xj )
0 -1 2
1 -0.9 1.63
2 -0.8 1.32
3 -0.7 1.07
4 -0.6 0.88
5 -0.5 0.75
6 -0.4 0.68
7 -0.3 0.67
8 -0.2 0.72
9 -0.1 0.83
10 0 1
Por lo que tenemos que:
f (a) = 2, f (b) = 1 y
n1
X 9
X
2 f (xj ) = 2 f (xj ) = 2[1,63+1,32+1,07+0,88+0,75+0,68+0,67+0,72+0,83] = 28,55 = 17,1
j=1 j=1
finalmente la aproximasion de la integracion:
) = 0,1 [2 + 17,1 + 1] = 1,005
I(f
2
por otro lado: Z 1
1 + 2x + 3x2 = x + x2 + x3 /01 = 1
0
Entonces la aproximasion de la aproximasion no esta lejos del valor real 1 es casi 1.005
pero para ver si esto es correcto calcularemos el E.
1 0,01 (2)
E= (0,1)2 f (2) () = f ()
12 12
= 0,00083f (2) ()
donde:
f 00 (x) = 6
luego:
|E| = |0,00083 6| = 0,005
entonces se cumple que el error verdadero no puede superar la cota del error, es decir:
0,005 |E|
Ejercicio 2:
(ba)
Para f C 2 [a, b], h = (2m+2)
, xj = a + (j + 1)h y para todo j = 1, 0, 1, ..., 2m + 1, se
tiene m
b
b a 2 00
Z X
f (x)dx = 2h f (x2j ) + h f () (a, b)
a j=0
6
Este metodo es llamado Regla del Punto Medio, con n = 2m subintervalos.
Definiendo Z b
I(f ) = f (x)dx
a
m
X
) = 2h
I(f f (x2j )
j=0
a)Determine el valor de n N tal que para un > 0 dado, se tenga
)| <
|I(f ) I(f
Indicacion: M = max x[a,b] |f 00 (c)|.
b)Dada la integral Z 1,5
1
I(f ) = dx
0 1+x
Determine el valor de n tal que la integral anterior difiera de su aproximacion sea a lo
mas de 102
),con n=10.
c)Calcule I(f
Desarrollo:
a)Tenemos que:
)| =
|I(f ) I(f ba 2 00
h f (), donde:
6
1 2
f 0 (x) = 2
, f 00 (x) =
(1 + x) (x + 1)3
Luego
b a 2 00
|E| = h f (), con M = f 00 ()
6
Se tiene
2
max = , a 6= 1
M [a,b] (a + 1)3
Entonces:
ba 2
|E| = h M <
6
ba 2 2 ba
= h < , para h =
6 (a + 1)3 n+2
(b a) (b a)2 1
= <
3 (n + 2) (a + 1)3
2
(b a)3 1 2
p
= < (n + 2) / () al ser natural no cambia la desigualdad
3 (a + 1)3
s
(b a)3 1
= <n+2
(a + 1)3 3
s
(b a)3
= 2<n
3(a + 1)3
Por lo tanto s
(b a)3
n= 2, con n N
3(a + 1)3
R 1,5 1
b)Tenemos 0 1+x con = 102 , a = 0, b = 1,5
Por lo que ahora lo reemplezaremos en el valor obtenido en a):
s
(1,5 0)3
2<n
3 102 (0 + 1)3
p
112,5 2 < n
8,606601 < n
Por lo tanto para que la diferencia de su aproximacion sea a lo mas de 102 , n debe
ser de(por redondeo):
n=9
c)Usaremos n = 2m, por lo que tenemos:
ba 1,5 0
h= = = 0,125
n+2 10 + 2
ademas
xj = a + (j + 1)h, j = 1, 0, 1...n + 1
xj = 0 + 0,125(j + 1), j = 1, 1..,11
Ahora tenemos la tabla:
j xj f (xj )
0 -1 2
-1 0 1
0 0.125 8/9
1 0.25 4/5
2 0.375 8/11
3 0.5 2/3
4 0.625 8/13
5 0.75 4/7
6 0.875 8/15
7 1 1/2
8 1.125 8/17
9 1.25 4/9
10 1.375 8/19
11 1.5 2/5
Entonces la:
n
2 5
X X 8 8 8 8 8 8
f (x2j ) = f (x2j ) = + + + + + = 3,6565204
j=o j=o
9 11 13 15 17 19
Ahora n
2
X
) = 2h
I(f f (x2j ) = 2 0,125 3,6565204 = 0,9141301
j=o
Por lo tanto la aproximasion a la integral es
) = 0,9141301
I(f
Por otro lado Z 1,5
1
= ln |1 + x|/01,5 = 0,91629073
0 1+x
Finalmente podemos decir que estamos frente a una aproximacion adecuada.
Ejercicios Diferenciacion Numerica
Reglas para la dfirenciacion
Regla de dos puntos
Si f C (2) [a, b],se tiene que existe (x0 , x0 + h) tal que
f (x0 + h) f (x0 ) h 00
f 0 (x0 ) = f () (1)
h 2
Se llama formula de la diferencia progresiva si h > 0 y formula de diferencia regre-
siva si h < 0
Regla de tres puntos
Si f c(3) (I), se tiene que existen 1 (x0 , x0 + 2h) y 2 (x0 , x0 2h)tal que
3f (x0 ) + 4f (x0 + h) f (x0 + 2h h2 (3)
f 0 (x0 ) = + f (1 ) (2)
2h 3
y
f (x0 + h) f (x0 h) h2 (3)
f 0 (x0 ) = f (2 ) (3)
2h 6
Regla de cinco puntos
Si f c(4) (I), se tiene que existen 3 (x0 2h, x0 + 2h) y 4 (x0 , x0 + 4h)tal
que
f (x0 2h) 8f (x0 h) + 8f (x0 + h) f (x0 + 2h h4 (5)
f 0 (x0 ) = + f (3 ) (4)
12h 30
y
25f (x0 ) + 48f (x0 + h) 36f (x0 + 2h) + 16f (x0 + 3h) 3f (x0 + 4h) h4 (5)
f 0 (x0 ) = + f (4 )
12h 5
(5)
Ejercicio 1
Con la regla de dos y tres puntos encuentre una a proximacion a la derivada de la
funcion f (x) = Ln(x) en x0 = 1,8 con h = 0,1 0,01 0,001
Para la regla de dos puntos (1) tenemos
f (1,8 + h) f (1,8)
f(1,8) =
h
con la cota del error
h 00
|f 0 (1,8) f(1,8)| = |f ()|
2
donde
1 1
f 0 (x) = f 00 (x) =
x x2
entonces
h 1
|f 0 (1,8) f(1,8)| = | 2 |
2
h 1
| 2|
2 1,8
pues tomando el maximo se tiene
1,8 < < 1,8 + h
1 1 1
< <
1,8 + h 1,8
Por lo tanto la cota del error para la regla de dos puntos es
h
6,48
Por otro lado la cota del error a partir de la regla de tres puntos (3) se obtiene con
el error
h2
|E| = f (3) (2 )
6
donde
1 1 2
f 0 (x) = f 00 (x) = 2 f (3) = 3
x x x
y
x [1,8 h , 1,8 + h]
es decir
1 1 1
1,8 h < x < 1,8 + h < <
1,8 + h x 1,8 h
1 1 1
( )3 < ( )3 < ( )3
1,8 + h x 1,8 h
2 2 2
( )3 < ( )3 < ( )3
1,8 + h x 1,8 h
por lo que el
2 2
max | 3
|=
x (1,8 h)3
luego la cota del error a partir de tres puntos es
h h2 2
|f 0 (1,8) f0 (1,8)| = max |f 00 ()| =
2 6 (1,8 h)3
Ahora finalmente obtenemos el valor de las derivadas por (1) y (3)
dos puntos error(1) tres puntos error(3)
f (1,8+h)f (1,8) h f (1,8+h)f (1,8h) h2 1
h h 2(1,8)2 2h 3 (1,8h)3
4
0.1 0.5406722 0.015434 0.556128 6.784710
0.01 0.5540180 0.0015432 0.555561 5.8117106
0.001 0.5554013 0.0001534 0.555555 5.7251108
por otro lado la derivada en el punto
f 0 (1,8) = 0.5
Finalmente podemos concluir que el metodo de tres puntos obtiene una aproximacion
mas rapida y con menor error que el de dos puntos.
Ejercicio 2
Los valores de f (x) = xex estan dados en la tabla siguiente
x 1.8 1.9 2.0 2.1 2.2
f (x) 10.889365 12.703199 14.778112 17.148957 19.855030
a)Aproxime el valor de f 0 (2,0) usando la regla de dos puntos.
b)Aproxime el valor de f 0 (2,0) usando regla de tres puntos.
c)Considere el verdadero valor de f 0 (2,0) y determine el error absoluto respecto a los
dos items anteriores y diga cual es mejor.
Desarrollo
Primero podemos ver que dados los intervalos tenemos h = 0,1, por lo que por la
regla de dos puntos (1)
f (2,1) f (2)
f0 (2) =
0,1
con los daos de la tabla
17,148957 14,778112
f0 (2) = = 23,708450
0,1
ahora con 3 puntos(3)
f (2,1) f (1,9)
f0 (2) =
0,2
17,148957 12,703199
=
0,2
= 22,228790
por otro lado el valor real de la derivada es
f 0 (x) = (x + 1)ex
f 0 (2) = 3e2 = 22,167168
lo que nos lleva a los errores
|E| = |f 0 (x) f0 (x)|
para dos puntos
|E| = |22,167168 23,708550| = 1,541282
para tres puntos
|E| = |22,167168 22,228790| = 0,061622
Ahora veremos el error absoluto:
Para dos puntos
0 0 h2 00
|f (x) f (x)| = f ()
3
donde
f 00 (x) = (x + 2)ex
luego
h2 00
|E| = f ()
3
h2 00
con M = 3
f (), por lo que
max |M | = 4,1e2,1 33,481210
x[2;2,1]
entonces
(0,1)2
|E| = 33,481210 0,11160403
3
Para tres puntos
h2
|E| = |f 0 (x) f0 (x)| = f (3) ()
6
donde
f (3) (x) = (x + 3)ex
con
M = f (3) ()
por lo que
max |M | = 5,1e2,1 41,647417
x0 [1,9;2,1]
entonces
0,12
|E| = 41,647417 0,06941245
6
A simple vista observamos que por el metodo de tres puntos obtenemos una aproxi-
macion mas certera, lo que ratificamos con su error absoluto y real, pues son menores
que los del metodo de dos puntos.
Rodrigo Melendez
Ayudantias Metodos Numericos
