Facultad de Ingeniería

Escuela de Ingeniería Civil

GUÍA DE APRENDIZAJE
SEMANA N° 09

Curso: Métodos Numéricos

Docente: Mg. Nilthon Arce Fernández

Jaén – Perú, diciembre 2023

0
 Facultad de Ingeniería
Escuela de Ingeniería Civil

INDICE

1. INTRODUCCIÓN .....................................................................................................................2
2. OBJETIVOS EDUCACIONALES Y RESULTADOS DEL ESTUDIANTE ..........................2
3. DESARROLLO .........................................................................................................................2
3.1 DIFERENCIACIÓN..............................................................................................................2
3.2 INTEGRACIÓN ....................................................................................................................3
3.3 FÓRMULAS DE INTEGRACIÓN DE NEWTON-COTES................................................4
3.4 LA REGLA DEL TRAPECIO .............................................................................................. 5
4. ACTIVIDAD Y EVALUACIÓN .............................................................................................10
4.1 Actividad ..............................................................................................................................10
4.2 Evaluación............................................................................................................................ 10
5. GLOSARIO .............................................................................................................................11
6. REFERENCIAS.......................................................................................................................11

1
 Facultad de Ingeniería
Escuela de Ingeniería Civil
1. INTRODUCCIÓN
El cálculo es la matemática del cambio. Como los ingenieros deben tratar en forma continua
con sistemas y procesos que cambian, el cálculo es una herramienta esencial en nuestra
profesión. En la esencia del cálculo están dos conceptos matemáticos relacionados: la
diferenciación y la integración (Chapra et al., 2007).

2. OBJETIVOS EDUCACIONALES Y RESULTADOS DEL ESTUDIANTE

Esta semana la guía de aprendizaje contribuye con el logro de las siguientes capacidades:
▪ Conocer y comprender las fórmulas de integración de Newton – Cotes.
▪ Valorar la importancia de la regla del trapecio en la integración numérica.

3. DESARROLLO
3.1 DIFERENCIACIÓN
Diferenciar significa “marcar por diferencias; distinguir; percibir la diferencia en o entre”. En
el contexto de las matemáticas, la derivada sirve como el principal vehículo para la
diferenciación, representa la razón de cambio de una variable dependiente con respecto a una
variable independiente. La definición matemática de la derivada empieza con una
aproximación por diferencias:
y f ( xi + x) − f ( xi )
=
x x
donde y y f ( x) son representaciones alternativas de la variable dependiente y x es la

variable independiente. Si se hace que x se aproxime a cero, el cociente de las diferencias

se convierte en una derivada
y f ( xi + x) − f ( xi ) dy f ( xi + x) − f ( xi )
lim = lim → = lim
x →0 x x → 0 x dx x → 0 x

Figura 1. La definición gráfica de una derivada: conforme x se aproxima a cero al ir de a) a c), la

aproximación por diferencias se va convirtiendo en una derivada.
Como se observa en la descripción visual de la Figura 1c), la derivada evaluada es la pendiente
de la recta tangente a la curva en xi .

2
 Facultad de Ingeniería
Escuela de Ingeniería Civil
3.2 INTEGRACIÓN
En cálculo, el proceso inverso de la diferenciación es la integración. Integrar significa “juntar
partes en un todo; unir; indicar la cantidad total”. Matemáticamente, la integración se
representa por:
b
I =  f ( x)dx
a

que representa la integral de la función f ( x) con respecto a la variable independiente x ,

evaluada entre los límites x = a y x = b . La función f ( x) se llama integrando.

Figura 2. Representación gráfica de la integral de f ( x) entre los límites x = a y x = b . La

integral es equivalente al área bajo la curva.

Como se dijo antes, la “distinción” o “discriminación” de la diferenciación y el “juntar” de la

integral son procesos estrechamente relacionados, de hecho, inversamente relacionados. Por
ejemplo, si se tiene una función dada y (t ) que especifica la posición de un objeto en función
del tiempo, la diferenciación proporciona un medio para determinar su velocidad
d
v(t ) = y (t ) = y ' (t )
dt
De manera inversa, si se tiene la velocidad como una función del tiempo, la integración se
utilizará para determinar su posición
t
y (t ) =  v(t )dt
0

3
 Facultad de Ingeniería
Escuela de Ingeniería Civil
3.3 FÓRMULAS DE INTEGRACIÓN DE NEWTON-COTES
Las fórmulas de Newton - Cotes son los tipos de integración numérica más comunes. Se basan
en la estrategia de reemplazar una función complicada o datos tabulados por un polinomio de
aproximación que es fácil de integrar:
b b
I =  f ( x)dx   f n ( x)dx
a a

Donde f n ( x) = un polinomio de grado " n " de la forma

f n ( x) = a0 + a1 x + a2 x 2 + ... + an−1 x n−1 + an x n (1)

Figura 3. La aproximación de una integral mediante el área bajo a) una sola línea recta y b) una
parábola.
En la Figura 3a), se utiliza un polinomio de primer grado (una línea recta) como una
aproximación. En la Figura 3b), se emplea una parábola con el mismo propósito.
La integral también se puede aproximar usando un conjunto de polinomios aplicados por
pedazos a la función o datos, sobre intervalos de longitud constante. Por ejemplo, en la Figura
4, se usan tres segmentos de línea recta para aproximar la integral. Aunque pueden utilizarse
polinomios de grado superior con los mismos propósitos.

Figura 4. La aproximación de una integral mediante el área bajo tres segmentos de línea recta.
4
 Facultad de Ingeniería
Escuela de Ingeniería Civil
3.4 LA REGLA DEL TRAPECIO
La regla del trapecio es la primera de las fórmulas cerradas de integración de Newton - Cotes.
Corresponde al caso donde el polinomio de la ecuación (1) es de primer grado:
b b
I =  f ( x)dx   f1 ( x)dx
a a

La forma más simple de interpolación consiste en unir dos puntos con una línea recta. Dicha
técnica, llamada interpolación lineal, se define como
f ( x1 ) − f ( x0 )
f1 ( x) = f ( x0 ) + ( x − x0 )
x1 − x0
Entonces, la línea recta entre los límites a y b se puede representar como
f (b) − f (a)
f1 ( x) = f (a) + ( x − a)
b−a
El área bajo esta línea recta es una aproximación de la integral de f ( x)

 f (b) − f (a)  f (a) + f (b)

b b
I   f1 ( x)dx =   f (a ) + ( x − a ) dx = (b − a )
a a 
b−a  2

f (a) + f (b)
 REGLA DEL TRAPECIO: I  (b − a)
2

Figura 5. Representación gráfica de la regla del trapecio.

Una forma de mejorar la precisión de la regla del trapecio consiste en dividir el intervalo de
integración de a a b en varios segmentos, y aplicar el método a cada uno de ellos (Figura
6). Las áreas de los segmentos se suman después para obtener la integral en todo el intervalo.
Las ecuaciones resultantes se llaman fórmulas de integración, de aplicación múltiple o
compuestas (Chapra et al., 2007).

5
 Facultad de Ingeniería
Escuela de Ingeniería Civil

Figura 6. Ilustración de la regla del trapecio de aplicación múltiple, a) dos segmentos, b) tres
segmentos, c) cuatro segmentos y d) cinco segmentos.

La Figura 7, muestra el formato general y la nomenclatura que usaremos para obtener

integrales de aplicación múltiple. Hay n +1 puntos igualmente espaciados
( x0 , x1 , x2 , x3 ,..., xn ) . En consecuencia, existen n segmentos del mismo ancho

Figura 7. Formato general y nomenclatura para integrales de aplicación múltiple.

6
 Facultad de Ingeniería
Escuela de Ingeniería Civil

Si x0 = a y xn = b , la integral completa se representará como

x1 x2 x3 xn−1 xn

I =  f ( x)dx +  f ( x)dx +  f ( x)dx +... +  f ( x)dx +  f ( x)dx

x0 x1 x2 xn−2 xn−1

Sustituyendo la regla del trapecio en cada integral se obtiene

f ( x0 ) + f ( x1 ) f ( x1 ) + f ( x2 ) f ( x2 ) + f ( x3 ) f ( xn −2 ) + f ( xn −1 ) f ( xn −1 ) + f ( xn )
I h +h +h + ... + h +h
2 2 2 2 2
h n −1
 b−a
I 
2
f ( x0 ) + 2 
i =1
f ( xi ) + f ( xn )  , h =
 n
, xi = x0 + ih

 n −1

 f ( x0 ) + 2  f ( xi = x0 + ih) + f ( xn ) 
I  (b − a)  i =1

 2n 
 

Se tiene un error con la regla del trapecio de aplicación múltiple al sumar los errores
individuales de cada segmento, así

(b − a)3 n ´´
Et = −  f (i )
12n3 i =1
donde f ´´ (i ) es la segunda derivada en un punto i , localizado en el segmento i . Este

resultado se simplifica al estimar la media o valor promedio de la segunda derivada en todo el

intervalo como
n
− f ´´
(i ) n −

f ´´  i =1
→  f (i )  n
´´
f ´´
n i =1

luego:
−
(b − a)3
Ea = − n f ´´
12n3
−
(b − a)3
Ea = −
12n 2 f ´´

donde:
b

f
´´
− ( x)dx
´´ ( x) = a
f b−a

7
 Facultad de Ingeniería
Escuela de Ingeniería Civil
Ejemplo 1. Use la regla del trapecio con dos segmentos para estimar la integral de
f ( x) = 0.2 + 25x − 200 x2 + 675 x3 − 900 x4 + 400 x5
desde a = 0 hasta b = 0.8 . Estimar el error aproximado, recuerde que el valor correcto para
la integral es I = 1.640533 .
Solución
a) Hallando la integral aproximada
xn − x0 b − a 0.8 − 0
h= = = = 0.4
n n 2
 n −1

 f (a) + 2 f ( xi = x0 + ih) + f (b)   f (0) + 2 f ( x1 = x0 + h) + f (0.8)   f (0) + 2 f ( x1 = 0.4) + f (0.8) 
I  (b − a)  i =1
 = (0.8 − 0)   = 0.8   (*)
 2n   4   4
 

Pero, f (0) = 0.2 , f (0.4) = 2.456 y f (0.8) = 0.232

 0.2 + 2(2.456) + 0.232 

Remplazando en (*) : I  0.8   → I  1.0688
 4
b) Hallando el error aproximado
− −
(b − a)3 0.512
Ea = −
12n 2 f ´´ = − 48 f ´´ (**)
f ´( x) = 25 − 400 x + 2025 x 2 − 3600 x3 + 2000 x 4
f ´´( x) = −400 + 4050 x − 10800 x 2 + 8000 x3
b 0.8 0.8

 f ´´ ( x)dx  f ´´ ( x)dx  (−400 + 4050 x − 10800 x + 8000 x 3 ) dx

2
−

f ´´ ( x) = a

b−a
= 0

0.8
= 0

0.8
= −60

0.512
Remplazando en (**) : Ea = − (−60) = 0.64
48
Error verdadero: Et = valor verdadero − valor aproximado = 1.640533 − 1.0688 = 0.57173

Ejemplo 2. Escribir el código en Matlab para el método del trapecio con la siguiente interfaz:
function[ A] = trapecio( fun, a, b, n)
donde fun es la función que se desea integrar, a y b son los límites inferior y superior de

integración, n es el número de divisiones del intervalo [a, b] que usa el método del trapecio.
4
Probarlo con la siguiente integral  xsen( x + 2)dx , para n = 8. (Díaz et al., 2015)
0

Solución

8
 Facultad de Ingeniería
Escuela de Ingeniería Civil
a) Código en Matlab para el método del trapecio
function[ y ] = fun( x)
y = x.*sin( x + 2);
end

function[ A] = trapecio( fun, a, b, n)

h = (b − a ) / n;
s = 0;
for i = 1: n − 1
s = s + feval ( fun, a + i * h);
end
A = ( feval ( fun, a) + 2* s + feval ( fun, b))* h / 2;
end
4
b) Probando el código con  xsen( x + 2)dx , para n = 8.
0

4
I =  xsen( x + 2)dx = −4.9741
0

9
 Facultad de Ingeniería
Escuela de Ingeniería Civil
2
Ejemplo 3. Calcular la siguiente integral  ln( x)dx , usando el método del trapecio con n = 6
1

Solución
xn − x0 b − a 1
h= = =
n n 6
 n −1

 f ( a ) + 2  f ( xi = x0 + ih) + f (b) 
 f (1) + 2[ f ( x1 ) + f ( x2 ) + f ( x3 ) + f ( x4 ) + f ( x5 )] + f (2) 
I  (b − a)  i =1
=  ()
 2 n   12 
 

f (1) = 0 y f (2) = 0.6931

xi = a + ih; i = 1, 2,3, 4,5
Luego, f (7 / 6) = 0.1542 , f (4 / 3) = 0.2877 , f (3 / 2) = 0.4055 , f (5 / 3) = 0.5108 y

f (11/ 6) = 0.6061
Remplazando en () :

 0 + 2[0.1542 + 0.2877 + 0.4055 + 0.5108 + 0.6061] + 0.6931 

I   → I  0.3851
 12

4. ACTIVIDAD Y EVALUACIÓN
4.1 Actividad
2
1. Hallar la siguiente integral  x ln( x + 1)dx , usando el método del trapecio con n = 8.
0

e
− x2
2. Usando el método del trapecio, hallar la integral dx , con n = 8.
−1

 [1 + sen( x )]dx , usando el método del trapecio con n = 4.

2
3. Hallar la siguiente integral
−

4. Verificar los resultados de los ejercicios anteriores utilizando códigos en Matlab

5. Hacer un resumen indicando las partes más importantes de la presente guía de aprendizaje

Nota: La fecha de presentación de las actividades es el día domingo 31 de diciembre de

8:00am a 11:59pm, utilizar el formato PDF para la presentación, no exceder los 3 MB de peso.
4.2 Evaluación
Para evaluar esta actividad se utilizará la siguiente Rúbrica.

10
 Facultad de Ingeniería
Escuela de Ingeniería Civil

5. GLOSARIO

6. REFERENCIAS
Chapra, S. C., Canale, R. P., Ruiz, R. S. G., Mercado, V. H. I., Díaz, E. M., & Benites, G. E.
(2007). Métodos numéricos para ingenieros (Vol. 5). McGraw-Hill México.
Díaz, L., Mestas, R., & Hancco, R. (2015). Métodos Numéricos con Matlab.

11