U.L.L. Facultad de Ciencias. Grado en Matemáticas. Métodos Numéricos II. Curso 2023-2024.

Métodos Numéricos II. Grado en Matemáticas. Curso 2023-2024.

Problemas del Tema 5: Integración numérica.

1. Supongamos que f (0,25) = f (0,75) = α. Determinar α sabiendo que la regla del Trapecio
R1
para n = 2 da el valor 2 para 0 f (x)dx y que para n = 4 da el valor 1,75.
Z 1
2. Deseamos estimar integrales del tipo I(f ) = f (x)dx mediante este tipo especial de
−1
0 0
fórmula de cuadratura: J(f ) = af (−1) + bf (1) + cf (−1) + d f (1).

a) Determinar el mayor grado de precisión polinómica de J(f ) con respecto a I(f ).

b) Encontrar una expresión del error, admitiendo que la función f es suficientemente
regular.

3. Dada la tabla
xk 1 1,1 1,2 1,3 1,4
f (xk ) 0,01 0,252 0,586 1,024 1,578
Estimar el valor de la integral Z 1,4
f (x)dx,
1

mediante la aplicación de las reglas del rectángulo (derecha e izquierda), del trapecio y
de Simpson 13 .

4. Sea f (x) = |x|, x ∈ [−1, 1]. Probar que si n es par, la regla de los trapecios la integra
correctamente. Calcular la estimación si n es impar.
Rb
5. Sea Tn (f ) la regla de los trapecios para aproximar a f (x)dx. Probar que

1
T2n (f ) = Tn (f ) + h(f1 + f3 + · · · + f2n−1 ),
2
b−a
donde fi = f (xi ), xi = a + ih y h = 2n
. Como aplicación, dada la integral
Z 1,4
ln(x)dx
0,6

y trabajando con redondeo a 6 cifras decimales, aplicar sucesivamente la regla del trapecio
Tn para encontrar una aproximación a la integral anterior buscando los valores para n=
1, 2, 4 y 8 intervalos.
R1 t
6. Se pretende estimar la integral I = 0 1+t dt.

Métodos Numéricos II. Prof. Ruymán Cruz Barroso. Problemas del Tema 5: Integración numérica.
U.L.L. Facultad de Ciencias. Grado en Matemáticas. Métodos Numéricos II. Curso 2023-2024. 2

a) Utilizar las reglas del trapecio y de Simpson 13 , considerando n = 10. Comparar el

valor obtenido con el valor exacto I = 1 − ln2.
b) Para la fórmula del trapecio, obtener la expresión del error en términos de alguna
derivada de la función integrando, obtener una cota superior del error obtenido, y
obtener la expresión asintótica del error.
c) Realizar los mismo para la fórmula de Simpson 31 , sabiendo que la expresión del error
4
viene dada por En (f ) = − (b−a)h
180
f (4 (ξ), ∃ξ ∈ [a, b], y que la expresión asintótica del
4
error es en este caso En (f ) ≈ −h
 (3 
180
f (b) − f (3 (a) , siempre que f ∈ C 4 ([a, b]).

7. Hallar la fórmula de cuadratura cerrada de Newton-Cotes para n = 4.

8. Hallar las fórmulas de cuadratura abiertas de Newton-Cotes para n = 2 y n = 3.

9. a) Hallar a0 y a1 para que la fórmula de cuadratura

Z 1
f (x)xα dx = a0 f (0) + a1 f (1) + E(f ), α > −1
0

sea de tipo interpolatorio.

b) Obtener una expresión para el error en términos de alguna derivada apropiada de f.
Rh
c) Utilizando los resultados anteriores, obtener una fórmula de integración para 0 g(t)tα dt,
h > 0, incluyendo una expresión para el error.

10. Supongamos que aproximamos

Z 1
f (x)dx ≈ af (1/4) + bf (1/2) + cf (3/4).
0

Determinar a, b y c para conseguir el mayor grado de precisión posible. Deducir que la

fórmula de cuadratura tiene grado de precisión igual a 3.

11. Determinar a0 y a1 para que la fórmula de cuadratura

Z 1
f (t) dt ≈ a0 f (0) + a1 f (1)
0

sea exacta para las funciones de la forma

π 
t
f (t) = ae + bcos t , siendo a y b constantes.
2
Z 1
12. Encontrar una expresión de la forma t · f (t) dt ≈ a0 f (x0 ) + a1 f (x1 ) que sea exacta
−1
en P3 .

Métodos Numéricos II. Prof. Ruymán Cruz Barroso. Problemas del Tema 5: Integración numérica.
U.L.L. Facultad de Ciencias. Grado en Matemáticas. Métodos Numéricos II. Curso 2023-2024. 3

Z 1
13. Encontrar una expresión de la forma t · f (t) dt ≈ a0 f (x0 ) + a1 f (x1 ) + a2 f (x2 ) que sea
−1
exacta en P5 .

14. Fórmula de cuadratura de Gauss-Legendre.- Sean Ln los polinomios ortogonales de

Legendre (ω(x) ≡ 1) en [−1, 1].

a) Deducir, a partir del Problema 25 del Tema 3, la ley de recurrencia para la corres-
pondiente familia ortonormal, y obtener de ahı́ la matriz de Jacobi.
b) Computar, a partir de la matriz de Jacobi, la fórmula de cuadratura de Gauss-
Legendre de 4 puntos.
R 1 dx
c) Aprovechar el apartado anterior para obtener una estimación de −1 1+x 2 , cuyo valor
π
exacto es 2 .

15. 
Fórmula de cuadratura de Gauss-Chebyshev.- Sean Tn los polinomios de Chebyshev
2 −1/2
ω(x) = (1 − x ) en [−1, 1].

a) Probar que los coeficientes de la fórmula de cuadratura de Gauss-Chebyshev vienen

dados por wj = πn , 1 ≤ j ≤ n.
Ayuda: utilizar la igualdad
Z π
cos (nθ) sin (nθ)
dθ = π · .
0 cos θ − cos φ sin φ

Nota: se sabe que la función peso de Chebyshev de primera especie es la única para
la cual la correspondiente fórmula Gaussiana tiene todos sus pesos iguales.
b) Describir dicha fórmula, dando además una expresión para el error.
R1
16. Obtener la fórmula de cuadratura Gaussiana con 2 nodos para I(f ) = 0 f (x) ln( x1 )dx,
considerando la función peso ω(x) = ln x1 . Dar una expresión del error.

17. Encontrar quiénes tienen que ser A1 , A2 , x1 y x2 para que la fórmula siguiente sea exacta
para polinomios de hasta grado cúbico (fórmula de Gauss-Laguerre con 2 nodos,
α = 0):
Z ∞
e−x g(x)dx ≈ A1 g(x1 ) + A2 g(x2 ).
0

Obtener una expresión del error y utilizar la fórmula para aproximar las integrales
Z ∞ −x Z ∞
e −x


I1 = dx, I 2 = log 1 + e dx.
0 1 + x2 0

Métodos Numéricos II. Prof. Ruymán Cruz Barroso. Problemas del Tema 5: Integración numérica.
U.L.L. Facultad de Ciencias. Grado en Matemáticas. Métodos Numéricos II. Curso 2023-2024. 4

18. Un notable resultado de la Teorı́a Analı́tica de Números establece que una estimación
para la cantidad de números primos que hay en (a, b) es
Z b
dt
.
a ln(t)
1 3
a) Utilizar las fórmulas de cuadratura simples de Simpson 3
y 8
para estimar la cantidad
de números primos en el intervalo (100, 200).
b) Repetir empleando ahora la fórmula de cuadratura de Gauss-Legendre considerando
tres nodos.
c) Dar cotas de los errores cometidos en los apartados anteriores. Comparar con el error
real sabiendo que la solución es 21.

19. Experimente con diferentes formas de calcular las integrales

Z 1
cos(t)
Ic = √ dt = 1,80904847580054...
0 t
Z 1
sen(t)
Is = √ dt = 0,62053660344676...
0 t
a) Utilice la regla de los trapecios con 4 subintervalos de igual longitud en [0,01, 1], con
el fin de evitar la singularidad en x = 0.
b) Utilice la regla de los trapecios con 3 subintervalos de longitud h = 1/4 en el intervalo
[h, 1], en combinación con la fórmula de Newton-Cotes ponderada en el intervalo [0, h]
con el peso ω(t) = t−1/2 .
c) Realice el cambio de variable t = x2 y aplique la regla de los trapecio con 4 intervalos
de longitud 14 en [0, 1].
d ) Utilice la cuadratura de Gauss-Legendre con 5 nodos para las integrales obtenidas
en el apartado (c).
e) Utilice la cuadratura de Gauss-Jacobi con 5 nodos y con parámetros α = 0 y β = − 12
directamente para las integrales Ic y Is.

Compare los resultados.

20. Hallar las fórmulas de cuadratura con cinco nodos

a) de Gauss-Chebyshev.
b) de Gauss-Laguerre (con parámetro α = 0).
c) de Gauss-Hermite.

Métodos Numéricos II. Prof. Ruymán Cruz Barroso. Problemas del Tema 5: Integración numérica.
U.L.L. Facultad de Ciencias. Grado en Matemáticas. Métodos Numéricos II. Curso 2023-2024. 5

Aplicar la fórmula que se considere más adecuada para estimar las integrales
1 2 ∞ ∞
et
Z Z Z
t−t2

2
I1 = √ dt, I2 = e dt, I3 = sin t2 e−t dt.
−1 1 − t2 0 −∞

21. Considere una fórmula de cuadratura del tipo

Z ∞
e−t · f (t)dt = af (0) + bf (c) + E(f ).
0

a) Encontrar a y b para que la fórmula tenga el grado de exactitud igual a 2.

b) Sea H2 (f ; 0, 2, 2; x)Zel polinomio de Hermite que interpola a la función f en los nodos
∞
{0, 2, 2}. Calcular e−t H2 (f ; 0, 2, 2; x)dt y comparar con el apartado anterior.
0
c) Encontrar una expresión del error en términos de la derivada tercera de f .

22. Usar la fórmula de cuadratura de Gauss-Laguerre para n = 2, 4, 6 y 8 nodos para

estimar las siguientes integrales:
Z ∞ √
−t2 π
a) e dt = .
0 2
Z ∞
tdt 1
b) 2 2
= .
0 (1 + t ) 2
Z ∞
sen(t) π
c) dt = .
0 t 2

Métodos Numéricos II. Prof. Ruymán Cruz Barroso. Problemas del Tema 5: Integración numérica.