Facultad de Ingeniería

EAP INGENIERÍA CIVIL

UNIDAD TEMÁTICA N° 6: CUADRATURA

1. Calcular usando el método de trapecio y Simpson con 5 cifras decimales

6.5  x ln x  ln x x   sen 2  x  

a.  5

 x 1
dx para n = 8

e.  0  1  cos2 x  dx para n = 8
 
2.5  tanx  1  
 sen x 2  
b.   dx para n = 8
 ln x  1.5 
f.  2  dx para n = 8
 sen  x  
1 e
0 0

 senx    esen x  x 2 
2


1
c.  dx para n = 8 g.  
1.7
dx para n = 8
0 4 
 1 x  0.2  x 
 
 ex 
2

d.   dx para n = 8  x  212  ln x  3 

1

01 x  h.   dx para n = 8

1
 1  1 x 2 
 

2. Indica todas las instrucciones y halla el valor de la integral en [-1, 1] de la función f x  e
 ,
cos x 2

utilizando el programa de la fórmula de Simpson compuesta con 20, 40 y 80 subintervalos. ¿Sería posible
calcularlo con 75 subintervalos? ¿Por qué?
3. Dada la función f  x   ex x1 , indica las instrucciones necesarias para calcular la integral en [-2, 2]
2

utilizando la regla trapezoidal y 50, 100 y 200 subintervalos.

5x
4. Dada la función f  x   3 , indica las instrucciones necesarias para calcular la integral en [1, 2]
x
utilizando la fórmula de la extrapolación de Richardson y 10, 20 y 40 subintervalos.
5. Utiliza el método de Simpson iterativo para aproximar la integral en [0, 4] de la función f  x   x 2  4x  3
con un error menor de 10-4. Indica las instrucciones, resultado y número de iteraciones que han sido
necesarias partiendo de 2 subintervalos.
6. Utiliza el método de trapecios iterativo para aproximar la integral en [/12,  /4] de la función
f  x   Cosec  2 x  Ctg  2 x  con un error menor de 10-4. Indica las instrucciones, resultado y número de
iteraciones que han sido necesarias partiendo de 2 subintervalos.
dy 1
7. Considere la ecuación diferencial 
dt 1  t4
Con la condición de contorno y(0)=0. Intégrese esta ecuación diferencial hasta t = 0.5 con pasos de
integración h = 0.1, usando el método de Euler. Compárese el resultado con la evaluación numérica de la
0.5


1
integral y  0.5  
0 1  t4

8. Dada la tabla de valores:

x -1 -0.9 -0.8 -0.7 -0.6 -0.5 -0.4 -0.3 -0.2 -0.1 0

y 0.37 0.49 0.64 0.83 1.04 1.28 1.55 1.84 2.14 2.44 2.72

Indica las instrucciones que utilizarías para hallar la integral usando la fórmula de Simpson simple con
los valores adecuados y después la compuesta con todos los valores. Calcula en ambos casos el resultado.

Prof. Carlos Deudor Gomez 1 Métodos Numéricos

9. A partir de la siguiente tabla de valores
x 1 2 3 4 5 6 7
f(x) 2.0000 4.2500 9.1111 16.0625 25.0400 36.0277 49.0204

7
Determinar la mejor aproximación posible al valor de la integral
1
f  x  dx

PROBLEMAS
1
10. Un cuerpo está sometido a la fuerza conservativa F  x   x  y se desplaza desde la posición x = 1 hasta
x
x = 1.8. Determinar el trabajo realizado utilizando la regla trapezoidal con pasos de h = 0.2, 0.4 y 0.8.
Obtener asimismo, una estimación del error del mejor resultado.

11. Se tiene un tanque esférico de radio r = 5m, la velocidad de salida por el orificio del fondo es 4.895
raíz(h)m/s, el diámetro de dicho orificio es 10 cm. Si el tanque tiene inicialmente un nivel de agua de 4
m, calcular el tiempo requerido para que el nivel de agua sea h=3.
dh 0.122375 h
Este problema se puede resolver a partir de la siguiente ecuación  .
dt 10h  h 2
12. Utiliza la fórmula de Simpson compuesta, con h=1, para aproximar el volumen del solido de revolución
cuyo radio R(x) viene dado por
x 0 1 2 3 4 5 6
R(x) 6.2 5.8 4.0 4.6 5.0 7.6 8.2

13. Suponga que se desea encontrar el área de un terreno de forma irregular acotado por un camino recto y
la orilla de un lago. Los límites del terreno se indican mediante las líneas trazadas en la figura (a). Suponga
que la frontera de 1 milla a lo largo del camino se divide, por ejemplo, en n = 8 subintervalos y que luego,
como se muestra en la figura (b)
Considere que
Mediciones
esta curva es la
en pies
gráfica de alguna
función y = f(x).

14. Use los datos mostrados en la figura y la regla de Simpson para encontrar una aproximación al área bajo
la gráfica de la función continua f sobre el intervalo 0, 4

Prof. Carlos Deudor Gomez 2 Métodos Numéricos

15. Use la regla trapezoidal con n=9 para encontrar una aproximación al área bajo la gráfica en la
figura ¿Proporciona esta regla el valor exacto del área?

16. El gran estanque para peces en forma irregular que se muestra en la figura contiene agua hasta una
profundidad uniforme de 4 pies. Use la regla de Simpson para encontrar una aproximación al número de
galones de agua que hay en el estanque. Las medidas están en pies; la separación vertical entre mediciones
horizontales es 1.86 pies. En 1 pie3 de agua hay 7.48 galones de agua.

17. El momento de inercia de I de la hélice con tres aspas de un barco cuya dimensión se proporcionan en la
3  3  4.5 2
g 1
figura está dado por I   r Adr , donde  es la densidad del metal, g es la aceleración de la
2g
gravedad y A es el área de una sección transversal de la hélice a una distancia de r pies del centro del

Prof. Carlos Deudor Gomez 3 Métodos Numéricos

 cubo. Si   570 lb/pie para el bronce, use los datos en la figura y la regla trapezoidal para encontrar
3

una aproximación a I.

r(pies) 1 1.5 2 2.5 3 3.5 4 4.5

A(pies) 0.3 0.5 0.62 0.7 0.6 0.5 0.27 0

18. En el siguiente gráfico se muestra delineada la zona de un derrame de petróleo ocurrido en cierta región.
Las mediciones han sido obtenidas a distancias de 4 Km.

Con la fórmula de Simpson, encuentre en forma aproximada el área total cubierta por el derrame de
petróleo.

1 x
19. La siguiente definición se denomina función error: erf  x  
 
2
e t dt .
0

Esta función se usa para definir la función de densidad de probabilidad de la distribución normal
estándar. Calcule erf(1) con la fórmula de Simpson con m=2,4,6. Con estos resultados estime la precisión
en la respuesta.
20. Una empresa está vendiendo una licencia para manejo de un nuevo punto de venta. La experiencia indica
que dentro de t años, la licencia generará utilidades según f(t) = 14000+490t dólares por año. Si la tasa de
interés r permanece fija durante los próximos n años, entonces el valor presente de la licencia se puede
n
 f t  e
 rt
calcular con V  dt
0

Calcule V si n = 5 años y r = 0.07. Use la fórmula de Simpson con m = 4,6,8. Con estos resultados estime
el error de truncamiento
21. Un comerciante usa el siguiente modelo para describir la distribución de la proporción x del total de su
mercadería que se vende semanalmente: f  x   20x 3  1  x  , 0x1.
El área debajo de f(x) representa entonces la probabilidad que venda una cantidad x en cualquier semana.

Prof. Carlos Deudor Gomez 4 Métodos Numéricos

 Se desea conocer la probabilidad que en una semana venda más de la mitad de su mercadería (debe
integrar f entre 0.5 y 1). Use la Fórmula de Simpson con m=4
22. Según la teoría de Kepler, el recorrido de los planetas es una elipse con el Sol en uno de sus focos. La


 a 1  k 2 sen 2  t dt
2
longitud del recorrido total de la órbita está dada por s  4
0

Siendo k: excentricidad de la órbita y a: longitud del semieje mayor.

Calcule el recorrido del planeta Mercurio sabiendo que k=0.206, a=0.387 UA (1 UA = 150 millones de
km). Use la fórmula de Simpson con m=4 (cantidad de franjas).
23. Una placa rectangular metálica de 0.45 m por 0.60 m pesa 5 Kg. Se necesita recortar este material para
obtener una placa de forma elíptica, con eje mayor igual a 50 cm, y eje menor igual a 40 cm. Calcule el
área de la elipse y determine el peso que tendrá esta placa. Para calcular el área de la elipse use la fórmula
de Simpson con m = 4. Finalmente, estime cual es el el error de truncamiento en el valor del área calculada.
24. En el techado de las casas se utilizan planchas corrugadas con perfil
ondulado: Cada onda tiene la forma f  x   sen  x  , con un

periodo de 2 pulgadas.
El perfil de la plancha tiene 8 ondas y la longitud L de cada onda se
2
1   f '  x  dx
2
la puede calcular con la siguiente integral: L   0

Este integral no puede ser calculado por métodos analíticos.

a) Use la fórmula de Simpson con m = 4, 6, 8, 10 para calcular L y estime el error en el último resultado.
b) Con el último resultado encuentre la longitud del perfil de la plancha.

25. Una curva C puede darse en forma paramétrica con las ecuaciones: x  f t  , y  g t  , t  a, b
Si no hay intersecciones entre f y g, entonces, la longitud del arco de C se puede calcular con la integral:
b
x '  t    y '  t  dt .
2 2
S
a
Use la fórmula de Simpson, m = 4, para calcular la longitud del arco de la curva dada con las siguientes
 
ecuaciones paramétricas x  2 cos  t  , y  3sen  t  , t  0,
 2 

26. Corriente arriba de una presa, el agua ejerce una presión p  z    .g  D  z  donde p  z  es la presión en
N/m2 ejercida a una elevación z metros por encima del fondo fluvial;   103 kg / m3 es la densidad del
agua; es la aceleración debida a la gravedad y D es la elevación en metros de la superficie del agua por
encima del fondo. Notemos que la presión aumenta con la profundidad. Si se omite la presión
atmosférica, la fuerza f puede ser determinada al multiplicar la presión por el área de la cara de la presa.
Esta área se obtiene integrando la función w(z) que proporciona el ancho del río a una altura z como se
muestra en la figura. Debido a que ambas, la presión y el área, varían con la elevación, la fuerza total se
obtiene al calcular:
D
f    . g.w  z  D  z  dz
0

Prof. Carlos Deudor Gomez 5 Métodos Numéricos

 a) Calcula la fuerza total ejercida por el agua mediante la fórmula compuesta de Trapecios para 6
subintervalos.
b) Calcula la fuerza total ejercida por el agua mediante la fórmula compuesta de Trapecios para 12
subintervalos. Para averiguar la anchura del río w a una profundidad que no aparezca en la figura,
utiliza la media de las anchuras adyacentes.
c) Calcula la fuerza total ejercida por el agua mediante la fórmula compuesta de Simpson para 6
subintervalos.
d) Implementa una función .m que aproxime una integral por el método de Simpson siendo los
parámetros de entrada los nodos x e y la evaluación de la función en los nodos.
e) Calcula la fuerza total ejercida por el agua mediante la fórmula compuesta de Simpson para 12
subintervalos. Para averiguar la anchura del río w a una profundidad que no aparezca en la figura,
utiliza la media de las anchuras adyacentes.
27. La fuerza total ejercida por el mástil de un velero se puede expresar como:
z
30 250 z 
f  .e 10 dz
0 6 z

donde z es la distancia vertical a la cubierta.

a) Calcula la fuerza F aplicando el método de Trapecios con 10 subintervalos.
b) Calcula F aplicando el método de Simpson con 10 subintervalos.
c) La fórmula de punto medio compuesta en [a,b] se obtiene dividiendo el intervalo en m subintervalos
y tomando como nodos los puntos medios de los mismos, x1 , x2 , , xm . Se evalúa la integral por punto
medio en cada subintervalo, multiplicando su longitud por el valor de la función en el nodo
correspondiente y se suman los resultados. La fórmula queda I M h  h  f  x1   f  x2   f  xm  
ba
donde h es la longitud de los subintervalos: h  . Modifica tu programa de trapecios para que
m
b
aproxime la integral I   f  x  dx aplicando la el método de punto medio.
a

d) Aplica el programa obtenido en el apartado (c) para calcular la fuerza F tomando 10 subintervalos.

28. Una partícula de masa m que se desplaza por un fluido está sujeta a una resistencia viscosa R, que a su
vez es función de la velocidad v. La relación entre la resistencia R, la velocidad v y el tiempo t está dada
por la ecuación:
v t  m
t du
 
v t0  R u

Supongamos que R  v    v v para determinado fluido, donde R se da en newtons y v se da en metros

por segundo. Si m = 10 kg y si v(0) = 10 m/seg, aproxima el tiempo que la partícula tarda en reducir su
velocidad a v = 5 m/s, utilizando para ello:

Prof. Carlos Deudor Gomez 6 Métodos Numéricos

 a) El método de trapecios con 10 subintervalos.
b) La función trapz.m de manera que el resultado coincida con el del apartado i).
c) El método de Simpson con 10 subintervalos.
d) El método de trapecios iterativo hasta alcanzar una tolerancia 10-6.
e) Compara los resultados obtenidos en los apartados anteriores con la solución exacta.
29. La función de densidad de una distribución normal de probabilidad es
2
1  x2
f x  e
2
b
La probabilidad de que x esté en el intervalo [a, b] es  f  x  dx
a

a) Calcula por el método de trapecios la probabilidad de que x esté en [-1,1].tomando 10 subintervalos.

b) Obtén la probabilidad de que x esté en [-1,1] utilizando el método de Simpson-Richardson iterativo
partiendo de 2 subintervalos hasta alcanzar una precisión de 10-6. ¿Qué número de subintervalos se
han tomado para alcanzar la precisión indicada?
c) Crea una función .m de manera que llamando a simpson .m y dado un valor a obtenga como salida la
probabilidad de que x esté en [-a, a]. Completar la siguiente tabla:

d) Modifica la función anterior para aplicarle el método de bisección paso a paso y obtener el valor de a
tal que la probabilidad de que x esté en [-a, a] sea del 80%. ¿Cuántas iteraciones se realizan para
alcanzar una tolerancia de 10-3 ?
30. Ejercicio
a) Calcula el área de la parte del círculo centrado en el origen y de radio r, situada en el primer cuadrante,
por el método de trapecios iterativos, con determinada tolerancia. Compara el resultado con el valor
exacto, πr2/4
b) Repite el cálculo con el sector de la figura, que es el mismo cuadrante girado 45º.
Compara el resultado número de iteraciones con los del apartado anterior.
Explica por qué se produce la diferencia.

El volumen del cuerpo de revolución engendrado por la superficie limitada entre la curva, el eje OX
y las ordenadas x = a y x = b, al girar entorno a dicho eje viene dado por la integral
b
V    f  x  dx
2
a

c) Halla el volumen de la semiesfera engendrada al girar el cuadrante de círculo del apartado a) entorno
al eje OX (Sol: 2/3πr3). Explica por qué aquí no surgen las dificultades encontradas en el cálculo del
área.

Referencia Bibliográfica

Prof. Carlos Deudor Gomez 7 Métodos Numéricos

 N° Autores Titulo Páginas
[1] Métodos numéricos para ingenieros Steven C. Chapra 965
[2] Análisis Numérico Richard L. Burden y J Douglas Faires 831
[3] Análisis Numérico Gutierrez, Olmos, Casillas 421
[4] Análisis Numérico Axel Osses A. 163
[5] http://www.robcas64.com/Numerico/Numerico1.html
[6] http://ocw.upm.es/matematica-aplicada/programacion-y-metodos-numericos

Prof. Carlos Deudor Gomez 8 Métodos Numéricos