Calculo II

Sumatorias, integrales definidas, cálculo

de áreas, métodos numéricos,
sumatorias e integrales dobles,
Interpolación.

Universidad Nacional Del Callao

 Calculo II 2012

Índice
1) Sumatorias 2
 Definición
 Propiedades
 Fórmulas de sumatoria
 Ejemplo
2) Calculo de la integral definida mediante sumatorias 9
 Definición
 Propiedades de la integral definida
 Ejemplos

3) Calculo de Áreas por sumatorias 17

 Introducción
 Partición de un intervalo
 Refinamiento de una partición
 Área de una región plana por sumatorias de rectángulos

4) Métodos numéricos para calcular integrales definidas 35

 Método del rectángulo 35
 Método del trapecio 37
 Método de Simpson 40
 Integración por series de Taylor 51

5) Sumatorias dobles 58
 Definición
 Integral doble sobre un rectángulo
 Integral doble como volumen
 Propiedades
 Ejercicios

6) Calculo de integrales dobles por sumatorias dobles 63

 Definición
 Integral doble sobre un rectángulo
 Integral doble como volumen
 Propiedades
 Ejercicios

7) Calculo de integrales por interpolación 72

 Introducción
 Interpolación de Lagrange
 Interpolación de Hermite
 Interpolación de trazadores cúbicos

1
 Calculo II 2012

SUMATORIAS

1.-DEFINICION
Consideremos “m” y “n” dos números naturales tal que m ≤ n, y “f” una función definida para cada
n
i Є N donde m ≤ i ≤ n, luego la notación  f  i  , nos representa la suma de los términos; f(m);
i m

f(m+1); f(m+2); …; f(n); es decir:

n

 f  i   f  m   f  m  1  f  m  2  ...  f  n 
i m

Donde “i” es el índice o variable, “m” es el límite inferior y “n” el límite superior.
Particularmente: a la suma de los n números a1, a2, …, an , representaremos por la notación:
n

ai 1
i  a1  a2  ...  an
n
En la sumatoria  f  i  , existen (n-m+1) términos los cuales son f(m), f(m+1), f(m+2), …, f(m+(n-
i m

m))

2.-PROPIEDADES

Siendo f,g funciones definidas  i Є Z, k contante.

n
a)  k  kn
i 1

DEMOSTRACION:

 k  k  k  ...  k (" n " veces)  kn

i 1

n
b)  k  (n  m  1)k
i m

DEMOSTRACION:
n

 k  k  k  ...  k (" n  m  1" veces)  (n  m  1)k

i 1

n n
c)  kf  i   k  f  i 
i 1 i 1

2
 Calculo II 2012

DEMOSTRACION:

n n

 k f i   kf 1  kf  2  ...  kf  n   k ( f 1  f  2  ... f  n )  k  f i 

i 1 i 1
n n n
d)  ( f i   g i )   f i    g i 
i 1 i 1 i 1

DEMOSTRACION:

 ( f  i   g i )  ( f 1  g 1)  ( f  2   g  2 )...  ( f  n   g  n )

i 1
n n
 ( f 1  f  2   ... f  n )  ( g 1  g  2   ...g  n )   f  i    g  i 
i 1 i 1

b bc
e)  f i    f i  c 
i a i a c

DEMOSTRACION:

 f  i   f  a   f  a  1  f  a  2   ...  f b 
i a

 f  (a  c)  c   f  (a  c  1)  c   f  (a  c  2)  c   ...  f  (b  c)  c 
bc
  f i  c 
i a c

b b c
f)  f i    f i  c 
i a i a c

DEMOSTRACION:

 f  i   f  a   f  a  1  f  a  2   ...  f b 
i a

 f  (a  c)  c   f  (a  c  1)  c   f  (a  c  2)  c   ...  f  (b  c)  c 
b c
  f i  c 
i a c

n
g)  ( f  i   f  i  1)  f  n   f  0 1ra Regla Telescópica
i 1

3
 Calculo II 2012

DEMOSTRACION:
n n n

 ( f  i   f  i  1)   f i    f i  1
i 1 i 1 i 1

 ( f 1  f  2   ...  f  n  1  f  n )  ( f  0   f 1  ...  f  n  2   f  n  1)

 f  n   f  0
n
h)  ( f i   f i  1)  f  n   f  k  1 1ra Regla Telescópica Generalizada
i k

DEMOSTRACION:

n n n

 ( f  i   f  i  1)   f  i    f  i  1
i k i k i k

 ( f  k   f  k  1  ...  f  n  1  f  n )  ( f  k  1  f  k   ...  f  n  2   f  n  1)

 f  n   f  k  1

n
i)  ( f  i  1  f  i  1)  f  n  1  f  n   f 1  f  0 2da Regla Telescópica
i 1

DEMOSTRACION:

n n n

 ( f  i  1  f  i  1)   f i  1   f i  1
i 1 i 1 i 1

 ( f  2   f  3  ...  f  n   f  n  1)  ( f  0   f 1  ...  f  n  2   f  n  1)

 f  n  1  f  n   f 1  f  0 

n
j)  ( f  i  1  f i  1)  f  n  1  f  n   f  k   f  k  1 2da Regla Telescópica
i k

Generalizada

DEMOSTRACION:

n n n

 ( f  i  1  f  i  1)   f i  1   f i  1
i k i k i k

 ( f  k  1  f  k  2   ...  f  n   f  n  1)  ( f  k  1  f  k   ...  f  n  2   f  n  1)

 f  n  1  f  n   f  k   f  k  1

4
 Calculo II 2012

3.-FORMULAS DE LA SUMATORIA
n
n(n  1)
a) i 
i 1 2
DEMOSTRACION:

Aplicando la regla telescópica

  i  1 
n
 i 2  f  n   f  0  , donde f (i)  (i  1)2
2

i 1

   i  1 
n
 i 2   n  1  1
2 2

i 1

 i  2i  1  i 2   n 2  2n
n
2

i 1
n n
2 i   1  n 2  2 n
i 1 i 1
n
2 i  n  n 2  2 n
i 1
n
2 i  n 2  n
i 1
n
n(n  1)
i 
i 1 2

n
n  n  1 2n  1
b) i
i 1
2

6
DEMOSTRACION:

  i  1 
n
 i 3  f  n   f  0  , donde f (i)  (i  1)3
3

i 1

   i  1 
n
 i 3   n  1  1
3 3

i 1

 i  3i 2  3i  1  i 3    n  1  1
n
3 3

i 1
n n n
3 i 2  3 i  1   n  1  1
3

i 1 i 1 i 1
n
3
3 i 2  n  n  1  n   n  1  1
3

i 1 2
n
3
3 i 2   n  1   n  1  n  n  1 
3

i 1 2
n
n(n  1)  2n  1
 i 
i 1 6

5
 Calculo II 2012

n 2  n  1
2
n
c)  i  3

i 1 4
DEMOSTRACION:

  i  1 
n
 i 4  f  n   f  0  , donde f (i)  (i  1)4
4

i 1

   i  1 
n
 i 4   n  1  1
4 4

i 1

 i  4i 3  6i 2  4i  1  i 4    n  1  1
n
4 4

i 1
n n n n
4 i 3  6 i 2  4 i  1   n  1  1
4

i 1 i 1 i 1 i 1
n
6 4
4 i 3  n  n  1 2n  1  n(n  1)  n   n  1  1
4

i 1 6 2
n
3 i 3  n  n  1 2n  1  2n(n  1)  n   n  1  1
4

i 1
n
3 i 3   n  1  (n  1)  2n(n  1)  n  n  1 2n  1
4

i 1

n 2  n  1
2
n

 i3 
i 1 4

n n  n  1  6n3  9n 2  n  1
d)  i4 
i 1 30
DEMOSTRACION:

  i  1 
n
 i 5  f  n   f  0  , donde f (i)  (i  1)5
5

i 1

   i  1 
n
 i 5   n  1  1
5 5

i 1

 i  5i 4  10i 3  10i 2  5i  1  i 5    n  1  1
n
5 5

i 1
n n n n n
5 i 4  10 i 3  10 i 2  5 i  1   n  1  1
5

i 1 i 1 i 1 i 1 i 1
n
10 2 10 5
5 i 4  n (n  1) 2  n  n  1 2n  1  n(n  1)  n   n  1  1
5

i 1 4 6 2
n
5 5 5
5 i 4  n 2 (n  1) 2  n  n  1 2n  1  n(n  1)  n   n  1  1
5

i 1 2 3 2
n
5 5 5
5 i 4   n  1  n(n  1)  n  n  1 2n  1  n 2 (n  1) 2
5

i 1 2 3 2
n
3 n2 n 6n 3  9n 2  n  1
5 i 4   n  1 [n 4  n3   ]  n  n  1 [ ]
i 1 2 6 6 6
n n  n  1  6n3  9n 2  n  1
i
i 1
4

30

6
 Calculo II 2012

n
r n 1
e) Si r≠1,  ar i  a.
i 1 r 1
DEMOSTRACION:

 ar
i 1
i
 a(r 0  r  r 2  ...  r n )  a(1  r  r 2  ...  r n )

1 rn r n 1
 a.  a.
1 r r 1

n
f) Si r=1,  ar
i 1
i
 n.a

Cuyo resultado se encuentra demostrado en la propiedad 1.

4.-EJEMPLOS

a) Calcular la suma de los n primeros términos de: 1.6+2.7+3.8+4.9+…

Solución:

Observando los términos de la suma nos damos cuenta que el término general es ai=i(i+5);
por lo tanto, la suma se expresa como:

 i  i  5    i 2  5i 
n n

i 1 i 1
n n
  i 2  5 i
i 1 i 1

n  n  1 2n  1  n  n  1 
  5 
6  2 
n3  9n 2  8n

3

n
1
b) Calcular   2i  1 2i  1
i 1

Solución:

Expresando en fracciones parciales:

1 A B A(2i  1)  B(2i  1)
  
 2i  1 2i  1 2i  1 2i  1 (2i  1)(2i  1)
1 2 A  A  2 Bi  B

 2i  1 2i  1 (2i  1)(2i  1)

7
 Calculo II 2012

1  i  2 A  2B   ( A  B) (por identidad de polinomios)

2 A  2B  0
A B 1
1 1
Entonces, A  ,B  
2 2
Reemplazando en la sumatoria:
1 1
n n n
2  2 ) 1 ( 1  1 )1 ( 1  1 )
 (
i 1 2i  1

2i  1 2 i 1 2i  1 2i  1

2 i 1 2i  1 2i  1

Por la 1ra Regla Telescópica:

1 1
   f  n   f  0   , donde f  i  
2 2i  1
1 1  1 1  1 2n n
   1  1    ( )
2  2n  1  2  2 n  1  2 2 n  1 2 n  1

n
1
c) Calcular  log  2  .log  2 
i 1
2i 2i  2
a a

Solución:
n n
1 1
 log  2  .log  2    2i.log  2  .  2i  2  .log  2 
i 1
2i 2i  2
i 1
a a a a

n n
1 1 1 1
  2i  2i  2    i  i  1
 log a 2   2 log a 2 
2 2
i 1 i 1

1 n
1 1 1  1  1 n
  ( i  i  1)  1 
 n  1   .
 log a 4   log a 4   log a 4  n  1
2 2 2
i 1

n
d) Calcular  sen
i 1
2i
ax

Solución:
n

 sen
i 1
2i
ax  sen 2 ax  sen 4 ax  sen6 ax  ...  sen 2 n ax

 sen 2 ax 1  sen 2 ax  sen 4 ax  ...  sen 2 n 2 ax 

 1  sen 2 n ax 
 sen ax 
2
  tg ax(1  sen ax)
2 2n

 1  sen ax 
2

8
 Calculo II 2012

CALCULO DE LA INTEGRAL DEFINIDA MEDIANTE SUMATORIA

Sea D el conjunto de todas las particiones posibles  del intervalo a, b . Si f es una función
acotada sobre a, b entonces existen números m y M tal que:
m  f ( x)  M , x   a, b
Se sabe que la siguiente desigualdad se cumple
m  b  a   L( f , P)  U ( f , P)  M (b  a)
Para toda partición P en D, esto asegura que el conjunto numérico L( f , P) / P  D es acotado

superiormente y el conjunto U ( f , P) / P  D es acotado inferiormente, luego el conjunto

L( f , P) / P  D tiene un supremo (la mayor cota inferior) y U ( f , P) / P  D tiene un ínfimo

(mínima cota superior) con estos valores supremo e ínfimo daremos la siguiente definición:

1. DEFINICION:

Si f es una función acotada en a, b , al número sup L( f , P) / P  D se llama integral inferior
de f en a, b y se indica.
b

 f ( x)dx  sup L( f , P) / P  D = integral inferior de f desde a hasta b.

a

Al número inferior U ( f , P) / P  D se llama integral superior de f en a, b y se indica

b

 f ( x)dx  inf U ( f , P) / P  D = integral superior de f desde a hasta b.

a

En muchas aplicaciones del Cálculo integral es preferible definir la integración como una
sumatoria. De hecho, el Cálculo integral se inventó con el fin de calcular el área de las superficies
limitadas por curvas, suponiéndose la superficie dividida en “un número infinito de partes
infinitamente pequeñas que se llama elementos, siendo la suma de las áreas de todos estos
elementos el área buscada”. Históricamente, el signo integral no es otra cosa que la “S” larga,
empleada por los primeros autores para indicar la palabra suma.
Esta nueva definición, que se desarrollara en el presente informe, es de importancia fundamental,
y es indispensable comprender a fondo lo que se quiere decir, para que podamos aplicar el Cálculo
integral a los problemas prácticos.

A. Teorema fundamental del Cálculo integral:

Si ( x) es la derivada de f ( x) , se sabe que el valor de la integral definida
b

 ( x)dx  f (b)  f (a)

a
………. (1)

9
 Calculo II 2012

Dada el área de la superficie limitada por la curva y  ( x) , el eje de las x y las ordenadas
correspondientes a x  a y x  b .
Ahora bien, a propósito de esta área hagamos la siguiente construcción. Dividamos el
intervalos desde x  a hasta x  b en un número “n” de partes iguales, tracemos las
ordenadas en los puntos de división y completemos los rectángulos trazando líneas
horizontales por las extremidades de las ordenadas, como se indica en la figura (1). Es
evidente que la suma de las áreas de estos “n” rectángulos (el área sombreada) es un valor
aproximado del área que consideramos. Además es también evidente que el límite de la
suma de las áreas de estos rectángulos, cuando se aumenta indefinidamente su número
“n”, será igual al área bajo la curva.
Ahora efectuamos la siguiente construcción más general. Dividimos el intervalo en “n”
partes, que no serán necesariamente iguales, y levantemos ordenadas en los puntos de
división. Elijamos de cualquier modo un punto dentro de cada parte,

Figura (1) Figura (2)

Levantemos ordenadas en estos puntos y tracemos por sus extremidades líneas
horizontales de manera de formar rectángulos, como se indica en la figura (2). Entonces,
como antes, la suma de las áreas de estos “n” rectángulos (el área sombreada) es igual,
aproximadamente, al área bajo la curva; y el límite de esta suma cuando “n” tiende al
infinito y cada parte tiende a cero es, precisamente, el área bajo la curva. De estas
consideraciones vemos que la integral definida (1) puede mirarse como el límite de una
suma. Ahora formulemos este resultado.

Figura (3)
 Designemos las longitudes de las divisiones sucesivas por

x1 , x2 , x3 ,..., xn .

 Designemos las abscisas de los puntos elegidos en cada una de las divisiones por

10
 Calculo II 2012

x1 , x2 , x3 ,..., xn .
Entonces las ordenadas de las curvas en esos puntos son
( x1 ), ( x2 ), ( x3 ),..., ( xn ).
 Las áreas de los rectángulos sucesivos son, evidentemente,
( x1 )x1 , ( x2 )x2 , ( x3 )x3 ,..., ( xn )xn .
 Por tanto, el área bajo la curva es igual a
lim[( x1 )x1  ( x2 )x2  ( x3 )x3  ...  ( xn )xn ].
n
b
Pero según (1) el área bajo la curva = a
( x)dx.
Luego,
b
(A)  a
( x)dx  lim[( x1 )x1 ( x2 )x2  ( x3 )x3  ...  ( xn )xn ].
n 
Esta igualdad se ha deducido sirviéndonos de la noción de área.
La intuición nos ha ayudado en establecer el resultado. Ahora consideremos (A)
simplemente como un teorema de Análisis matemático, que se puede formular como
sigue:
Teorema Fundamental del Cálculo integral.- Sea ( x) una función continua en el intervalo
desde x  a hasta x  b . Dividiéndose este intervalo en “n” subintervalos cuyas
longitudes son x1 , x2 ,..., xn , y elijándose puntos, uno en cada subintervalo, que
tengan las abscisas x1, x2,..., xn, respectivamente. Considere la suma
n

(2) ( x1 )x1  ( x2 )x2  ...  ( xn )xn   (xi )xi .

i 1
Entonces, el valor límite de esta suma cuando “n” tiende al infinito, y cada subintervalo
tiende a cero, es igual al valor de la integral definida
b
 a
( x)dx.
La igualdad (A) se puede abreviarse como sigue:
n
( x)dx  lim  ( xi )xi .
b
(3) 
a n 
i 1
La importancia de este teorema resulta del hecho que si podemos calcular, por
integración, una magnitud que sea el límite de una suma de la forma (2).
Puede observarse que cada término de la suma (2) es una expresión diferencial, puesto
que las longitudes x1 , x2 ,..., xn , tienden a cero. Además, cada termino se llama en
elemento de la magnitud que se trata de calcular.
La siguiente regla será muy útil en la aplicación de este teorema a los problemas prácticos.

11
 Calculo II 2012

B. Regla para Aplicar el teorema Fundamental:

PRIMER PASO: se divide la magnitud buscada en partes semejantes de manera que sea
claro que el resultado deseado se encuentra tomando el límite de una suma de esas
partes.
SEGUNDO PASO: para las magnitudes de estas partes se hallan expresiones tales que su
suma sea de la forma (2).
TERCER PASO: elegidos los limites apropiados, x  a y x  b , se aplica el teorema
fundamental.
n
lim  ( xi )xi   ( x)dx
b

n  a y se integra.
i 1

C. Propiedades de la Integral Definida

 Si ‘f’ y ‘g’ son integrables en a, b y ‘c’ es una constante, entonces:

b b

 cf ( x)dx  c  f ( x)dx
a a
b b b

 [ f ( x)  g ( x)]dx   f ( x)dx   g ( x)dx

a a a

 Si x está definido para x=a, entonces:

a

 f ( x)dx  0
a

 Si f es integrable en a, b , entonces:

b a

 f ( x)dx   f ( x)dx
a b

 Propiedad de aditividad del valor: si f es integrable en los dos intervalos cerrados definidos
por a, b y c entonces
c b c

 f ( x)dx   f ( x)dx   f ( x)dx

a a b

 Si en el intervalo [a,b] la función f es mayor o igual que la función g, entonces:

b b

 f ( x)dx   g ( x)dx
a a

 Si f ( x)  0 , x  a, b ,entonces:
b

 f ( x)dx  0
a

12
 Calculo II 2012

 Si 0  f ( x), x   a, b 
b
0   f ( x)dx
a

 Si en el intervalo a, b la función ‘f’ es mayor que la función ‘g’, entonces:

b b

 f ( x)dx   g ( x)dx
a a

 Si f ( x)  0, x   a, b  , entonces:
b

 f ( x)dx  0
a

 Si 0  f ( x), x   a, b  , entonces:

b
0   f ( x)dx
a

 Si f es una función par sobre a, b , f ( x)  f ( x), x   a, b , y a, b es un intervalo

simétrico respecto del origen, entonces:
b b

 f ( x)dx  2 f ( x)dx
a 0

 Si f es una función impar sobre a, b , f ( x)   f ( x), x   a, b ,y a, b es un intervalo

simétrico respecto al origen, entonces:
b

 f ( x)dx  0
a
b
bk
 
a
f ( x)dx  
ak
f ( x  k )dx

 Si f es una función continua en el intervalo , a entonces:

a a

  f ( x)dx   f (a  x)dx
0 0

13
 Calculo II 2012

Ejemplo 1:
1
Halle: 
2
x3dx
que se muestra en la figura4

Figura4
Como f ( x)  x es continua en el intervalo [-2,1] sabemos que es integrable.
3

b  a 1  ( 3
Dividimos el intervalo en ‘n’ subintervalos de igual longitud xi   
n n n
Y para el cálculo de la integral consideremos el extremo derecho de cada intervalo
3i
ti  2 
n
n n
3i 3 3 3 n 36i 54i 2 27i 3
 f (ti )xi  lim  (2  )  lim  (8 
1
2 x dx  lim  2  3 )
3
n  n  n n n n i 1 n n n
i 1 i 1

1 3 36n(n  1) 54n(n  1)(2n  1) 27n2 (n  1)2

2     
3
x dx lim [ 8n ]
n  n 2n 6n 2 4n 3
1 54 81 27 81 81 81 81 15

2
x3dx  lim[   
n 
  2 
n n n
 2 ]  24   
4 2n 4n 4 4
1 15

2
x3dx  
4
Observación: Esta integral definida es negativa, no representa el área graficada. Las integrales
definidas pueden ser positivas, negativas o nulas.

Ejemplo 2:
Hallar el área limitada (A) limitada por la gráfica f(x)=x+1 en el intervalo [0,4].
Solución:
Dado que a=0 y b=4, se tiene que la longitud de cada subintervalo es:
40 4
x  
n n
Sustituyendo en la siguiente formula

14
 Calculo II 2012

n
ba ba
A  lim  f (a  k ).
n  n n
k 1

Se tiene:
n
4 4 4 n 4
A  lim  f (0  k ).  lim  f ( )
n  n n n  n n
k 1 k 1

4 n 4k n
4 4 n n

 lim
n  n

k 1
f (  1)     k  1
n k 1 n  k 1
n k 1 

4  4  n(n  1)   n 16  n(n  1)  

 lim     n       4
  2   k 1  2  n  
n  n n 2

  1   1
 lim 8 1    4  8lim 1    4  8  4  12
n 
  n  n 
 n
Por lo tanto, A=12u²

Figura5
Grafica realizada en Logger Pro 3.2.1 Español (figura5)

Ejemplo 3:
2

 x  x  3dx
3
Hallar:
1

Como f ( x)  x  x  3 es continua en el intervalo [-1,2] sabemos que es integrable.

3

b  a   ( 3
Dividimos el intervalo en ‘n’ subintervalos de igual longitud xi   
n n n
Y para el cálculo de la integral consideremos el extremo derecho de cada intervalo
3i
ti  
n
  3i   3 3 n 27i 3 27
 x  x  3dx  lim  f (ti )xi  lim  ( )3   1    3  lim  ( 3  2 i 2  i  1)
n n
2
3 3i 12
 n 
i 1
n 
i 1  n  n   n n n i 1 n n n

3 27n2 (n  1)2 12n(n  1) 27n(n  1)(2n  1)

   n n [ 4n3  2n 
2
x 3
 x  3 dx  lim  n]
6n 2
27  1  9  1  1  1
 x  x  3 dx  lim3[
2 57
1    1   2    6 1    1] 
3
 n  4  n  2  n  n  n 4

15
 Calculo II 2012

Ejemplo 4:

x , x  3
2

Calcular el área de la región sombreada acotada por la curva, f ( x)   el eje X y


 6 x  x 2
, x  3
las rectas x=1, x=7. De la figura6

Figura6

Calcularemos el R:
3 1 2
x  
n n
2i
ti   
n
2
 2   2  4 4
 f 1  i   1  i   1  i  2 i 2
 n   n  n n
n
 4 4 2
 lim  1  i  2 i 2 
i 1  n n
n  n
2n 8n(n  1) 8n(n  1)(2n  1)
 lim[   ]
n  n 2n 2 6n 3
 1  4  1  1  26
 lim[2  4 1    1   2  ] 
n 
 n  3  n  n 3
Calcularemos el R1:
63 3
x  
n n
3i
ti  3 
n
2
 3   2   2  9
 f  3  i   6 1  i   1  i   9  2 i 2
 n   n   n  n
entonces se tiene:


n
9 23 n
 1 i2 
 lim   9  2 i  27 lim    3 
i 1  n
n  n 
n i 1  n n 

16
 Calculo II 2012

1 n(n  1)(2n  1) 1 1  1
 lim[1  ]  lim[1  (1  )  2  ]  18
n  3 6 n  6 n  n
80
Entonces el área de la región sombreada será R+R1 lo que igual a
3

Ejemplo 5:
1
Hallar  senh( x)dx
0

Como f ( x)  senh( x) es continua en el intervalo [0,1] sabemos que es integrable.

Dividimos el intervalo en ‘n’ subintervalos de igual longitud xi  b  a  1    1
n n n
Y para el cálculo de la integral consideremos el extremo derecho de cada intervalo
i
ti 
n
i
 senh( )
n
n
i 1
 lim  senh( )
n  n n
i 1

 1 1 1 
 cosh(n  1) n  cosh(n n )  cosh n  1  1
 lim  
n 
 1 n
2senh( )
 n 

 
 
 
1 1
 cosh(1  )  cosh(1)  cosh  1 
 lim  n n   2 cosh(1)  2  cosh(1)  1
n   1  2
 senh( ) 
 2 n 
 1 
 n 

17
 Calculo II 2012

CÁLCULO DE ÁREAS USANDO SUMATORIAS

1.-INTRODUCCIÓN

En este capítulo desarrollaremos el cálculo de áreas de regiones planas usando aproximaciones

mediante la sumatoria de áreas de rectángulos, para lo cual es necesario hacer uso de algunos
conocimientos previos al tema, los cuales se desarrollará a continuación.

2.-Partición de un intervalo
Dado un intervalo cerrado  a, b de ℝ, se llama partición de  a, b a cualquier conjunto finito

P   x0 , x1 ,..., xn  de puntos de  a, b tal qué a  x0  x1  ...  xn  b .Se denotará por P  a, b

al conjunto de todas las particiones de  a, b .

Entonces podemos decir que  xk 1 , xk  es, para k  1, 2, . . . , n , el sub intervalo k-ésimo de la

partición y llamamos amplitud del sub intervalo a :
xk  xk  xk 1
(Propiedad 1)
Se llama diámetro de la partición a la mayor de las amplitudes, es decir:

P  max x1 , x2 ,..., xn  , k  1, 2, . . . , n .

(Propiedad 2)

3.-Refinamiento de una partición:

Dado una partición de  a, b esta se puede hacer más fina usando puntos. Si P1 se obtiene de P
añadiendo por lo menos un punto, entonces P  P1 .
Al número P  max xi  xi 1 / i  1,2,..., n le llamamos norma o diámetro de3 la partición P y
que es la mayor de las longitudes  i x (Propiedad 3).

Ejemplo 1:
Para el intervalo 0, 4 los siguientes conjuntos representan particiones:
P1  0,1, 2,3, 4
P2  0,1,3 / 2, 4
P3  0, 2,5 / 2,3,7 / 2, 4
P4  0, 2 / 3,1,3 / 2, 2,5 / 2, 4
Siendo P4 un refinamiento de P1 .

Ejemplo 2:
Dado un intervalo  0, 4 y la partición P  0,2,5 / 2,3,7 / 2,4 hallar la norma de la partición P .

18
 Calculo II 2012

Calculamos las longitudes de  i x :

1 x  2  0  2
5 1
2 x  2
2 2
5 1
3 x  3  
2 2
7 1
4 x   3 
2 2
7 1
5 x  4  
2 2
Luego, nos damos cuenta que la norma de la partición P es: P  2 .

4.-Área de una región plana por sumatoria de rectángulos

“Si en cualquier figura delimitada por rectas y por una curva se inscriben y circunscriben
rectángulos en número arbitrario, y si la anchura de tales rectángulos se va disminuyendo a la par
que se aumenta su número hasta el infinito, afirmo que las razones entre las figuras inscrita y
circunscrita y la figura curvilínea acabarán siendo razones de igualdad”´(Isaac Newton).
i. Sumas inferiores y sumas superiores
Consideremos una región plana acotada por la gráfica de una función y  f  x  continua y no
negativa (ver fig. 1) la región está limitada por el eje x y las rectas x  a, x  b .

Fig. 1
Para aproximar su área comenzamos dividiendo el intervalo  a, b en sub intervalos, cada uno de
ba
longitud x  (Propiedad 4),(ver fig. 2), sus puntos terminales son:
n
a   x   a  1 x   a  2  x   ...  a  n  x 
(Propiedad 5)

19
 Calculo II 2012

Fig. 2

Por ser f continua, el teorema de los valores extremos verifica asegura la existencia de un mínimo
y de un máximo de f  x  en cada sub intervalo.

f  m  =valor mínimo de f  x  en el i-esimo sub intervalo.

f  M  = valor máximo de f  x  en el i-esimo sub intervalo.

Ejemplo 3:
Calcular las sumas inferiores y superiores para la región limitada por la gráfica de f  x   x y el
2

eje  .
Solución:
a) El primer paso es realizar la partición del intervalo  0, 2 en sub intervalos, cada uno de
longitud.
ba 20 2
x   
n n n
Como f es creciente en el intervalo  0, 2 , en cada sub intervalo el valor mínimo de f
ocurre en el punto terminal de la izquierda y el máximo en el punto terminal de la
derecha.
b) El segundo paso es calcular los puntos terminales en ambos lados:
 2  2  i  1
Puntos terminales izquierdos: mi  0   i  1   
n n
 2  2i
Puntos terminales derechos: M i  0   i    
n n
c) El tercer paso es usar los valores de ambos lados.
*Usando los de la izquierda, tenemos la siguiente suma inferior:

20
 Calculo II 2012

n n
 2  i  1   2 
s  n    f  mi x   f   
i 1 i 1  n  n 
 2  i  1   2 
2
n
    
i 1  n  n
2
8
   3   i 2  2i  1
n

i 1  n 

8 n 2 n n

3   
 i  2 i  1
n  i 1 i 1 i 1 

8  n  n  1 2n  1  n  n  1   
 3 
2   n
n  6  2   

3 
2n3  3n2  n 
4

3n
8 4 4
   2
3 n 3n

* Usando los de la derecha, tenemos la siguiente suma superior:

n n
 2i   2 
S  n    f  M i x   f    
i 1 i 1  n  n 
2
n
 2i   2 
    
i 1  n   n 
2
n
8
   3  i2
i 1  n 

8  n  n  1 2n  1 
  
n3  6 

3 
2n3  3n2  n 
4

3n
8 4 4
   2
3 n 3n
La conclusión que se puede rescatar de este ejemplo es que se observa claramente algunas
propiedades de las sumas inferiores y superiores.
a. Para cualquier valor de n la suma inferior es menor o igual que la suma superior.
8 4 4 8 4 4
s  n    2    2  S  n
3 n 3n 3 n 3n
(Propiedad 6)

21
 Calculo II 2012

b. La diferencia entre estas disminuye conforme n crece. De hecho, si tomamos el límite

8
n   , tanto las sumas inferiores como las sumas superiores tienden a .
3
8 4 4  8
lim s  n   lim    2   Límite de las sumas inferiores
n 

n  3 n 3n  3
8 4 4  8
lim S  n   lim    2   Límite de las sumas superiores
n 

n  3 n 3n  3

ii. Límite de las sumas inferiores y superiores

Sea f contínua y no negativa en un intervalo  a, b . Los límites cuando n   , de las sumas
inferiores y de las sumas superiores son idénticos .Es decir:
n
lim s  n   lim  f  mi  x
n n
i 1
n
 lim  f  M i  x
n
i 1
n
 lim  S  n 
n 
i 1

(Propiedad 7)
ba
Con x    donde f  mi  , f  M i  son los valores mínimo y máximo de f en el
 n 
i  ésimo sub intervalo.

Ahora supongamos que estamos interesados en calcular el área que esta debajo de una curva
y  f  x  en el intervalo cerrado  a, b (para facilitar el ejemplo, consideremos que f es
positiva). Una forma sencilla de aproximar dicha área es dividir el intervalo  a, b en pequeños sub
intervalos y sumar las áreas de los rectángulos que tienen como base los sub intervalos y como
altura tienen el valor de la función f en un punto de dicho sub intervalo.
Nota: cuanto más pequeña sea la base de los rectángulos mejor será la aproximación.

iii. Definición de áreas mediante rectángulos

Sea f una función continua no negativa en el intervalo  a, b .El área de la región limitada por la
gráfica de f ,el eje x  a y x  b es la siguiente:
n
xi 1  ci  xi A  lim  f  ci x
n 
i 1

22
 Calculo II 2012

xi 1  ci  xi
(Propiedad 8)
ba
Donde x    (ver fig. 3).
 n 

Fig. 3
 La aproximación para calcular el área de este tipo de regiones se puede hallar usando una
serie de rectángulos inscritos.
Si A es el área de la región pedida, se cumple que :
n
A   f  ci i x Aproximación por defecto
i 1
n
A  lim  f  ci i x Valor exacto
n 
i 1

 En forma similar se puede aproximar el área por exceso , usando la serie de rectángulos
circunscritos.
n
A   f  ci i x Aproximación por exceso (Propiedad 8)
i 1
n
A  lim  f  ci i x Valor exceso (Propiedad 9)
n 
i 1

5.-EJEMPLOS:
3
1. Hallar el área de la región A acotada por la gráfica de y  x  4 al eje X , y a las rectas
2
x 0,x  4

x  4  Donde x   0, 4
3
y  f ( x) 
2
Primero calculamos la partición del intervalo  0, 4
40 4 4
x    x 
n n n

23
 y

12 Calculo II 2012

10

A
2

x
-7 -6 -5 -4 -3 -2 -1 1 2 3 4 5 6 7 8 9 10 11

-2

f(x)=3x/4+4
Relleno 1 -4
x=4 Luego calculamos el punto terminal (ci )
4 4 -6
ci  a  ix ci  0  i  i
n n
Como:
3 3 4i 6
f ( x)  x4 f (ci )  ( )  4  f (ci )  i  4
2 2 n n

Ahora procedemos a hallar el área de la región acotada

n
Area  lim  f (ci )x
n 
i 1

6 4 n n
 24 16 
A  lim   i  4   lim   2 i  
i 1  n  n n i 1  n n
n 

 24 n n
16   24n  n  1 16n 
A  lim  2  i     lim  2
 
 i 1
n  n
i 1 n 
n 
 2n n 
  1 
A  lim 12 1    16 
n 
  n 
A  28u 2

2. Hallar el área de la región A comprendida entre y  x 2 , y  4  3x 2

6 y

a
2

x
-4.5 -4 -3.5 -3 -2.5 -2 -1.5 -1 -0.5 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5

-1

-2

-3

f(x)=x^2
f(x)=4-3x^2
Relleno 1
24
 Calculo II 2012

y  f ( x)  (4  3x2 )  x2  4  4 x2  Donde x   1,1

1  (1) 2
Primero calculamos: x  
n n
2
Ahora: ci  a  ix  ci  1  i
n
2 16 2 16
Y como: f ( x)  4  4 x 2  f (ci )  4  4( i  1) 2   i  i
n n2 n
Entonces recién ahora procedemos a hallar el área de la región acotada
n
Area  lim  f (ci )x
n 
i 1
n
16 2 16 2
A  lim  ( i  i)
n 
i 1 n2 n n
n
32 2 32
A  lim  ( i  2 i)
n 
i 1 n3 n
 n 32 n
32 
A  lim   3 i 2   2 i 
n 
 i 1 n i 1 n 
 32n(n  1)(2n  1) 32n( n  1) 
A  lim   
n 
 6n 3 2n 2 
 16 1 1 1 
A  lim   (1  )(2  )  16(1  ) 
n 
 3 n n n 
16
A  u2
3
3. Hallar el área de la región A acotada por la gráfica de y  6 x  x 2  x3 al eje  y a las
rectas x  1, x  3

10 y

a
3

1
x
-6 -5 -4 -3 -2 -1 1 2 3 4 5 6 7 8 9

-1

-2

-3

f(x)=6x+x^2-x^3 -4
Relleno 1
x=-1; y<0 and y>-4 -5
25
 Calculo II 2012

y  f ( x)  6 x  x 2  x 3  Donde x   1,3
3  (1) 4
Primero calculamos: x  
n n
4
Ahora: ci  a  ix  ci  1  i
n
Y como: f ( x)  6 x  x 2  x 3
4 4 4 64 64 4
f (ci )  6( i  1)  ( i  1) 2  ( i  1)3   3 i 3  2 i 2  i  4
n n n n n n

Entonces recién ahora procedemos a hallar el área de la región acotada

n
Area  lim  f (ci )x
n 
i 1
n n
64 3 64 2 4 4 256 3 n 256 2 n 16 n
16
A  lim  ( 3
i  2
i  i  4)  lim[   4
i   3
i   2
i   ]
n  n n n n n i 1 n
i 1 i 1 n i 1 n i 1 n

256n 2 (n  1) 2 256n(n  1)(2n  1) 16n(n  1) 16n

A  lim[    ]
n  4n 4 6n 3 2n 2 n
1 128 1 1 1
A  lim[64(1  )  (1  )(2  )  8(1  )  16]
n  n 3 n n n
40 2
A u
3 8 y

4. Hallar el área de la región A acotada por la gráfica de y  3 x al eje  y a las rectas

7
x  0, x  4
y
6 6
X
5.5

5 5

4.5

4 4

3.5

A
3

a
2.5

2 2

1.5

1 1

0.5
x x
-1.5 -1 -0.5 0.5 1 1.5 2 2.5 3 3.5 4-1.5 4.5-1 5-0.5 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 6.5 7 7.5
-0.5

-1 -1

26
 Calculo II 2012

Para hacerlo más fácil elevamos al cuadrado a la función y tomamos como variable independiente
Y
y2
a y es decir: f ( y )  y como vemos en el grafico #2 la región está limitada entre las curvas
9
y2
f ( y)  , g ( y)  4 y las rectas y  0, y  6
9
Por ende mí f (ci )  [ g (ci )  f (ci )]
60 6
Primero calculamos: y  
n n
6
Ahora: ci  a  ix  0  i
n
y2 4
Y como: f ( y )  y g ( y)  4  g (ci )  4, f (ci )  2 i 2
9 n
Entonces recién ahora procedemos a hallar el área de la región acotada
n
Area  lim [ g (ci )  f (ci )]y
n 
i 1
n
4 2 6
A  lim  [4  i ]
n 
i 1 n2 n
n
24 n 24 2
A  lim[  3 i ]
n 
i 1 n i 1 n

24n 24n(n  1)(2n  1)

A  lim[  ]
n  n 6n 3
1 1
A  lim[24  4(1  )(2  )]
n  n n
A  16u 2

( x  4) 2 ( x  4)2 ( x  4) 2
5. Dada la región A acotada por la curvas y  ,y ,y ,
4 4 4
( x  4)2
y hallar su área.
4
y
6

a1 1
a2 x
-5 -4 -3 -2 -1 1 2 3 4 5 6

a4 a3 -1

-2
=((x-4)^2)/4
eno 1
=((x+4)^2)/4 -3
eno 2
=-((x-4)^2)/4
eno 3 -4
=-((x+4)^2)/4
eno 4
-5 a=a 1+a 2 +a 3+a 4

27
 Calculo II 2012

Si observamos la gráfica podemos percibir que existe simetría a los ejes y al origen de
coordenadas, entonces encontramos el área de una región y la multiplicamos por cuatro.

( x  4) 2
y  f ( x)   Donde x [0, 4]
4
40 4
Primero calculamos: x  
n n
4
Ahora: ci  a  ix  0  i
n
4
( i  4)2
( x  4) 2
4 8
Y como: f ( x)   f (ci )  n  2 i2  i  4
4 4 n n
Entonces recién ahora procedemos a hallar el área de la región acotada
n
Area  4 A2  4lim  f (ci )x
n 
i 1
n n
4 2 8 4 16 2 n 32 n
16
A  4 lim  ( 2
i  i  4)  4 lim[  3
i   2
i   ]
n  
i 1 n n n i 1 n i 1 n i 1 n
n

16n(n  1)(2n  1) 32n(n  1) 16n 8 1 1 1

A  4 lim[ 3
 2
 ]  4 lim[ (1  )(2  )  16(1  )  16]
n  6n 2n n n  3 n n n

64 2
A u
3

6. Hallar el área de la región A acotada por la gráfica de y  1  x 2  2 x 4 al eje  , al eje Υ, y

a la recta x  1
y

3.5

2.5

1.5

0.5
a
x
-3 -2.5 -2 -1.5 -1 -0.5 0.5 1 1.5
1+x^2+2x^4
no 1
y<4 and y>0 -0.5

28
 Calculo II 2012

y  f ( x)  1  x2  x4  Donde x [0,1]
1 0 1
Primero calculamos: x  
n n
1
Ahora: ci  a  ix  0  i
n
2 4
i i
Y como: f ( x)  1  x  x  f (ci )  1     2  
2 4

n n
Entonces recién ahora procedemos a hallar el área de la región acotada
n
Area  lim  f (ci )x
n 
i 1

  i 2
n
 i  1
4

A  lim  1     2   
 n  n   n
n 
i 1 

n
 1 i2 i4 
A  lim    2  2 4 
n 
i 1  n n n 
1n 1n(n  1)(2n  1) 2n(n  1)(6n3  9n 2  n  1) 
A  lim    
n  n
 6n 3 30n5 
 1  1  1  1  1  9 1 1 
A  lim 1  1   2    1   6   2  3  
n 
 6  n  n  15  n  n n n 
26
A  u2
15


x , x  2
2

7. Dada la región A acotada por las curva f  x    ,el eje  y las rectas
4 x  x , x  2

2

x  1, x  6 , calcular su área.
y
6

2
a2
a1 x
-11 -10 -9 -8 -7 -6 -5 -4 -3 -2 -1 1 2 3 4 5 6 7

-2

a3
-4

-6

-8

-10

-12

29
 Calculo II 2012

Según la gráfica: A  A1  A2  A3

a) Si x  2  y  x 2
2 1 1
Primero calculamos: x  
n n
1
Ahora: ci  a  ix  1 
n
2
 i
Y como: f  x   x  f  ci   1 
2

 n
n
A1  lim  f (ci )x
n 
i 1

 1 i 2 2i 
2
n
 i 1 n
A1  lim  1    lim    3  2 
i 1  n  n n i 1  n n n 
n 

 1n n  n  1 (2n  1) 2n  n  1 
A1  lim    
n 
n 6n 3 2n 2 
7
A1  u 2
3
b) Si x  2  y  4 x  x 2
42 2 2
Primero calculamos: x   Ahora: ci  a  ix  2  i
n n n
2
 2i   2i 
Y como: f  x   f  ci   4  ci    ci   4 2     2  
2

 n  n
n
A2  lim  f (ci )x
n 
i 1

n   2i   2i   2
2 n
 8 8i 2 
A2  lim   4  2     2     lim    3 
n 
 
i 1  n  n   n n i 1  n n 
 8n 8n  n  1 2n  1   4  1  1 
A2  lim     lim 8  1   2   
 3  n  n 
n  3 
n 6n  n

16
A2  u 2
3

c) Como A3 está debajo del eje  , se toma su valor absoluto.

64 2
Primero calculamos: x  
n n
2
Ahora: ci  a  ix  4  i
n

30
 Calculo II 2012

2
 2i   2i 
Y como: f  x   f  ci   4  ci    ci   4  4  4 
2

 n  n
n
A3  lim  f (ci )x
n 
i 1

n   2i   2i   2
2

A3  lim   4  4     4   
n 
 
i 1  n  n   n
n
 16i 8i 2 
A3  lim    2  3 
n 
i 1  n n 
 n 16i n 8i 2 
A3  lim   2   3 
n 
 i 1 n i 1 n 

16n  n  1 8n  n  1 2n  1 
A3  lim  2
 
n 
 2n 6n 3 
32 2
A3  u
3
Entonces el área pedida es:

7 16 32
A  A1  A2  A3   
3 3 3
55 2
A u
3
x
8. Calcular el área de la región A limitada por las gráficas : y  e 2
, al eje  y las rectas
x  0, x  1
y

1.5

0.5
A
x
-2.5 -2 -1.5 -1 -0.5 0.5 1 1.5

-0.5

-1

31
-1.5
 Calculo II 2012

x
y  f ( x)  e 2  Donde x   0,1
1 0 1
Primero calculamos: x  
n n
1 1
Ahora: ci  a  ix  ci  0  i  i
n n
x i
Y como: f ( x)  e 2 f (ci )  e 2 n

Entonces recién ahora procedemos a hallar el área de la región acotada

n
Area  lim  f (ci )x
n 
i 1

Por propiedad telescópica de la sumatoria:

n
 i2n1 2in 
  e  e   f  n   f  0
i 1  

 21n  n  2in  n2n1 21n

 e  1   e   e  e
  i 1  
1
 1 
e 2 n  e 2  1
1 n 2in  
e
n i 1

 21n 
 e  1 n
 
1
 1 
e  e 2  1
2n

i n 2in  
A  lim
n  n

i 1
e  lim
n   1 
 e  1 n
2n

 

1
Si z  ,n    z  0
n
z
 1 
e 2  e 2  1
A  2 lim  z z
z /20   2
 e  1
2

 

32
 Calculo II 2012

Aplicando la regla de L´Hospital:

ze z
lim 1
z 0 e z  1

Aplicamos esta igualdad a nuestra ecuación:

z
z
e2  
 1
  1 
A  2  e 2  1 lim z   2  e 2  1 u 2
2
  z /20  e 2  1  
 
 

  
9. Calcular el área de la región bajo la gráfica de y  cos  x  en   , 
 2 2

1.5

0.5

x
-4π/5 -7π/10 -3π/5 -π/2 -2π/5 -3π/10 -π/5 -π/10 π/10 π/5 3π/10 2π/5 π/2 3π/5

-0.5

-1

  
y  f ( x)  cos  x   Donde x    , 
 2 -1.52 
 
 ( )
cos(x) Primero calculamos: x-2 2 2 
no 1 n n
 
Ahora: ci  a  ix -2.5 ci    i
2 n
   
Y como: f ( x) -3cos  x  f (ci )  cos    i 
 2 n 

33
 Calculo II 2012

Entonces recién ahora procedemos a hallar el área de la región acotada

n
Area  lim  f (ci )x
n 
i 1
n
    
A  lim  cos    i 
n 
i 1 n  2 n 

     
Donde  cos    i   sen  i  y por telescópica obtenemos:
 2 n  n 

     n   

 1  cos  n   cos  n   cos  n  1 n  
     
A  lim 
n   n   
 2 sen   
 n 
 
 
  n   
 1  cos    cos    cos  n  1  
n  n 
A  lim  
n 
n    
 2 sen   
 n 
  
 n 
1111
A  2u 2
2

34
 Calculo II 2012

METODOS NUMERICOS PARA CALCULAR INTEGRALES DEFINIDAS

1.-Metodo Rectángulo
Si una función ƒ es continua en el intervalo *ɑ,Ƅ+ , entonces ƒ es integrable en [ɑ,Ƅ+ . Sin
embargo, de la proposición ya mencionada son sumamente generales, no se tiene garantía
de que, al aplicar los métodos usualmente conocidos para resolver integrales, se pueda
encontrar la antiderivada de una función ƒ(x) cualquiera, necesariamente para obtener la
integral definida. Para poder calcular estas dichas integrales se usa una de las técnicas
básicas que permiten resolver dicha situación, a través de la denominada “INTEGRACION
APROXIMADA, POR EL METODO DE LOS RECTANGULOS”.
Método Rectangular: consiste en dividir el área que se desea encontrar en n sub-áreas en
forma de rectángulos. Para el desarrollo del modelo se toman como referencia las
siguientes variables:
n : Numero de sub-áreas en las cuales se divide el área a calcular
 ó dx : Ancho o base de cada sub-área
li ó a : límite inferior definido para el cálculo del área
ls ó b : límite superior definido para el cálculo del área

Integración numérica de una función por el método de rectángulos:

La integral definida entre los puntos a y b de una función continua y acotada ƒ(x) representa el
área comprendida debajo de esa función. En ocasiones es necesario calcular integrales (áreas) de
modo numérico, es decir, sin conocer la integral explicita de la función ƒ(x). Existen varios posibles
métodos para calcular esta área. Quizás el más sencillo sea sustituir el área por un conjunto de n
sub-áreas donde cada sub-área semeja un pequeño rectángulo elemental de base dx = (b-a)/n y
altura h, el área seria:

Figura (1)

35
 Calculo II 2012

De la Figura (1)
Área=h*  donde h es el valor de la función calculada en el punto medio del área, ósea

f  a  dx / 2  y  ó dx es el ancho definido para dicha sub-área. Nótese que entre más grande

es dx entonces mayor será el área que se quita y pone al área real que se desea calcular (área
colocada entre la función y las línea azul y área quitada entre la función y la línea roja). Si se toma
li como límite inferior definido para el cálculo de la integral entonces el punto sobre el eje x para el
cálculo de h será:
li  x / 2 .
Teniendo en cuenta lo anterior el área será:
area  x * f (li  dx / 2) .

Si el área que se desea calcular se divide entre n sub-áreas, donde cada una de ellas representa un
pequeño rectángulo, entonces el área total será:
Área de primer rectángulo:
A1  h1 * dx  h1 * x

Donde h1 será la función evaluada en la mitad de la sección del primer rectángulo, se podría decir

entonces en términos generales que h1 es igual a la función evaluada en

x1 , h1  f ( x1 )

Teniendo en cuenta la anterior se deduce que x1  (li  dx / 2) y por lo tanto el área de ese

primer rectángulo será:

x
A1  x * f ( x1 )  x * f (li  )
2
Del mismo modo se puede Decir que el área del segundo rectángulo es:
x
A2  x / f ( x2 )  x * f (li  3* )
2
Área del tercer rectángulo es:
x
A3  x * f ( x3 )  x * f (li  5* )
2
Área del i-esimo rectángulo es:
x
Ai  x * f ( xi )  x * f (li  (2* i  1) )
2

36
 Calculo II 2012

Área total que será la sumatoria de todas las áreas parciales y quedara así:

x
ls n
At   f ( x)dx  x *  f (li  (2* i  1)* ) (1)
li i 1 2
La representación gráfica de esta forma de aproximar la integral se presenta en la Fig.2. Resulta
que si n se hace muy grande ( dx muy pequeño) el error Sera pequeño.

Figura 2

2.- Método del trapecio

A menudo es necesario evaluar la integral definida de una función que no tiene una antiderivada
explicita, o cuya antiderivada no es fácil de obtener, por eso el método para obtener una integral
definida simple recibe el nombre de cuadratura numérica y emplea una suma de tipo:
n

 f (c )  x
i 1
i i

La regla del trapecio nos da se aplica cuando el f(x) es una función continua en [a, b] la integral
definida está dada por:
n
f ( x)dx  lim  f (ci )i x
b
 a x 0
i 1

Para aproximar la medida del área de una región, usaremos trapecios en vez de rectángulos, para
este caso también usaremos particiones regulares y evaluaremos la funciones en los puntos cuyas
distancias sean las mismas

37
 Calculo II 2012

 f  x  dx , al intervalo [a, b] dividiremos en n sub-intervalos es entonces

b
En la integral definida
a

que; X1 = a , X2 =a + ∆x , ………. , Xn = b :

 f  x  dx expresaremos como la suma de n integrales definidas:

b
Luego la integral definida
a

Xn 1
 f  x  dx   f  x  dx   f  x  dx  .....   f  x  dx   f  x  dx
b X1 X2 X

a a X1 Xn  2 Xn 1

f  x  dx , da la medida del área de la región acotada por las rectas X=a, X=X1 y la
X1
La integral  a

f  x  dx , se la puede aproximar por la

X1
porción de la curva PoP1 . Luego a la integral definida a

medida del trapecio formado por las rectas X=a, X=X1, PoP1 y el eje X donde la medida de este
trapecio es:
1
 f  Xo   f  X 1  x
2
Teniendo en cuenta que las demás integrales se pueden aproximar la medida del área de un
trapecio; entonces para la n-esima integral se tiene:
1
f  X  dx   f  Xi  1  f  Xi   X
Xi
Xi 1 2

38
 Calculo II 2012

 f  x  dx , se tiene:
b
Por lo tanto para la integral definida
a

X
 f  x  dx  2  f  Xo   2 f  X 1  2 f  X 2  ....  2 f  Xn 1  f  Xn 
b

Esta sería la regla del trapecio:

Una integral definida por la regla del trapecio tiene las siguientes características:
 Exactitud
 Error

2.1.- Exactitud.- la exactitud de una integral definida por la regla del trapecio se obtiene cuando ∆X
se aproxima a cero (∆X  0) y n crece sin límite, el límite de la aproximación por la regla del
trapecio es el valor exacto de la integral definida; es decir:

1
T   f  X 1  f  X 2   ....  f  Xn   X   f  Xo   f  Xn  X
2
n
1
T   f  Xi X   f  a   f  b   X
i 1 2
n
1
lim T  lim  f  Xi X  lim  f  a   f  b   X
x 0 x  0 x  0 2
i 1

Por lo tanto la exactitud queda definido como:

lim   f  x  dx
b

x 0 a

2.2.- Error.- Al aplicar la ley de los trapecios es posible que se cometan errores, que denotaremos
por “ε” y que pueden ser hallados mediante la siguiente condición:

Sea f una función continua en el intervalo [a,b+ y que, f‘ , f‘‘ existen en *a,b+:
Si:

 T   f  x  dx  T
b

 f  x  dx , que se encontró mediante la Regla del Trapecio,

b
Donde T es el valor aproximado de
a

entonces existe un número “ŋ” en *a, b+ tal que:

 b  a  f ''    X 
1
T  
2

12

39
 Calculo II 2012

3.- Método de Simpson

Otro método para aproximar el valor de una integral definida es proporcionado por la regla de
Simpson también conocida con el nombre de regla parabólica. Para una partición dada del
intervalo cerrado [a,b], la regla de Simpson usualmente da una aproximación mejor que la regla
trapecial. Sin embargo, la regla de Simpson requiere más esfuerzo para ser aplicada. En la regla
trapecial los puntos sucesivos en la gráfica de y  f (x) son conectados por segmentos de líneas
rectas, mientras que en la regla de Simpson los puntos se conectan por segmentos de parábolas.

La regla de Simpson se divide en dos:

3.1.-Regla de Simpson 1/3

Como ya mencionamos anteriormente la regla de Simpson nos da una aproximación mejor que
la regla del trapecio ya que al utilizarla nos da un resultado con menor error, pero antes de
desarrollar la regla de Simpson es necesario conocer y demostrar el siguiente teorema.

3.1.1.-Teorema 1: Si P0 ( x0 , y0 ), P1 ( x1 , y1 ) y P2 ( x2 , y2 ) son tres puntos no colineales en la

parábola que tiene la ecuación y  Ax 2  Bx  C , donde y0  0, y1  0, y2  0, x1  x0  h
y x2  x0  2h , entonces la medida del área de la región acotada por la parábola , el eje x y
las rectas x1  x0 y x2  x1 está dada por:

1
h( y 0  4 y1  y 2 ) ………….…(1)
3

DEMOSTRACIÓN: La parábola cuya ecuación es y  Ax 2  Bx  C tiene un eje vertical (ver

Fig. 1)

( x2 , y2 )
y
( x1 , y1 ) P2
P1
( x0 , y 0 )
P0

Fig. 1: parábola acotada

x
h h
40
 Calculo II 2012

La figura muestra la región acotada por la parábola, el eje x y las rectas x1  x0 y x2  x1 .

Ya que P0 , P1 y P2 son puntos en la parábola, sus coordenadas satisfacen la ecuación de la

parábola. Así cuando reemplazamos x1 por x0  h y x 2 por x0  2h tenemos:

y0  Ax 0  Bx 0  C
2

y1  A( x0  h) 2  B( x0  h)  C  A( x02  2hx0  h 2 )  B( x0  h)  C

y 2  A( x0  2h) 2  B( x0  2h)  C  A( x02  4hx0  4h 2 )  B( x0  2h)  C

Por lo tanto:

y0  4 y1  y2  A(6 x02  12hx0  8h 2 )  B(6 x0  6h)  6C ……....(2)

Si S (en unidades cuadradas) es el área de la región acotada por la parábola, entonces S puede
ser calculada por el límite de una suma de Riemann, de la siguiente manera:

n
S  lim  ( A i2  B i  C )x
x 0
i 1

Aplicamos el límite:

x0  2 h
S ( Ax 2  Bx  C )dx
x0

x 2h

0
1 3 1 2
S  Ax  Bx  Cx 
3 2  x0

1 1 1 1
S A( x0  2h) 3  ( x0  2h) 2  C ( x0  2h)  ( Ax 03  Bx 02  Cx0 )
3 2 3 2

1
 
S  h A(6 x02  12hx0  8h 2 )  B(6 x0  6h)  6C …….…(3)
3

Reemplazando (2) en (3) se tiene:

S  h y 0  4 y1  y 2 
1
3

Con esta última ecuación queda demostrado el Teorema 1.

41
 Calculo II 2012

Ahora desarrollaremos la regla de Simpson.

Sea una función f continua en el intervalo cerrado [a,b]. Consideramos una partición regular del
intervalo [a,b] de 2n subintervalos (2n se usa en lugar de n porque queremos que sea un número
(b  a)
par de subintervalos ). La longitud de cada subintervalo está dada por x  . Denotemos a
2n
los puntos en la curva que tienen estos puntos de la partición como abscisas por P0 ( x 0 , y 0 ),
P1 ( x1 , y1 ) ,..., P2n ( x2n , y 2n ) ; ver Fig. 2 donde f ( x)  0 para toda x en [a,b].

y  f (x)

P1 P2 n1 P2 n
P0
P2 n2
P2 P5
P3 P4

............................
x
a  x0 x1 x2 x3 x4 x5 ............................ x2 n2 x2 n1 x2 n  b
Fig. 2

Aproximamos el segmento de la curva y  f (x) de P0 a P2 por el segmento de la parábola con

su eje vertical a través de P0 , P1 y P2 . Entonces por el Teorema 1 la medida del área por la región
acotada por esta parábola, el eje x y las rectas x  x0 y x  x2 con h  x , esta dad por:

x y0  4 y1  y2  x f ( x0 )  4 f ( x1 )  f ( x2 )
1 1
o
3 3

En forma análoga aproximamos el segmento de la curva y  f (x) de P2 a P4 por el segmento de

la parábola con su eje vertical a través de P2 , P3 y P4 . La medida del área de la región acotada
por esta parábola, en el eje x y las rectas x  x2 y x  x4 esta dad por:

x y2  4 y3  y4  x f ( x2 )  4 f ( x3 )  f ( x4 )
1 1
o
3 3

42
 Calculo II 2012

Este procedimiento se continúa hasta que tengamos n regiones y la medida del área de la última
región está dada por:

x y2 n2  4 y2 n1  y2 n  x f ( x 2 n  2 )  4 f ( x 2 n 1 )  f ( x 2 n )

1 1
o
3 3

La suma de las medidas de las áreas de estas regiones aproxima la medida del área de la región
acotada por la curva cuya ecuación es y  f (x) el eje x y las rectas x  a y x  b . La medida
b
del área de la región está dada por la integral definida a
f ( x) dx . Así tenemos como una
aproximación de la integral definida.

x f ( x0 )  4 f ( x1 )  f ( x 2 )  x f ( x 2 )  4 f ( x3 )  f ( x 4 )  ...  x f ( x 2 n  2 )  4 f ( x 2 n 1 )  f ( x 2 n )

1 1 1
3 3 3

Agrupando todas las funciones se tiene:

f ( x)dx  x f ( x 0 )  4 f ( x1 )  2 f ( x 2 )  4 f ( x3 )  2 f ( x 4 )  .....  2 f ( x 2 n  2 )  4 f ( x 2 n 1 )  f ( x 2 n )

b 1
a 3

1
La cual se le denomina regla de Simpson 1/3 por el factor x
3

ba
Dónde: x 
2n

Anteriormente quedo demostrado que la regla de Simpson es un método más preciso que otros,
pero eso no quiere decir que no tenga un error de aproximación.

El siguiente teorema que se demuestra en el cálculo avanzado, da un método para determinar el

error al aplicar la regla de Simpson denotado por:  s

3.1.2.-Teorema 2: Sea la función f continua en el intervalo cerrado a, b , y f ' , f ' ' , f ''' y f iv
existe en a, b . Si

b
 s   f ( x)dx  Rs
a

b
Donde Rs es el valor aproximado de a
f ( x)dx que se encontró por la regla de Simpson,
entonces existe algún número n en a, b tal que

s  
1
b  a  f iv (n)x 4
180

43
 Calculo II 2012

OBSERVACIÓN: Cuando f (x) sea un polinomio de tercer grado o menos, entonces f iv ( x)  0

por lo tanto  s  0 , lo que ocasionará utilizar la siguiente fórmula práctica:

b ba  ab 
 a
f ( x)dx 
6  f ( a )  4 f (
2
)  f (b)


A continuación veremos algunos ejemplos de integración numérica utilizando la regla de Simpson.

1

2
Ejemplo 1: Utilizar la regla de Simpson para aproximar la siguiente integral e x dx con 2n  4
0
y dar el resultado con cuatro cifras decimales.

Tengamos en cuenta que el valor real de la integral con cuatro decimales es:

1
 e dx  1.4626
2
x
0

SOLUCIÓN: Aplicando la regla de Simpson con 2n  4 , tenemos

x 
1
1  0  1
4 4

i
f ( x)  e x , x0  0, xi  x0  ix 
2

1
 f ( x0 )  4 f ( x1 )  2 f ( x2 )  4 f ( x3 )  f ( x4 ).....(1)
1
 e x dx 
2

0 12
i) Realizamos nuestra tabla de valores

Tabla 1

i xi f ( xi ) k k  f ( xi )

0 0 1 1 1

1 0.25 1.0644 4 4.2576

2 0.5 1.2840 2 2.568

3 0.75 1.7550 4 7.02

4 1 2.7182 1 2.7182

suma 17.5638

44
 Calculo II 2012

ii) Reemplazamos en (1)

1
17.5638  1.4636
1
e dx 
x2
0 12

2
 xe x dx
3
Ejemplo 2: Utilizar la regla de Simpson para aproximar la siguiente integral con
1
2n  4 y dar el resultado con cuatro cifras decimales.

Tengamos en cuenta que el valor real de la integral con cuatro decimales es:

1
 xe x dx  1.3763
3
0

SOLUCIÓN: Aplicando la regla de Simpson con 2n  4 , tenemos

x 
1
1  0  1
4 4

f ( x)  3 x e x

i
x0  0, xi  x0  ix 
4

1
 f ( x0 )  4 f ( x1 )  2 f ( x2 )  4 f ( x3 )  f ( x4 ).....(2)
1
 xe x dx 
3
0 12

i) Realizamos nuestra tabla de valores

Tabla 2

i xi f ( xi ) k k  f ( xi )

0 0 0 1 0

1 0.25 0.8088 4 3.2352

2 0.5 1.3085 2 2.617

3 0.75 1.9234 4 7.6936

4 1 2.7182 1 2.7182

suma 16.264

45
 Calculo II 2012

ii) Reemplazamos en (2)

1
16.264  1.3553
1
e dx 
x2
0 12

2 dx
Ejemplo 3: Utilizar la regla de Simpson para aproximar la siguiente integral 1
1  x3
con

2n  4 y dar el resultado con cuatro cifras decimales.

Tengamos en cuenta que el valor real de la integral con cuatro decimales es:

2 dx

1
1  x3
 0.4925

SOLUCIÓN: Aplicando la regla de Simpson con 2n  4 , tenemos

x 
1
2  1  1
4 4

1
f ( x) 
1  x3

i
x0  1, xi  x0  ix  1 
4

dx 1
 f ( x0 )  4 f ( x1 )  2 f ( x2 )  4 f ( x3 )  f ( x4 ).........(3)
2

1
1 x 3

12

i) Realizamos nuestra tabla de valores

Tabla 3

i xi f ( xi ) k k  f ( xi )

0 1 0.7071 1 0.7071

1 1.25 0.5819 4 2.3277

2 1.5 0.4781 2 0.9562

3 1.75 0.3965 4 1.5862

4 2 0.3333 1 0.3333

suma 5.9105

46
 Calculo II 2012

ii) Reemplazamos en (3)

dx 1
5.9105  0.4925
2

1
1 x3

12
Lo que nos da un valor exacto con cuatro cifras decimales de la integral. Esto sucede porque
f (x) es un polinomio de tercer grado o menor.
3 ln( x)
Ejemplo 4: Utilizar la regla de Simpson para aproximar la siguiente integral 1
x3
.dx con

2n  4 y dar el resultado con cuatro cifras decimales.

Tengamos en cuenta que el valor real de la integral con cuatro decimales es:
3 ln( x)

1
x3
.dx  0.4220

SOLUCIÓN: Aplicando la regla de Simpson con 2n  4 , tenemos

x 
1
3  1  1
4 2

ln( x) i
f ( x)  x0  1, xi  x0  ix  1 
x3 , 2

ln( x) 1
 f ( x0 )  4 f ( x1 )  2 f ( x2 )  4 f ( x3 )  f ( x4 ).........(4)
3

1
x 3
.dx 
6

i) Realizamos nuestra tabla de valores

Tabla 4

i xi f ( xi ) k k  f ( xi )

0 1 0.0000 1 0.0000

1 1.5 0.2207 4 0.8828

2 2 0.2451 2 0.4901

3 2.5 0.2318 4 0.9272

4 3 0.2114 1 0.2114

suma 2.5116

47
 Calculo II 2012

ii) Reemplazamos en (4)

ln( x) 1
2.5116  0.4186
3

1
x3
.dx 
6
1 dx
Ejemplo 5: Utilizar la regla de Simpson para aproximar la siguiente integral 0
1 x2
con 2n  4

y dar el resultado con seis cifras decimales. Tengamos en cuenta que el valor real de la integral con
seis decimales es:
1 dx

0
1 x2
 0.881373

SOLUCIÓN: Aplicando la regla de Simpson con 2n  4 , tenemos

x 
1
1  0  1
4 4

1
f ( x) 
1 x2

i
x0  0, xi  x0  ix 
4

dx 1
 f ( x0 )  4 f ( x1 )  2 f ( x2 )  4 f ( x3 )  f ( x4 ) .........(5)
1

0
1 x2

12

i) Realizamos nuestra tabla de valores

Tabla 5

i xi f ( xi ) k k  f ( xi )

0 0 1 1 1

1 0.25 0.970142 4 3.880570

2 0.5 0.894427 2 1.788854

3 0.75 0.800000 4 3.200000

4 1 0.707106 1 0.707106

suma 10.576531

48
 Calculo II 2012

ii) Reemplazamos en (5)

dx 1
10.576531  0.881377
1

0
1 x 2

12

iii) Calculamos el error de la regla de Simpson  s

s  
1
b  a  f iv (n)x 4
180
9
 f ( x)  105 x (1  x ) , Como x  0;1
1 
f ( x)  iv 4 2 2

1  x2

f iv (0)  3 , f iv (1)  1.01646

4
1 1
Para: n  0, Es   (1)(3)   6.5104110 5
180 4
4
1
(1) 1.01646   2.20587  10 5
1
n  1, Es  
180 4

0.00006  Es  0.00002

1 dx
Ejemplo 6: Utilizar la regla de Simpson para aproximar la siguiente integral 1 x
0
con 2n  4 y

dar el resultado con cinco cifras decimales. Tengamos en cuenta que el valor real de la integral con
cinco decimales es:

1 dx
 1  x  0.69314
0

SOLUCIÓN: Aplicando la regla de Simpson con 2n  4 , tenemos

x 
1
1  0  1
4 4

1
f ( x) 
1 x

i
x0  0, xi  x0  ix 
4

 1  x  12  f ( x )  4 f ( x )  2 f ( x )  4 f ( x3 )  f ( x4 ).........(6)
1 dx 1
0 1 2
0

49
 Calculo II 2012

i) Realizamos nuestra tabla de valores

Tabla 6

i xi f ( xi ) k k  f ( xi )

0 0 1.00000 1 1.00000

1 0.25 0.80000 4 3.20000

2 0.5 0.66667 2 1.33333

3 0.75 0.57143 4 2.28571

4 1 0.50000 1 0.50000

suma 8.31905

ii) Reemplazamos en (6)

 1  x  12 8.31905  0.69325
1 dx 1
0

iii) Calculamos el error de la regla de Simpson  s

s  
1
b  a  f iv (n)x 4
180

 f iv(x)  1201  x  , Como x  0;1

1 6
f(x) 
1 x

3
f iv (0)  24 , f iv (1) 
4
4
1 1
Para: n  0, E s   (1)(24)   0.00052
180 4
4
1  3  1 
n  1, E s   (1)    0.00002
180  4  4 

 0.00052  Es  0.00002

50
 Calculo II 2012

4.- Integración por series de Taylor

Esta serie permite obtener una muy buena aproximación polifónica de una función en un entorno
de cierto punto en que la función sea diferenciable. Además la serie permite acotar el error
obtenido mediante dicha estimación. Por supuesto, para hacer esta aproximación sólo se pueden
tomar unas cuantas expresiones de esta serie, por lo que el resto resulta en un error conocido
como el término residual.

Dicha serie es la siguiente:

f '(a) f ''(a) f ( n ) (a)

f ( x)  f (a)  ( x  a)  ( x  a)  ... 
2
( x  a) n  Rn ( f )
1! 2! n!

O en forma compacta:


f ( n ) (a)
f ( x)   ( x  a)n  Rn ( f )
n o n!

Dónde:

n! : factorial de n.

f (n)(a) : La n-ésima derivada de f en el punto a.

Rn(f): Residuo.

f ( n1) ( )
Rn ( f )  ( x  c) n1  [c, x]
(n  1)!

Demostración de la serie:

La representación en serie de potencias de f(x).

f ( x)  a0  a1 ( x  c)  a2 ( x  c)2  ...  an ( x  c) n

La primera derivada:

f '( x)  a1  2a2 ( x  c)1  3a3 ( x  c)2  ...  nan ( x  c)n1

Evaluamos en x=c:

f '(c)  a1

La segunda derivada:

f ''( x)  2a2  3.2a3 ( x  c)1  4.3a4 ( x  c)2 ...  n(n  1)an ( x  c)n2

51
 Calculo II 2012

Evaluamos en x=c:

f '(c)  2a2

La tercera derivada:

f '''( x)  3.2a3  4.3.2a4 ( x  c)1  5.4.3a4 ( x  c) 2 ...  n(n 1)(n  2)an ( x  c) n3

Evaluamos en x=c:

f '''(c)  3.2a3

La cuarta derivada:

f 'v ( x)  4.3.2a4  5.4.3.2a4 ( x  c)1...  n(n 1)(n  2)(n  3)an ( x  c) n4

Evaluamos en x=c:

f 'v (c)  4.3.2a4

En general:

f ( n ) (c)  n!an

Despejamos an:

f ( n ) (c )
an  …(*)
n!

f ( x)  a0  a1 ( x  c)  a2 ( x  c) 2  ...  an ( x  c) n   an ( x  c) n …(**)
n 0

Reemplazamos (*) en (**) y nos queda la serie de Taylor centrada en c:


f ( n ) (c )
f ( x)   ( x  c) n
n 0 n!

Veamos un ejemplo donde nos muestre como la gráfica de la serie se acerca a la función:

52
 Calculo II 2012

Ejemplo 01:

f ( x)  cos( x) , f ´( x)   s en( x) , f ´´( x)  co s( x) , f ´´´( x)  s en( x) , f ´v ( x)  co s( x) ,

f v ( x)   s en( x)

Si evaluamos en c=0 para centrar la función en 0:

f (0)  1 , f ´(0)  0 , f ´´(0)  1 , f ´´´(0)  0 , f ´v (0)  1 , f v (0)  0

Nos quedaría:

1 0 0 1 1 2 0 3 1 4
cos( x)  * x  * x  * x  * x  * x  ....
0! 1! 2! 3! 4!

f ( n ) (o ) n
cos( x)   *x
n 0 n!

Ahora veamos cómo funcionan los polinomios de Taylor:

0
f ( n ) (0) n
cos( x)   *x
n 0 n!

2
f ( n ) (0) n
cos( x)   *x
n 0 n!

53
 Calculo II 2012

4
f ( n ) (0) n
cos( x)   *x
n 0 n!

6
f ( n ) (0) n
cos( x)   *x
n 0 n!

8
f ( n ) (0) n
cos( x)   *x
n 0 n!

54
 Calculo II 2012

10
f ( n ) (0) n
cos( x)   *x
n 0 n!

Obs:

 Bueno al hacer la serie de Taylor se nos va a hacer mucha más fácil integrar algunas
funciones que no se puedan integrar con los métodos básicos.
 Como pudimos observar en el primer ejemplo mientras más grande sea la sumatoria más
cerca estaremos de la función.
 Al centrar la función en 0, estamos aplicando la serie de Maclaurin, la cual nos puede ser
útil para encontrar integrales las cuales no podamos sacar primitivamente.

Ahora veamos algunos ejemplos más:

Ejemplo 02:

Evalué la integral utilizando series.

e
x3
dx

Solución:

Primero hallemos la serie para ex y luego reemplazamos x por x3.

f ( x)  e x  f (0)  e0  1
f ´( x)  e x  f ´(0)  e0  1
f ´´( x)  e x  f ´´(0)  e0  1
f ´´´( x)  e x  f ´´´(0)  e0  1

55
 Calculo II 2012

Quedaría:

1 0 1 1 1 2 1 3
f ( x)  x  x  x  x  ....
0! 1! 2! 3!
 n
x
f ( x)  
n 0 n !

Ahora reemplacemos x por x3:


( x3 ) n  x3n
f ( x)   
n 0 n ! n 0 n !

Integrando:


x3n
 f ( x)dx    dx
n 0 n!


x3n 1
 f ( x)dx   e x dx  
3

n  0 (3n  1) n !

Ejemplo 03:

Integrar f(x) = x2e-x

Solución:

Del ejemplo 1 se sabe que:


xn
e x

n 0 n !

Ahora le damos forma:


(  x) n 
x 2 ( x)n  x 2 n (1) n
e x   x 2e x   
n 0 n! n 0 n! n 0 n!

Integrando:


x 2 n (1)n 
x3n (1)n
 x e dx    dx  seria : 
2 x

n 0 n! n  0 ( n  3) n !

56
 Calculo II 2012

Ejemplo 04:

Utilice la serie de Taylor para senx en pi/4 para calcular el sen(47°) con una exactitud de 4 cifras
decimales.

Solución:

   
f( ) f ´( ) f ´´( ) f ´´´( )
s en( x)  4 * ( x   )0  4 *( x   )1  4 *( x   ) 2  4 *( x   )3  ...
0! 4 1! 4 2! 4 3! 4

47 1 1
x   x   
180 4 90

47 1 1  1 1  1 1 
s en( )  2 2 *( )  2 * ( )2  2 * ( )3  ...
180 2 2 90 2 2 90 2 6 90
47 1
s en( )  2(1  0.03490  0.00061  0.000002  ...)
180 2

sen(47)  0.73136

Ejemplo 05:

1 senx
Integrar: 
0.5 x
dx

Solución:

No puede obtenerse una antiderivada de la integral en términos de funciones elementales. Sin

embargo de la serie Taylor con centro en c=0 para sen(x) se tiene:

senx 1
 senx
x x
senx 1 x3 x5 x 7 x9
 ( x      ...)
x x 3! 5! 7! 9!
2
senx x x 4 x 6 x8
 1      ...
x 3! 5! 7! 9!
Lo cual es cierto para todo x diferente de 0. Al emplear la integración término a término resulta:

senx
1 x3 x5 x7 x9
0.5 x dx  x  3.3!  5.5!  7.7!  9.9!  ...|0.5
1

=(1-0.0555555+0.0016667-0.0000283+0.0000003)-(0.5-0.0069444+0.0000521-0.0000002)
1 senx
 0.5 x
dx  0.45298

57
 Calculo II 2012

SUMATORIAS DOBLES

1.-DEFINICION

Una sumatoria doble es un caso particular de la suma, en el cual el termino general ak es a su vez
una sumatoria, para cada k. veremos entonces como sumar más de un índice.

Si tenemos la siguiente sumatoria:

b
i 0
i …(1)

bi A su vez es otra sumatoria, donde:

m
bi   aij … (2)
j 0

De (1) y (2):

n m

 a
i 0 j 0
ij

Dónde: aij podrá depender de ambos índices y los limites inferior y superior de a
j 0
ij podrán

depender del índice i.

Para un mejor entendimiento analizaremos los términos usados pero acomodados de la siguiente
manera:

a00 a01 a02 ... a0 m

a10 a11 a12 ... a1m
.
. . . .
.
an 0 an1 an 2 ... anm

Si sumamos los miembros de la primera fila obtendremos: a

j 0
0j

58
 Calculo II 2012

Si sumamos los miembros de la segunda fila obtendremos: a

j 0
1j

m
 n  m 
a
n m
Si sumamos los miembros de la n-esima fila obtendremos:  ai b j    ai    b j  nj
i 0 j 0  i 1   j 1  j 0

n  m 
Y si sumamos todas estas sumas obtendremos:    aij 
i 0  j 0 
m
 n 
Análogamente con las columnas obtendremos el siguiente resultado:    aij 
j 0  i 0 
Observamos que los resultados anteriores son iguales:
n  m  m  n 
   aij  =    aij 
i 0  j 0  j 0  i 0 
Notamos entonces que la suma doble es el resultado de sumar cada fila a la vez, obviamente que
también es igual a sumar cada columna a la vez.

Con lo mencionado anteriormente tenemos ya un concepto claro de la definición de sumatoria

doble de las cuales desprenderemos las siguientes propiedades.

2.-PROPIEDADES DE LA SUMATORIA DOBLE

n  m  m  n 
   aij  =    aij  ; demostrada anteriormente.
i 0  j 0  j 0  i 0 
n m

a)  k  nmk , k  R
i 0 j 0
n m n m

b)  ka
i 0 j 0
ij  k  aij , k  R
i 0 j 0

  a  b   n b  m a
n m m n

c) i j j i
i 0 j 0 j 0 i 0

n m
 n  m 
d) 
i 0 j 0
ai b j    ai    b j 
 i 1   j 1 

59
 Calculo II 2012

2.1.-DEMOSTRACION DE LAS PROPIEDADES DE LA SUMATORIA DOBLE

n m
 n
 n
Demostrar: 
i 0 j 0
k    k  k  ...  k     mk 
i  0  `m `TERMINOS  i 0

 k  k  k  ...  k  mk
j 0 `m `TERMINOS
n mn
  n

 k   
i 0 
k  k  ...  k    mk 

i 0 j 0 " m "ter min os  i 0

  mk   mk  mk  ...  mk  nmk
i 0 " m"ter min os

n m n m
Demostrar:  ka
i 0 j 0
ij  k  aij , k  R
i 0 j 0

m m

 kaij  k  aij
j 0 j 0

n m n  m 
 kaij    k  aij 
i 0  j 0
i 0 j 0 
n

k a
i 0
i0  ai1  ...  aim    k  a00  a01  ...  a10  a11  ...  a20  a21  ...  anm 
n m
k  a00  a01  ...  a10  a11  ...  a20  a21  ...  anm   k  aij
i 0 j 0

  a  b   n b  m a
n m m n
Demostrar: i j j i
i 0 j 0 j 0 i 0

n  m 
 i j      ai  b j  
n m

 a  b 
i 0  j 0
i 0 j 0 
 m  n  
   ai  b j      mai   b j 
n m


i 0  j 0  i 0  j 0 
n  m  n m

  i  j
i 0 
ma  b  m  ai  n  bj
j 0  i  0 j  0

60
 Calculo II 2012

n m
 n  m 
Demostrar: 
i 0 j 0
a b
i j    ai    b j 
 i 1   j 1 

n m n  m 
 a b
i j   ai   b j  ai es una constante para la segunda sumatoria
i 0 j 0 k 0  j 0 
n m
 n  m 

i 0 j 0
a b
i j    ai    b j 
 k 0   j 0 

EJEMPLOS

1)
4  2 
Hallar     i  2 j  
i 1  j 1 
primero la sumatoria sobre j
4 4

  i  2 1    i  2  2      i  2    i  4 

i 1 i 1
4 4

  i  2    i  4     2i  6 

i 1 i 1
4

  2i  6    2 1  6    2  2   6    2 3  6    2  4   6 

i 1
4  2 
    i  2 j    44
i 1  j 1 
2)

Hallar   i.3 j 
5 4

i 1 j 1

operamos sobre j primero

5 5

 i.31  i.32  i.33  i.34    3i  9i  27i  81i 

i 1 i 1
5

 120i   120 1  120  2   120  3  120  4   120  5

i 1

  i.3   1800
5 4
j

i 1 j 1

61
 Calculo II 2012

3)
2
 ji 
3 4
Hallar   
j 1 i 1  j  i 

operando sobre i primero

3  j  2  2 j  2  3 j  2  4 j  2  3  j 2 4 j2 9 j2 16 j 2  
            
j 1   j  1
   

j 1  j  1             
2 2 2 2
  j 2   j 3   j 4  
  j 2 j 3 j 4 
operando sobre j
 12 4 1
2
9 1
2
16 1   1049 
2

          7.2    13.762
 1  1 1  2  1  3 1  4    225 
2 2 2 2


4)
4 5

  1 . j  j 
i 2
Hallar
i 1 j  2

operando sobre j

       
4 5 4

  1 . j  j
i 2      1i .  2    2 2   1i .  3   32   1i .  4    4 2   1i .  5    5 2 
i 1 j  2
 i 1  

   1 .  2    2     1 .3  3    1 . 4    4     1 . 5  5     2  .5  35
4 4
i 2 i 2 i 2 i 2 i

i 1 i 1

ahora sobre i

  2  .5  35    2  .5  35    2      

4
.5  35   2  .5  35   2  .5  35 
i 1 2 3 4

i 1

4 5

  1 . j  j   190

i 2

i 1 j  2


5)

3 9
Hallar   i  2  j  1
i 1 j  7

operando sobre j
3 9 3

  i  2  j  1     i  2  7  1     i  2 8  1     i  2  9  1 

i 1 j  7 i 1
3 3

   i  2  7  1     i  2 8  1     i  2  9  1     21i  42

i 1 i 1

operando sobre i
3

  21i  42   211  42    21 2   42    213  42 

i 1
3 9

  i  2  j  1  252
i 1 j  7

62
 Calculo II 2012

CALCULO DE LA INTEGRAL DOBLE CON SUMATORIAS DOBLES

1.-DEFINICIÓN

Ahora vamos a estudiar las integrales dobles donde nos ubicaremos en un plano xy
(bidimensional). Lo que sigue ahora es dar una notación nueva .En una sucesión ordinaria [a1],
cada término ha indizado por un entero i. La suma de todos los a desde i=1 hasta i=m se denota
por:

Pero vemos que aparecen 2 índices como por ejemplo:

aij= aij= o aij=

Para eso se nota la notación de doble sumatoria para denotar la suma de todos los términos
indizados. Mediante:

∑∑

Representaremos la suma de todos los términos donde i varía desde 1 hasta m y j varias desde
1 hasta n .Por ejemplo:

∑ ∑ =2.5+2. + .5+ =420

Dado que la adición es asociativa y conmutativa, podemos sumar los términos en cualquier orden.
Usualmente se hace:

∑ ∑ =∑ ∑

Podemos evaluar la expresión de la derecha desarrollando primero respecto de i y luego respecto

de j:

∑ ∑ =∑ ∑ +…+∑

= +

O podemos desarrollar primero respecto de j y luego respecto de i:

∑ ∑ =∑

=( +…+(

63
 Calculo II 2012

El resultado es el mismo

Por ejemplo, podemos escribir:

∑ ∑ =∑ +∑ ∑ =

Dados que las constantes se pueden sacar como factores de las sumas simples, también se puede
sacar como factores de las sumas dobles, concretamente,

∑ ∑ =α∑ ∑

Además:

∑ ∑ =∑ ∑ +∑ ∑

Las sumas dobles mas fáciles de manejar son aquellas en las cuales cada termino aparece
como un producto de , en la cual cada factor solo lleva un índice. En este caso podemos
expresar la suma doble como un producto de dos sumas simples:

∑ ∑ ∑ )(∑

2.-INTEGRAL DOBLE SOBRE UN RECTÁNGULO

Empezamos con una función f continua en un rectángulo

R: a c

Nuestro propósito es definir la integral doble de f sobre R:

∫∫ ,

64
 Calculo II 2012

Luego

(P)= (área de
)+ +

= ( ) ( ) ( ) ( ) ( ) =

En cada rectángulo ,f alcanza su valor mínimo en el punto ( que es el vértice

mas próximo al origen: , =

Luego

(P)= (area de

Volvamos ahora a la situación que es general. Como en el caso de una variable, se puede
demostrar que si f es continua, existe un número I y solo uno que satisface la desigualdad.

Para todas las particiones P de R

INTEGRAL DOBLE SOBRE UN RECTANGULO

Sea f una situación continua en un rectángulo cerrado R. El único número I que satisface la
desigualdad.

Para todas las particiones P de R

Se llama integral doble de f sobre R y se designa por

∫∫ ,

LA INTEGRAL DOBLE COMO VOLUMEN

Si f es continua y no negativa en un rectángulo R, la ecuación

Z=f(x,y)

Representa una superficie que está por encima de R. En este caso, la integral doble

∫∫ ,

65
 Calculo II 2012

Da el volumen del solido delimitado inferiormente por R y delimitado superiormente por la

superficie z=f(x, y)

Para verlo hay que considerar una partición P de R. Entonces P divide R en subrectángulos y
por tanto divide el solido Ten partes . Dado que contiene un paralelepípedo rectangular de
base y altura debe cumplirse

3.-PROPIEDADES

3.1.- Si f(x,y) es continua en la región R={a x b; g(x) y h(x)}, en la integral iterada doble de
f(x,y) se cumple:
b d n m

 f ( x, y )dxdy   dx  f ( x, y )dy  lim  f ( x , y j )xi y j

R n 
a c i 1 j 1
m 

Donde (xi , y j ) ϵ a la sub-región Rij.

3.2.-Para cada sub-región Rij = {(x,y) / x ϵ [xi-1 , xi], y ϵ [yj-1 , yj]} existe un supremo Mij y un ínfimo mij,
tal que:

Mij = sup {f (x, y) / (x,y) ϵ R ij} ; mij = inf {f(x,y) / (x, y) ϵ R ij}

Además de las suma de todos los productos Mij.∆xi ∆yj Y de todos los productos mij∆xi∆yj, se
obtiene respectivamente:

66
 Calculo II 2012

n m
U f ( P)   M ij xi y j
i 1 j 1
n m
L f ( P)   mij xi y j
i 1 j 1

Siendo Uf (P) la suma superior y Lf (p) la suma inferior para la partición P = P1×P2 = {x0, x 1,, x 2, x 3, …x
n}× {y0,y 1, y 2, y 3,…, y n}

3.3.-La función z=f(x,y) es integrable sobre una región R del plano XY, si “f” es acotada sobre R, y si
existe un numero I que satisface:

Lf (P) ≤ I ≤ Uf (P) ; I= 
R
f ( x, y)dxdy ; Para toda partición P de R.

3.4.- Para cada una de las subregiones Rij de una partición P de la región R, si f ( xi* , y*j ) es un
punto en Rij cualquiera sea la elección de f ( xi* , y*j ) en Rij , se tiene que:

n m
L f ( P)   f ( xi* , y*j )xi y j  U f ( P)
i 1 j 1

También vale decir:

n m

 f ( x, y)dxdy  f ( xi* , y*j )xi y j  U f  P   L f  P 

R
i 1 j 1

4.-EJERCICIOS

1.- Calcular:

R
7dxdy 

Dónde: R= {(x,y) / x ϵ [1,8], y ϵ [1,6]}

67
 Calculo II 2012

Solución:

 Para este caso:

Si (x,y) ϵ R, entonces f(x,y) = 7
 Además considérese la partición P1 = {x0=1, x 1,, x 2, x 3, …x n=8} en [1,8], la partición P2 =
{y0=1,y 1, y 2, y 3,…, y n=6} en [1,6] y la partición arbitraria P=P1×P2 en la región R.

 Así mismo en cada sub-región Rij, f(x,y)=7. Por lo que:

Mij = mij = 7

 Luego la suma superior es igual a la suma inferior:

n m m n
U f ( P)  L f ( P)   7xi y j  7( y j )( xi )  7(6  1)(8  1)  245
i 1 j 1 i 1 j 1

 Finalmente por propiedad:

Uf (P) ≤  R
f ( x, y)dxdy ≤ Lf(P)

245 ≤  R
f ( x, y)dxdy ≤ 245

 La integral es igual a:

 7dxdy =245
R

2.-Evaluar :

 R
f ( x, y)dxdy

68
 Calculo II 2012

Si f(x,y) = 8 - xy esta acotado en el rectangulo R de coordenadas (1;1), (1;4), (5;1) y (5;4).

-Solucion :

ym=4

yj ∆Aij
∆yj
yj-1

y0=1
∆xi

x0=1x0=1
. . .. . . xi-1 xi-1 xi . . .xix. n.=5
. xn=5 X

 Con las particiones [1,5] y [1,4] :

n m

 (8  xy)dxdy  lim  f ( xi* , y*j )Aij

R n 
i 1 j 1
m

P1 : xi = 1 + 4i/n , ∆xi = 4/n , i=0 ; 1 ; 2 ; 3… ; n

P2 : yj = 1 + 3j/m , ∆yj =3/m , j=0 ; 1 ; 2 … ; m

Apartir de P1 Y P2 obtenemos P = P1 × P2 ; entonces :

4i 3j
f ( xi* y*j )  f ( xi y j )  8  xi y j  8  (1  )(1  )
n m

Luego :

69
 Calculo II 2012

m m
4i 3j 3 4i 3 m(m  1) 3 4 9(m  1)
  f ( x,*i y*j )y j  [8  (1 
j 1 j 1 n
)(1  )]  [8m  (1  )(m  
m m n m 2
)]  24  (1  )(3 
m n 2m
)

9(m  1) 4 n(n  1) 4 9 8(n  1) 9 1

 [24n  (3  )(n   )]  96  [(3  4.5  )(4  )]  96  [(7.5  )(4  8 )]
2m n 2 n 2m n 2m n

Finalmente, como m y n tienden al infinito :

 R
(8  xy)dxdy =96-(7.5+0)(4+8+0) = 96-90 = 6

3.- Vamos a ver un ejemplo

∫∫

Donde R es el rectángulo 1 ,

Solución Con , , , como partición arbitraria de [1,4] y , , , como

partición arbitraria de [1,3] tenemos

={( , }

Como posición arbitraria de R. En cada rectángulo : , , la función

F(x,y)=x+y-2

Alcanza un máximo y un mínimo Luego:

∑∑

∑∑

Para cada índice i, j:

( )

Esto significa

∑ ∑ ( ) (P)

70
 Calculo II 2012

La suma doble en este conjunto de desigualdades se puede escribir como

∑ ∑ +∑ ∑ -∑ ∑ ,

La primera suma doble se reduce a

∑ ∑ = ∑ ( ) ∑ )=

La segunda suma doble se reduce a,

∑ ∑ = ∑ ∑ )=

La suma de estos números 15+12-12=15 satisface la desigualdad

(P) (P) para P arbitraria

Luego

∫∫

71
 Calculo II 2012

CALCULO DE INTEGRALES POR INTERPOLACIÓN

1.- INTRODUCCIÓN:
Sabemos que para desarrollar una integral debemos calcular la primitiva de la función pero no
siempre se puede hacer ello, en otros cosos no conocemos la función pero si tenemos datos en
forma de tablas o gráficas. En estos casos necesitamos un método para poder resolver estos
problemas, en este caso estudiaremos a Lagrange y posteriormente como calcular una integral con
este Lagrange.

2.-INTERPOLACIÓN DE LAGRANGE:

Sea ⌊ , ⌋ y {x0, x1,….., xn} n+1 puntos distintos del intervalo [a, b], que llamaremos nodos,
entonces existe un único solo polinomio de grado “n” a lo más Pn (x) tal que:

Pn ( xi )  f ( xi ) Para i  0,1,...., n

Este polinomio se le conoce como “enésimo polinomio interpolante de Lagrange”

n
Pn ( x)   f ( xi ) Li ( x) , i 0,1, 2,..., n
i 0

De donde:
n x  xj
Li ( x)  
i o xi  x j
i j

EJEMPLO 1:

Si queremos utilizar los números (o nodos), x0 = 2, x1 = 2.5 y x2 = 4 para obtener el segundo

1
polinomio interpolante para f ( x) 
x
Primero procedemos a determinar los coeficientes de L0 ( x), L1 ( x), L2 ( x) :
( x  2.5)( x  4)
L0 ( x)   x( x  6.5)  10
(2  2.5)(2  4)
( x  2)( x  4) x(24  4 x)  32
L1 ( x)  
(2.5  2)(2.5  4) 3
( x  2)( x  2.5) x( x  4.5)  5
L2 ( x)  
(4  2)(4  2.5) 3
Sabiendo que:
f ( x0 )  f (2)  0.5 , f ( x1 )  f (2.5)  0.4 , f ( x2 )  f (4)  0.25

72
 Calculo II 2012

Entonces tendríamos:
n
P( x)   f ( xi ) Li ( x)
i 0

x(24  4 x)  32 x( x  4.5)  5
P( x)  0.5( x( x  6.5)  10)  0.4  0.25
3 3

Por lo tanto este sería nuestro polinomio interpolante.

2.1.-ERROR DE INTERPOLACION:

Como ya tenemos nuestro polinomio, nos gustaría calcular el valor de la función en cualquier
punto del intervalo [a,b], si el punto que escogemos coincide con alguno de los nodos {x0, x1,…..,
xn} entonces:
Pn ( xi )  f ( xi )
Sin embargo sí x   a, b y es distinto a los nodos de interpolación entonces:
Pn ( xi )  f ( xi )
Es ahí donde se produce un error al que llamaremos error de interpolación denotado:

f ( n1) ( x ) n
En ( x)  f ( x)  Pn ( x) 
(n  1)! i 0
 ( x  xi )

La expresión nos permite obtener la cota, para cada x   a, b

 (x  x ) i

f ( x)  Pn ( x)  i 0
max f ( n 1) ( )
(n  1)! a  x b

EJEMPLO 2:

Suponga que debe preparar una tabla de la función f ( x)  e x , para x   0,1 . Suponga, además,
el número de cifras decimales de cada entrada o valor es de d  8 y que h , el tamaño del paso es
la diferencia entre los valores adyacentes de x . ¿Cuál debe ser el valor de h para que la
interpolación lineal (es decir, el polinomio de grado 1 de Lagrange) arroje un valor absoluto a lo
6
máximo de 10 ?

Sean x   0,1 , x0 , x1 ,...., xn los números en los que se evalúa f y suponga que j satisface
x j  x  x j 1 .

73
 Calculo II 2012

f (2) ( )
f ( x)  P( x)  ( x  x j )( x  x j 1 )
2!

f (2) ( )
f ( x)  P( x)  ( x  x j )( x  x j 1 )
2!

Por ser h el tamaño del paso, se deduce que x j  jh , x j 1  ( j  1)h , entonces

f (2) ( )
f ( x)  P( x)  ( x  jh)( x  ( j  1)h)
2!

Por tanto,

1
f ( x)  Pn ( x)  max e max ( x  jh)( x  ( j  1)h)
2 0 1 x j  x x j1

1
f ( x)  Pn ( x)  e max ( x  jh)( x  ( j  1)h)
2 x j  x x j1

Al considerar g ( x)  ( x  jh)( x  ( j  1)h) para jh  x  ( j  1)h y al aplicar los métodos del

cálculo, encontramos que

max ( x  jh)( x  ( j  1)h)

x j  x  x j 1

1
max g ( x)  g (( j  )h)
x j  x  x j 1 2

h2
max g ( x) 
x j  x  x j 1 4

Entonces el error de la interpolación lineal estaría acotado de la siguiente forma:

2
h
f ( x)  Pn ( x)  e
8

h  1.72 103

(1  0)
Ya que n  debe ser un entero, una elección lógica del tamaño del paso es h  0.001 .
h

74
 Calculo II 2012

2.2.-CALCULO DE UNA INTEGRAL USANDO LAGRANGE:

Como se mencionó anteriormente hay veces es necesario evaluar una integral definida de una
función que no tiene antiderivada explicita, o cuya antiderivada no es fácil de obtener. El método
b
básico con el que se aproxima  f ( x)dx recibe el nombre de cuadratura numérica y emplea una
a

suma de tipo:

 a f (x )
i 0
i i

b
Para poder aproximar  f ( x)dx el método de la cuadratura que se presenta se basa en los
a

polinomios interpolantes.

Primero seleccionamos un conjunto de nodos distintos  x0 ,...., xn  del intervalo  a, b . Después

integramos el polinomio interpolante de Lagrange.

n
Pn ( x)   f ( xi ) Li ( x)
i 0

Y su termino de error en  a, b para obtener

b b n b n f ( n1) (( x ) )
 f ( x)dx    f ( xi ) Li ( x)dx    ( x  xi ) dx
a a i 0 a i 0
(n  1)!

b n b
n
1
 f ( x)dx   ai f ( xi )    ( x  xi ) f ( n 1) (( x ) )dx
a i o (n  1)! a i 0

Donde  ( x ) se encuentra en  a, b para cada x

b
ai   Li ( x)dx Para cada i  0,1,....., n
a

Por lo tanto la fórmula de la cuadratura es:

b n

 f ( x)dx   ai f ( xi )
a i o

75
 Calculo II 2012

b
n
1
Con un error: E( f )   ( x  xi ) f (n1) (( x) )dx
(n  1)! a i 0

b
EJEMPLO 3: se desea aproximar  f ( x)dx
a

Sean x0  a , x1  b , h  b  a , usaremos el polinomio lineal de Lagrange:

( x  x1 ) ( x  x0 )
P1 ( x)  f ( x0 )  f ( x1 )
( x0  x1 ) ( x1  x0 )

Entonces:

 ( x  x1 ) ( x  x0 ) 
b 1x x
1 1

a
f ( x)dx   
x0 
( x0  x1 )
f ( x0 ) 
( x1  x0 )
f ( x1 ) dx   f ( ( x))( x  x0 )( x  x1 )dx
 2 x0

Dado que ( x  x0 )( x  x1 ) no cambia de signo en  x0 , x1  , podemos aplicar el teorema del valor

promedio ponderado de las integrales al término de error para poder obtener:

x1 x1

 f ( ( x))( x  x )( x  x )dx  f ( )  ( x  x )( x  x )dx

x0
0 1
x0
0 1

x1
 x3 x  x1 2 
 f ( )   x  x0 x1 x 
3 2  x0

h3
  f ( )
6

En consecuencia:

x1
b
 ( x  x1 )2 ( x  x0 )2  h3

a
f ( x)dx  
 2( x0  x1 )
f ( x0 ) 
2( x1  x0 )
f ( x1 )  
 x0 12
f ( )

( x1  x0 )
b
h3
 f ( x)dx   f ( x0 )  f ( x1 )  f ( )
a
2 12

Pero x1  x0  h

b
h h3
 f ( x)dx   f ( x0 )  f ( x1 )  f ( )
a
2 12

76
 Calculo II 2012

3.-INTERPOLACIÓN DE HERMITE

3.1.-Definición

Sea x0,x1,...,xn, n+1 números distintos en [a,b] y mi un entero no negativo asociado a xi para
i=0,1,...,n. Supónganse que f [a,b] y que m=max0<i<n mi. El polinomio oscilante que aproxima f
es el polinomio P(x) de menor grado tal que

d k p  xi  d k f ( xi )
 para cada i=0,1,..., n y k=0,1,...,mi .
dx k dx k

Nótese que, cuando mi= 1 para cada i=0,1,..., n, se produce una clase de polinomios denominados
polinomios de Hermite. En una función dada f, estos últimos concuerdan con f en x0,x1,...,xn.
Además, como sus primeras derivadas concuerdan con las de f; tendrán la misma “forma” que la
función en (xi, f (xi)), en el sentido de que las líneas tangentes del polinomio coinciden con las de la
función.

3.2.-Teorema

Si f [a,b] y si x0,x1,...,xn [a,b] son distintos, el polinomio único de menor grado que
concuerda con f y f’ en x0,x1,...,xn es el polinomio de Hermite de grado a lo más 2n+1 que está
dado por:

n n
H 2 n1 ( x)   f ( x j ) H n. j ( x)   f ( x j ) Hˆ n, j ( x)
j 0 j 0

Donde:

H n, j  x   1  2  x  x j  Ln, j  x j  L2n, j  x 

Hˆ n, j  x   ( x  x j ) L2n, j  x 

Más aún, si f [a,b], entonces para x [a,b]

( x  x0 )2 ( x  xn )2 (2 n 2)
f ( x)  H 2 n1 ( x)  f ( ) ; para algún con a< <b.
(2n  2)!

77
 Calculo II 2012

Ejemplos:

1) Determinar el polinomio interpolador de Hermite de la función f(x)=cos(x) en el soporte

{0, /2, } para los enteros {2,0,1}

Solución:
n
Grado del polinomio interpolador: m  n  m
i 0
i  2  2  0 1  5

p5 ( x)  m0  m1.x  m2 .x 2  m3 .x3  m4 .x 4
p5 (0)  cos(0)  m0  1
p5 (0)   sen(0)  m1  0
p5 (0)   cos(0)  2.m2  1
p5 ( / 2)  cos( / 2)  m0  m1.( / 2)  m2 .( / 2) 2  m3 .( / 2) 3  m4 .( / 2) 4  m5 .( / 2) 5  0
p5 ( )  cos( )  m0  m1.( )  m2 .( ) 2  m3 .( )3  m4 .( ) 4  m5 .( )5  1
p5 ( )   sen( )  m1  2.m2 .( )  3.m3 .( ) 2  4.m4 .( )3  5.m5 .( ) 4  0

Por sistemas de ecuaciones se tiene:

 1 
 
 0 
 0
m
   1 2 
 m1   
 m2   2.   10 
2

  3 
 3 
m
2 
 m4   5 . 12   
   2  4 
 m5   3 
   12 
 5 
Se tiene:
 2  10 3 5 12   2 4  2  12 5
p5 ( x)  1  (1/ 2).x 2  2( ).x  .( ).x  ( ).x
3 2 4 5

2) Utilizar el método de Hermite para hallar un polinomio P(x) de grado 2 que satisfaga:
p 1  0, p 1  7, p  2   10 .

Solución:

78
 Calculo II 2012

Como existe la derivada del polinomio p(x), quiere decir que el método a utilizar es el de
Hermite (en el caso de que no nos dijeran el método a utilizar), entonces la tabla quedaría
de la forma

X P(x)
10
10 7
2 10 p[1,2]=10 p[1,1,2]=3

p [1,2]= 10/1
p [1, 1,2]= (10-7)/1 = 3
El polinomio de interpolación quedaría expresado de la forma
p  x   7  x  1  3  x  1
2

p  x   3x 2  x – 4 Satisface las condiciones de

p 1  0, p 1  7, p  2   10 .

3) Sea la función: f(x)= , calcúlese el polinomio de interpolación de Hermite de f en x0=1 y

x1=2.

Solución:
Construimos la tabla de diferencias divididas que permite incluir los datos de interpolación
f(1), f’(1), f(2), f’(2), trabajando con dos dígitos decimales:

D.D.1 D.D.2 D.D.3

x0=1 f(1)=0 f’(1)=1 f[1,1,2]=-0.31 f[1,1,2,2]=0.12
x0=1 f(1)=0 f[1,2]=0.69 f[1,2,2]=-0.19
x1=2 f(2)=0.69 f*2,2+= f’(2)=0.5
x1=2 f(2)=0.69

Sea p(x) el polinomio de interpolación buscado, utilizando la formula de Newton

generalizada

p  x   f 1  f 1,1  x – 1  f 1,1, 2  x – 1  f 1,1, 2, 2   x – 1  x – 2 

2 2

  x  1 – 0.31 x – 1  0.12  x – 1  x – 2 
2 2

4) Supongamos que queremos encontrar el polinomio p(x) que satisfice p(x0)=c00, p(xi)=c10 y
p’(x0)=c01.

79
 Calculo II 2012

Solución:
Puesto que las condiciones son de interpolación de Hermite, con m + 1=2 + 1=3 (2 es el
número de condiciones sobre el nodo x0 y 1 es el número de condiciones sobre el nodo x1),
entonces existe un único polinomio p(x) ∏ que satisface las condiciones dadas.
Consideremos la siguiente tabla de diferencias dividas extensivas:
zi D.D.0 p(zi) D.D.1 D.D.2
z0=x0 c00=p(z0) p[z0,z1]= c01 p[z0,z1 ,z2]= c02
z1=x0 c00=p(z1) p[z1,z2]
z2=x1 c10=p(z2)

Nota: En esta tabla se repite el nodo x tantas veces como condiciones hayan impuestas
sobre ese nodo. Cada uno de los cálculos que aparecen e la tabla se realiza usando la
definición de diferencia dividida con o sin repetición, según sea el caso.
p  z0. z1   p  x0, x0   p  x0   c01 (Definición)
p( x1 )  p( x0 ) c10  c00
p  z1 , z2   p  x0 , x1   
x1  x0 x1  x0
c10  c00
 c01
p[ z1 , z2 ]  p[ z0 , z1 ] x1  x0 c c c
p[ z0 , z1 , z2 ]    10 002  01  c02
z 2  z0 x1  x0 ( x1  x0 ) x1  x 0
El polinomio
p( x)  p[ z0 ]  p[ z0 , z1 ]( x  z0 )  p[ z0 , z1, z2 ]( x  z0 )( x  z1 )
Es decir,
p( x)  c00  c01 ( x  x0 )  c02 ( x  x0 )2
Satisface las condiciones dadas. En efecto:
p( x)  c00
p( x1 )  c00  c01 ( x1  x0 )  c02 ( x1  x0 )
p( x1 )  c00  c01 ( x1  x0 )  [c10  c00  c01 ( x1  x0 )]
p( x )  c10
p( x)  c01  2c02 ( x  x0 ), p( x0 )  c01
En conclusión p(x) es el polinomio de interpolación de Hermite para las condiciones dadas.

5) Use el método de diferencias divididas extensivas de Newton para obtener el polinomio

de Hermite que interpola la siguiente tabla:
xi f(xi) f’(xi) f’’(xi)
0 2 3
1 6 7 8

80
 Calculo II 2012

Solución:
Como k0=2 (número de condiciones sobre el nodo x0=1) y k1=3 (número de condiciones
sobre el nodo x1=2), entonces en la tabla se repite el nodo x0=1 dos veces y el nodo x1=2
tres veces. Así que se construye la siguiente tabla de diferencias divididas:

zi D.D.0 D.D.1 D.D.2 D.D.3 D.D.4

z0=1 2 f[1,1]=3 f[1,1,2]=1 f[1,1,2,2]=1 f[1,1,2,2,2]=-1
z1=1 2 f[1,2]=4 f[1,2,2]=3 f[1,2,2,2]=1
z2=2 6 f[2,2]=7 f[2,2,2]=4
z3=2 6 f[2,2]=7
z4=2 6

En esta tabla los cálculos ser realizan como se indican a continuación:

f (1) f (2)
f [ z0 , z1 ]  f [1,1]   3; f [ z2 , z3 ]  f [2, 2]  7
1! 1!
f (2)
f [ z, z ]  f [2, 2]  7
1!
f (2) 8
f [ z2 , z3 , z4 ]  f [2, 2, 2]   4
2! 2
Luego el polinomio

p( x)  2  3( x  z0 )  1( x  z0 )( x  z1 )  2( x  z0 )( x  z1 )( x  z2 )  (1)( x  z0 )( x  z1 )( x  z2 )( x  z3 )
p( x)  2  3( x  1)( x  1) 2  2( x  1) 2 ( x  2)  ( x  1) 2 ( x  2) 2

Satisface las condiciones dadas. En efecto:

p(1)  2; p(2)  2  3(2  1)  (2  1) 2  2(2  1) 2 (2  2)  (2  1) 2 (2  2) 2  5  1  6

p( x)  3  2( x  1)  2( x  1) 2  4( x  1)( x  2)  2( x  1)( x  2) 2  2( x  1) 2 ( x  2)
p(1)3 ; p(2)  3  2(2  1)  2(2  1) 2  3  2  2  7
p( x)  2  4( x  1)  4( x  1)  4( x  2)  2( x  2) 2  4( x  1)( x  2)  4( x  1)( x  2)  2( x  1) 2
p(2)  2  8(2  1)  2(2  1) 2  8

81
 Calculo II 2012

4.-INTERPOLACIÓN DE TRAZADORES CÚBICOS

4.1.-Definición
Un trazador cúbico S es una función a trozos que interpola a f en los n+1 puntos (x 0,y0), (x1,y1),
(x2,y2),..., (xn,yn) con a=x0<x1<...<xn=b. S es definida de la siguientes manera:

 S0 si x  [ x0 , x1 ]
S x  [ x1 , x2 ]
 si
S ( x)   1

 Sn 1 x  [ xn 1 , xn ]
Dónde:
1. Si ( x)  ai  bi ( x  xi )  ci ( x  xi )2  di ( x  xi )3 para i=0,1,...,n+1
2. S ( xi )  yi , i  0,1, , n Para efectos prácticos, S j ( x j )  y j , j  0,1, , n 1 y
Sn1 ( xn1 )  yn1 , Sn1 ( xn )  yn . El formulario asegura que S j ( x j  1)  y j  1.
Si ( xi  1)  Si 1 ( xi  1)
3. para i=0,1,...,n-2

4. Si ( xi  1)  Si 1 ( xi  1) para i=0,1,...,n-2

5. Si ( xi  1)  Si 1 ( xi  1) para i=0,1,...,n-2
6. Se satisface una de las dos condiciones que siguen
 S ( x0 )  S ( xn )  0 (frontera libre o natural)
 S ( xi  1)  f ( x0 ) y S ( xn )  f ( xn ) (frontera sujeta)

De acuerdo a la condición 1, Si ( xi )  yi  ai  yi .

Haciendo algunas manipulaciones algebraicas en el sistema, se obtiene

ci 1  ci y  yi hi (2ci  ci 1 )
di  y bi  i 1 
3hi hi 3

La condición de frontera natural hace que c0  cn  0 .

Ahora todo depende del cálculo de los ci’s. Estos se calculan resolviendo el sistema (n+1)x(n+1).

82
 Calculo II 2012

 0 
1 0 0 0   c0   
   c   3( f [ x2 , x1 ]  f [ x1 , x0 ]) 
 h0 2(h0  h1 ) h 0  1 
 3( f [ x3 , x2 ]  f [ x2 , x1 ]) 
   
    
0 hn 3 2(hn 3  hn  2 ) hn  2   cn 1 
 3( f [ xn , xn 1 ]  f [ xn 1 , xn  2 ]) 
0 1  
 0  cn   
 0 

Si: f [ xi , x j ]  ( yi  y j ) / ( xi  x j )

4.2.-EJERCICIOS

1) Determinar el trazador cúbico (frontera libre) para la siguiente tabla:

xi yi=cos(3xi2).ln(xi3+1)
x0=0 0
x1=0.75 -0.04099838
x2=1.5 1.31799

Solución:
El trazador es:
S0 ( x)  a0  b0 ( x  x0 )  c0 ( x  x0 ) 2  d0 ( x  x0 )3 si x [ x0 , x1 ]
S ( x)  
 S1 ( x)  a1  b1 ( x  x1 )  c1 ( x  x1 )  d1 ( x  x1 ) x  [ x1 , x2 ]
2 3
si
Hay que determinar los coeficientes de S0 y S1 resolviendo el sistema 8x8,
 S0 ( x0 )  y0
 S (x )  y  a0  0
 1 1 1  a1  0.0409838
 S1 ( x2 )  y2 
 a1  0.75b1  0.5625c1  0.421875d1  1.31799
S
 0 1 ( x )  S ( x ) 
 a0  0.75b0  0.5625c0  0.421875d 0  a1
1 1

 S0 ( x1 )  S1 ( x1 )  
  b0  1.5c0  1.6875d 0  b1
 S0 ( x1 )  S1 ( x1 )  2c0  4.5d 0  2c1
 
 S0 ( x0 )  0  2c0  0
  2c1  4.5d1  0
 S1 ( x2 )  0 

La solución de este sistema es:

a0  0, b0  0.521299, c0  0, d0  0.829607, a1  0.0409838, b1  0.
878663, c1  1.86662, y d1  0.829607.

83
 Calculo II 2012

Es decir:
 S0  0.521299 x  0.829607 x3 si x  [0,0.75]
S ( x)  
S1 ( x)  0.0409838  0.878663( x  0.75)  1.86662( x  0.75)  0.829607( x  0.75) si x  [0.75,1.5]
2 3

2) Determinar el trazador cúbico (frontera libre) para la siguiente tabla:

xi yi=xi4-4xi3
x0=-2 48
x1=-1 5
x2= 0 0
x0= 1 -3
x1= 2 -16
Solución:

Observe que la función x 4  4 x3 tiene un punto de inflexión en x=0

El trazador es:
 S0 ( x)  9.85714( x  2)3  52.8571( x  2)  48 si x [ 2, 1]

 S ( x)  11.2857( x  1)  29.5714( x  1)  23.2857( x  1)  5 x [ 1, 0]
3 2
si
S ( x)   1
 S2 ( x)  0.714286 x3  4.28571x 2  2 x si x  [0,1]
 S3 ( x)  2.14286( x  1)3  6.42857( x  1) 2  8.71429( x  1)  3 si x  [1, 2]

3) Interpolar los siguientes datos mediante un trazador cúbico:

x 2 3 5
y -1 2 -7

Solución:
Definimos un polinomio cúbico en cada uno de los intervalos que se forman:
S0 ( x)  a0  b0 x  c0 x 2  d0 x3 si x  [2,3]
S  x  
 S1 ( x)  a1  b1 x  c1 x  d1 x x  [3,5]
2 3
si
A continuación, hacemos que se cumpla la condición de que el trazador debe pasar por los
puntos dados en la tabla. Así, tenemos que:
S (2)  1  a0  2b0  4c0  8d 0  1
S (3)  2  a0  3b0  9c0  27d0  2
S (5)  7  a1  5b1  25c1  125d1  7
Ahora calculamos la primera derivada de S(x):
b  2c0 x  3d0 x 2 si x  [2,3]
S ( x)   0
 b1  2c1 x  3d1 x x  [3,5]
2
si

84
 Calculo II 2012

Al igual que en el caso de los trazadores cuadráticos, se presentan ecuaciones que pueden
presentar discontinuidad en los cambios de intervalo; las posibles discontinuidades son los
puntos donde se cambia de intervalo, en este caso x=3. Para evitar esta discontinuidad,
evaluamos x=3 en los dos polinomios e igualamos:
b0  6c0  27a0  b1  6c1  27a1
Análogamente procedemos con la segunda derivada:
2c  6d0 x si x  [2,3]
S ( x)   0
 2c1  6d1 x si x  [3,5]
Para lograr que S’’(x) sea continua:
18d0  2c0  18d1  2c1
En este punto contamos con 6 ecuaciones y 8incognitas, por lo tanto tenemos 2 grados de
libertad; en general; se agregan las siguientes 2 condiciones:
S ( x0 )  0
S ( xn )  0
De lo cual vamos a obtener:
S (2)  0  2c0  12d0
S (5)  0  2c1  30d1
Con la cual, hemos completado un juego de 8 ecuaciones y 8 incógnitas, el cual es el
siguiente:
a0  2b0  4c0  8d 0  1
a0  3b0  9c0  27d 0  2
a1  3b1  9c1  27d1  2
a1  5b1  25c1  125d1  7
a0  3b0  9c0  27d 0  a1  3b1  9c1  27d1
2c0  18d 0  2c1  18d1
2c0  12d 0  0
2c1  30d1  0
Calculando las soluciones mediante matrices obtenemos:
a0  1.25 a1  0.625
b0  7.5 b1  9.375
c0  10.7 c1  39.875
d0  0.5 d1  50.125
Sustituyendo estos valores en nuestra función inicial, vemos que el trazador cúbico para la
tabla de datos dada, queda definida como sigue:
 S0 ( x)0.5  10.75 x  7.5 x 2  125 x3 si x  [2,3]
S ( x) 
S1 ( x)  50.125  39.875 x  9.375x  0.625x x  [3,5]
2 3
si

85
 Calculo II 2012

4) Interpolar los siguientes datos utilizando trazadores cúbicos:

x -1 1 2 4
y -1 1 5 -2

Solución:
Nuevamente, definimos un polinomio cúbico en cada uno de los intervalos:
 S0 ( x)  a0  b0 x  c0 x 2  d 0 x 3 si x  [1,1]

S ( x)  S1 ( x)  a1  b1 x  c1 x 2  d1 x3 si x  [1, 2]
 S ( x)  a  b x  c x 2  d x 3 si x  [2, 4]
 2 2 2 2 2

Después, hacemos que el trazador pase por los puntos dados en la tabla. Así tenemos:
S (1)  1 implica que, a0  b0  c0  d0  1
a0  b0  c0  d 0  1
S (1)  1 Implica que,
a1  b1  c1  d1  1

a1  2b1  4c1  8d1  5

S (2)  5 Implica que,
a2  2b2  4c2  8d 2  5

S (4)  2 Implica que, a2  4b2  16c2  64d2  2

Enseguida calculamos la primera derivada:
b0  2c0 x  3d 0 x 2 si x  [1,1]

S ( x)   b1  2c1 x  3d1 x 2 si x  [1, 2]
b  2c x  3d x 2 si x  [2, 4]
 2 2 2

Vemos entonces, que las posibles discontinuidades de S’(x) son x=1 y x=2 . Por lo tanto,
para hacer que S’(x) sea continua, igualamos las ecuaciones correspondientes en ambos
valores:
b0  2c0  3d0  b1  2c1  3d1
b1  4c1  12d1  b2  4c2  12d 2
Ahora procedemos a calcular la segunda derivada:
2c0  6d 0 x si x  [1,1]

S ( x)   2c1  6d1 x si x  [1, 2]
2c  6d x si x  [2, 4]
 2 2

Nuevamente, las posibles discontinuidades son x=1 y x=2 . Por tanto, para que S’’(x) sea
continua, se igualan las ecuaciones en ambos valores:
2c0  6d0  2c1  6d1
2c1  12d1  2c2  12d 2
Finalmente, se agregan las condiciones de que la segunda derivada se anule en los puntos
inicial y final de la tabla. En este caso,

86
 Calculo II 2012

S (1)  0  2c0  6d0  0

S (4)  0  2c2  24d 2  0
Con esto tenemos un juego de 12 ecuaciones con 12 incógnitas:
a0  b0  c0  d 0  1
a0  b0  c0  d 0  1
a1  b1  c1  d1  1
a1  2b1  4c1  8d1  5
a2  2b2  4c2  8d 2  5
a2  4b2  16c2  64d 2  2
b0  2c0  3d 0  b1  2c1  3d1
b1  4c1  12d1  b2  4c2  12d 2
c0  3d 0  c1  3d1
c1  6d1  c2  6d 2
c0  3d 0  0
c2  12d 2  0
Calculando las soluciones mediante matrices obtenemos:
153 89 153 51
a b c d
40 140 140 40
48 473 297 21
a b c d 
35 70 35 10
732 1867 288 24
a b c d
35 70 35 35
Por lo tanto, el trazador cúbico es:
 153 89 153 2 51 3
 S0 ( x)   40  140 x  140 x  140 x si x  [1,1]

 48 473 297 2 21 3
S ( x)  S1 ( x)   x x  x si x  [1, 2]
 35 70 35 10
 732 1867 288 2 24 3
 S2 ( x)   35  70 x  35 x  35 x si x  [2, 4]


5) Determinar el trazador cúbico (frontera natural) para la siguiente tabla:

xi yi=cos(3xi2).ln(xi3+1)
x0=0 0
x1=0.5 0.0861805
x2=1 -0.686211
x3=1.5 1.31799

Solución:
El trazador es:

87
 Calculo II 2012

 S0 ( x)  a0  b0 ( x  x0 )  c0 ( x  x0 ) 2  d0 ( x  x0 )3 si x  [ x0 , x1 ]

S ( x)  S1 ( x)  a1  b1 ( x  x1 )  c1 ( x  x1 ) 2  d1 ( x  x1 )3 si x  [ x1 , x2 ]
 S ( x)  a  b ( x  x )  c ( x  x ) 2  d ( x  x ) 3 si x  [ x2 , x3 ]
 2 2 2 2 2 2 2 2

Hay que determinar los coeficientes de S0 y S2 .Iniciamos calculando los ci’s. Resolvemos el
sistema 4x4. Recordemos que hi  xi 1  xi .
 0 
 
1 0 0 0   c0  y2  y1 y1  y0 
   c   3( x  x  )
x1  x0 
 h0 2(h0  h1 ) h1 0  1 

2 1

0 h1 2(h0  h1 ) h2   c2   y3  y2 y2  y1 
    3(  )
0 0 0 1  c3   x3  x2 x2  x1 
 
 0 
 1 0 0 0   c0   0 
  c   
 0.5 2 0.51 0   1    5.15143 
 0 0.5 2 0.5   c2   16.6596 
    
 0 0 0 1   c3   0 

Obtenemos c0  0, c1  4.96871, c2  9.57196 y por convenio el comodín c3  0 .

Ahora podemos obtener el resto de coeficientes:
ai  yi
b0  1.00048 b1  1.48387 b2  0.8177508
d 2  3.31247 d 2  9.69378 d 2  6.38131
Finalmente el trazador cúbico es:
 S0 ( x)  1.00048 x  3.31247 x3 si x  [0, 0.5]

S ( x)  S1 ( x)  0.0861805  1.48387( x  0.5)  4.96871( x  0.5) 2  9.69378( x  0.5)3 si x  [0.5,1]
 S ( x)  0.686211  0.8177508( x  1)  9.57196( x  1) 2  6.38131( x  1)3 si x  [1,1.5]
 2

88
 Calculo II 2012

BIBLIOGRAFÍA

Autor – editorial – nombre del libro – paginas

Espinoza Ramos Eduardo, Editorial: EDUARDO ESPINOZA RAMOS, Análisis Matemático II.
Pág. 269-279.
Carreño Campos, Ximena; Cruz Schmidt, Ximena., Editorial: ARRAYAN EDITORES S.A.,
Algebra. Pág. 473-480.

Granville, Smith, Longley, Editorial: UTHEA, Calculo diferencial e integral, pág. 309 -312.
Autor: Eduardo Espinoza Ramos, Editorial: EDUKPERU, Análisis Matemático II, pág. 323,
324,343-345.

Espinoza ramos, Eduardo – Análisis matemático II cuarta edición (pág. 280-314).

Jorge Gárate, Máximo – Cálculo integral (pág. 6-9).
Larson, Hostetler, Edwards – Cálculo y geometría analítica sexta edición (pág. 291-305).

Espinoza Ramos Eduardo, Editorial: EDUARDO ESPINOZA RAMOS, Análisis Matemático ll,
pág. 630-633
Juan de Burgos, 2º Edición, Cálculo infinitesimal de una variable, pág. 369.

Máximo Mitac-Luis Toro, Editorial San Marcos, Tópicos de cálculo volumen 2. pág. 178-
182
Louis Leithold. Editorial Iberoamérica, El cálculo con geometría analítica segunda edición
pág. 549-555.

Leithold, Editorial Iberoamérica, el cálculo 7ma edición, pág. 639 y 738

Integrales Múltiples, Armando Venero, Editorial Gemar, Matematicas III Transformaciones,
Pag. 73-95.
Heraldo Gonzales Serrano, Universidad de Santiago de Chile, Matemática General.
Tom Apostol, Editorial Reverté, Calculus Multivariable Vol. II Pág. 430-455.
Erwin Kreyszig, Grupo Noriega Editores, Matemáticas Avanzadas Para Ingeniera Vol. II Pág.
524-527.

Richard L. Burden, International Thomson Editores, Análisis Numérico, Pág. 105-115, 134-
157, 188-190 y 190-237.

89
 Calculo II 2012

INTEGRANTES CON SUS RESPECTIVOS TEMAS

Delegado: Reyna Córdova Marco Antonio

Yslache Galvan Miguel Ángel – Sumatorias Definición, propiedades, ejemplos.

Machaca Paitan Pedro Javier y Cancho Puse Erick Eduardo – Calculo de la integral definida
mediante sumatorias, definición, propiedades, ejemplos.

Sánchez Linares Johana Allisson y Yupanqui Rosas Anthony Juvenal – Calculo de Áreas por
sumatorias, definición, propiedades, Ejemplos

Salazar Leguía Adolfo Eduardo – Método del rectángulo.

Palomino Olivares Diego Oswaldo – Método del trapecio.

Sánchez Hervacio Jhonatan –Método de Simpson.

Mendoza Ttito Frank Jonathan – Integración por series de Taylor.

Peña Landeo Victor Daniel – Sumatorias dobles, definición, propiedades, ejemplos.

Bazan Pretell Antony Nilo y Licla Gutiérrez Luis Felipe – Calculo de integrales dobles por
sumatorias dobles.

Ruiz Rodríguez Omar y Reyna Córdova Marco Antonio – Calculo de integrales por
interpolación.

90