UNIVERSIDAD CENTROAMERICANA
JOS SIMEN CAAS
TAREA INTEGRACION NUMERICA
MATEMATICAS II
BILLY BRIAN HENRIQUEZ ALVAREZ 00061511
EDUARDO BALMORE MENJIVAR RIVERA 00085411
NESTOR JAVIER ALVARADO VENTURA 00095512
NOVIEMBRE 2013
ANTIGUO CUSCATLN
EL SALVADOR, C.A.
INDICE
PARTE I. ............................................................................................................................................... 3
Integracin Numrica...................................................................................................................... 3
Introduccin. ................................................................................................................................... 3
Aproximacin Trapezoidal. IT .......................................................................................................... 3
Aproximacin de Simpson. S2N ........................................................................................................ 5
Sumatoria de Riemann(Resea histrica de orgenes.) .................................................................. 6
PARTE II: ............................................................................................................................................ 10
PARTE III ............................................................................................................................................ 46
Densidad constante....................................................................................................................... 46
SECCIN P.................................................................................................................................. 46
SECCIN Q ................................................................................................................................. 46
EN SUPERFICIE TOTAL ............................................................................................................... 47
Densidad variable .......................................................................................................................... 47
SECCIN P.................................................................................................................................. 47
SECCIN Q ................................................................................................................................. 47
EN SUPERFICIE TOTAL ............................................................................................................... 48
Densidad constante....................................................................................................................... 48
MASAS ....................................................................................................................................... 49
SECCIN P.................................................................................................................................. 49
SECCIN Q ................................................................................................................................. 50
SECCIN R ................................................................................................................................. 50
SECCIN S .................................................................................................................................. 50
CONCLUSIONES ................................................................................................................................. 51
RECOMENDACIONES ......................................................................................................................... 51
REFERENCIAS. .................................................................................................................................... 51
PARTE I.
Integracin Numrica.
Introduccin.
La necesidad de aproximar numricamente el valor de una integral surge por dos motivos
fundamentalmente; la dificultad o imposibilidad en el clculo de una primitiva, la funcin a
integrar solo se conoce por una tabla de valores.
la dificultad o imposibilidad en el clculo de una primitiva, la funcin a integrar solo se
conoce por una tabla de valores.
Es una amplia gama de algoritmos para calcular el valor numrico de una integral definida,
muchas veces tambin se utiliza para resolver ecuaciones diferenciales. Cuadratura
numrica es ms o menos sinnimo de integracin numrica, especialmente si se aplica a
integrales de una dimensin, pero tambin se utiliza para integrales mltiples. La
integracin numrica es una aproximacin de a integral definida.
Aproximacin Trapezoidal. IT
La regla del trapecio es un mtodo de integracin numrica, y para obtener el valor de la
integral definida se hacen aproximaciones mediante trapecios, este mtodo se basa en la
aproximacin de la integral definida f(x) por medio de una funcin lnea (Leithold, 1973).
De la figura podemos deducir que
x=(b/a)/n. Si n es suficientemente
grande el (x seria suficientemente
pequeo), el rea de los trapecios
ser aproximadamente el rea
pedida.
Deduccin del mtodo.
Corresponde al caso donde n=1, es
3
decir:
datos:
( )
( )
x
y
donde F1(x) es un polinomio de grado 1 para los
A
f(a)
B
f(b)
Este polinomio de interpolacin puede expresarse como
Integrando se obtiene:
Necesitamos la suma de las sub reas de los trapecios donde x= (b-a)/n que ser el ancho
de las sub reas que comprenden la curva. Si llamamos a las ordenadas (i=0,1, 2,3,..., n) el
rea ser:
El rea comprendida entre X=a Y=b est dada por:
4
Si sustituimos (1) en la expresin (2) obtenemos:
Lo que se conoce como frmula trapezoidal que puede expresarse as:
(Leithold, 1973).
Aproximacin de Simpson. S2N
Este mtodo es ms exacto para obtener el valor de la integral, que el mtodo del trapecio,
para este mtodo utilizamos un polinomio de orden superior para conectar los puntos.
Este mtodo consiste en conectar grupos sucesivos de tres puntos, sobre la curva mediante
parbolas de segundo grado y se suman las reas bajo las parbolas. Para as obtener el rea
aproximada bajo la curva (Leithold, 1973).
Por ejemplo un polinomio de interpolacin cuadrtico. Sea p2(x) el polinomio de grado (a
lo ms) dos que interpola a f(x) en x=a, x=(a+b)/2, x=b. Este polinomio se puede escribir
como:
Tenemos ahora que
Pero con h=(b-a)/2 y u=x-a tenemos que
5
En forma similar se obtiene que
Tenemos pues que
Argumentando en forma similar a en mtodo del trapezoide, tenemos que si n es un entero
par; entonces
Usando la frmula (**) podemos aproximar
Ahora
Esta frmula se conoce como la regla (compuesta) de Simpson para aproximar a I (f)
(Leithold, 1973).
Sumatoria de Riemann(Resea histrica de orgenes.)
Riemann publica en 1854 su obra Ueber die Darstellbarkeit einer Function
durch eine trigonometrische Reihe para poder acceder al cargo de profesor auxiliar en la
Universidad de Gottingen. Se define por primera vez el concepto de integral de Riemann
y se inicia la teora de funciones de una variable real.
Este tipo de integral se define como la suma finita de reas de rectngulos. Estos se definen
en base a una particin de un intervalo [a; b] donde se quiere calcular la integral, como
altura de los rectngulos hay varias opciones, el mximo o mnimo de la funcin en el subintervalo o cogiendo el valor de la funcin en los extremos del mismo.
Entonces la integral S es igual a esta suma, si para cualquier particin con sub-intervalos
ms pequeos que
por lo tanto se cumple que:
Es decir al aumentar el nmero de rectngulos la suma de sus reas aproxima mejor el rea
de la integral.
La aplicacin de la integral definida de Riemann es principalmente el clculo de reas
planas delimitadas por una funcin y el eje x
El clculo de volmenes de revolucin, es decir un recinto limitado por las rectas
x = a, x = b y la grfica de la funcin de forma que:
Tambin se puede calcular la longitud de una curva:
1) Aproximacin por el punto de izquierda. IL
n-1
7
f(xk)x
n=0
= f(x0)x + f(x1)x + ... + f(xn-1)x
= [f(x0) + f(x1) + ... + f(xn-1)]x
Se define por:
El lado izquierdo coincide con el lado izquierdo de la grafica.
2) Aproximacin por el punto de Derecha. IR
Se puede asumir que el rea de toda la regin comprendida en el intervalo [a,b] se
aproxima con la suma de las reas de estos rectngulos.
n
f(xk)x
n=1
f(x1)x + f(x2)x + ... + f(xn)x
[f(x1) + f(x2) + ... + f(xn)]x
Se tiene que:
3) Aproximacin por el punto Medio. IM
Se toman como puntos de clculo los centros de los segmentos de la subdivisin regular.
La suma es :
9
PARTE II:
Con cada uno de los mtodos numricos (LI , RI , MI , TI , S2n) deber de resolver las siguientes
integrales:
1-
METODO DEL TRAPECIO
FUNCION =
f(x)= sqrt(tan(x))
Valor de pi =
3,141592654
Limite superior de la Integral (b) =
1,047197551
Limite inferior de la Integral (a) =
0,523598776
Valor de sub-intervalos (n) =
con n = 5
h = 0,104719755
n
f(x)
0,523598776
0,759835686
0,628318531
0,852374641
0,733038286
0,948896224
0,837758041
1,053856022
0,942477796
1,173193045
1,047197551
1,316074013
SUMATORIA
4,028319932
INTEGRAL
0,530539054
10
METODO DE SIMPSON
FUNCION
=
f(x)= sqrt(tan(x))
Valor de pi
=
3,141592654
Limite superior de la Integral (b) =
1,047197551
Limite inferior de la Integral (a) =
0,523598776
Valor de sub-intervalos (n) =
El metodo de aproximacion de Simpsom no acepta valores de subintervalos(n) impares
INTEGRACION NUMERICA A PARTIR
DE LA SUMA DE RIEMANN
FUNCION =
f(x)= sqrt(tan(x))
Valor de pi =
3,141592654
Limite superior de la Integral (b)
=
Limite inferior de la Integral (a)
=
Valor de sub-intervalos (n) =
1,04719755
0,52359878
5
con n = 5
Aproximacion por la izquierda
h = 0,10471976
n
f(x)
0,52359878 0,75983569
0,62831853 0,85237464
0,73303829 0,94889622
0,83775804 1,05385602
11
0,9424778 1,17319304
1,04719755 1,31607401
SUMATORIA
4,78815562
INTEGRAL
0,50141448
con n = 5
Aproximacion por la derecha
h = 0,10471976
n
f(x)
0,52359878 0,75983569
0,62831853 0,85237464
0,73303829 0,94889622
0,83775804 1,05385602
0,9424778 1,17319304
1,04719755 1,31607401
SUMATORIA
5,34439394
INTEGRAL
0,55966363
con n = 5
Aproximacion por el punto medio
h = 0,10471976
n
f(x)
0,52359878
0,8058583
0,62831853
0,89988001
0,73303829
0,83775804
1,11125927
0,9424778
1,24091296
1,04719755
SUMATORIA
5,05791053
INTEGRAL
0,52966315
12
2-
METODO DEL TRAPECIO
FUNCION =
f(x)= sin(x2)
Valor de pi =
3,141592654
Limite superior de la Integral (b) =
1,772453851
Limite inferior de la Integral (a) =
Valor de sub-intervalos (n) =
con n = 6
h = 0,295408975
n
X
0
f(x)
0
0,295408975 0,087155743
0,59081795 0,342020143
0,886226925 0,707106781
1,181635901 0,984807753
1,477044876 0,819152044
1,772453851 1,45478E-15
SUMATORIA
2,940242465
INTEGRAL
0,868574013
METODO DE SIMPSON
FUNCION =
f(x)= sin(x2)
Valor de pi =
3,141592654
13
Limite superior de la Integral (b) =
1,772453851
Limite inferior de la Integral (a) =
Valor de sub-intervalos (n) =
con n = 6
h = 0,295408975
n
f(x)
f(xi ) * K
0,295408975
0,087155743
0,348622971
0,59081795
0,342020143
0,684040287
0,886226925
0,707106781
2,828427125
1,181635901
0,984807753
1,969615506
1,477044876
0,819152044
3,276608177
1,772453851
5,66604E-16
5,66604E-16
SUMATORIA
9,107314066
INTEGRAL
0,896794105
INTEGRACION NUMERICA A PARTIR DE LA
SUMA DE RIEMANN
FUNCION =
f(x)= sin(x2)
Valor de pi =
3,141592654
Limite superior de la Integral (b) =
1,77245385
Limite inferior de la Integral (a) =
Valor de sub-intervalos (n) =
con n = 6
Aproximacion por la izquierda
h = 0,295408975
14
X
0
f(x)
0
0,295408975 0,087155743
0,59081795 0,342020143
0,886226925 0,707106781
1,181635901 0,984807753
1,477044876 0,819152044
1,772453851 1,45478E-15
SUMATORI
A
2,940242465
INTEGRAL
0,868574013
con n = 6
Aproximacion por la derecha
h = 0,295408975
n
X
0
f(x)
0
0,295408975 0,087155743
0,59081795 0,342020143
0,886226925 0,707106781
1,181635901 0,984807753
1,477044876 0,819152044
1,772453851 1,45478E-15
SUMATORI
A
2,940242465
INTEGRAL
0,868574013
con n = 6
Aproximacion por el punto medio
h = 0,295408975
n
X
0
f(x)
0
0,02181489
15
0,295408975
0,19509032
0,59081795
0,51877326
0,886226925
0,87672676
1,181635901
0,98078528
1,477044876
0,48098877
1,772453851
SUMATORI
A
2,593190501
INTEGRAL
0,766051748
3-
METODO DEL TRAPECIO
FUNCION =
f(x)= cos(x2)
Valor de pi =
3,141592654
Limite superior de la Integral (b) =
1,253314137
Limite inferior de la Integral (a) =
-1,25331414
Valor de sub-intervalos (n) =
con n = 4
h = 0,626657069
n
f(x)
-1,253314137
2,83302E-16
-0,626657069
0,923879533
0,626657069
0,923879533
1,253314137
2,83302E-16
16
SUMATORIA
2,847759065
INTEGRAL
1,784568348
METODO DE SIMPSON
FUNCION =
f(x)= cos(x2)
Valor de pi =
3,141592654
Limite superior de la Integral (b) =
1,2533141
37
1,2533141
4
Limite inferior de la Integral (a) =
Valor de sub-intervalos (n) =
con n = 4
h = 0,626657069
n
f(x)
-1,253314137
2,83302E-16
-0,626657069
0,923879533
0,626657069
0,923879533
1,253314137
2,83302E-16
SUMATORI
A
9,39103626
INTEGRAL
1,961653085
f(xi ) * K
2,83302E16
3,6955181
3
2
3,6955181
3
2,83302E16
17
INTEGRACION NUMERICA A PARTIR DE LA SUMA
DE RIEMANN
FUNCION =
f(x)= cos(x2)
Valor de pi =
3,141592654
Limite superior de la Integral (b) =
1,2533141
4
1,2533141
4
Limite inferior de la Integral (a) =
Valor de sub-intervalos (n) =
con n = 4
Aproximacion por la izquierda
h = 0,626657069
f(x)
-1,253314137
2,83302E-16
-0,626657069
0,923879533
0,626657069
0,923879533
1,253314137
2,83302E-16
SUMATORI
A
2,847759065
INTEGRAL
1,784568348
con n = 4
Aproximacion por la derecha
h = 0,626657069
n
f(x)
-1,253314137
2,83302E-16
-0,626657069
0,923879533
0,626657069
0,923879533
18
1,253314137
SUMATORI
A
2,847759065
INTEGRAL
1,784568348
2,83302E-16
con n = 4
Aproximacion por el punto medio
h = 0,626657069
n
X
0
-1,253314137
-0,626657069
0,626657069
1,253314137
SUMATORI
A
3,259156022
INTEGRAL
2,042373159
f(x)
0,6343932
8
0,9951847
3
0,9951847
3
0,6343932
8
4-
METODO DEL TRAPECIO
FUNCION =
f(x)= 1/ln(x)
Limite superior de la Integral (b) =
10 =
10
Limite inferior de la Integral (a) =
Valor de sub-intervalos (n) =
2
8
19
con n = 8
h= 1
n
f(x)
2 1,442695041
3 0,910239227
5 0,621334935
6 0,558110627
7 0,513898342
8 0,480898347
9 0,455119613
10 0,434294482
SUMATORIA
4,260948611
INTEGRAL
5,261351445
0,72134752
METODO DE SIMPSON
FUNCION =
f(x)= 1/ln(x)
Limite superior de la Integral (b) =
10 =
10
Limite inferior de la Integral (a) =
Valor de sub-intervalos (n) =
con n = 8
h= 1
n
f(x)
f(xi ) * K
20
2 1,442695041
1 1,44269504
3 0,910239227
4 3,64095691
0,72134752
2 1,44269504
5 0,621334935
4 2,48533974
6 0,558110627
2 1,11622125
7 0,513898342
4 2,05559337
8 0,480898347
2 0,96179669
9 0,455119613
4 1,82047845
10 0,434294482
1 0,43429448
SUMATORIA
15,40007098
INTEGRAL
5,133356993
INTEGRACION NUMERICA A PARTIR DE LA
SUMA DE RIEMANN
FUNCION =
f(x)= 1/ln(x)
Limite superior de la Integral (b) =
10
10
Limite inferior de la Integral (a) =
Valor de sub-intervalos (n) =
con n = 8
Aproximacion por la izquierda
h= 1
f(x)
2 1,44269504
3 0,91023923
4 0,72134752
5 0,62133493
6 0,55811063
7 0,51389834
8 0,48089835
21
9 0,45511961
10 0,43429448
SUMATORIA
5,703643652
INTEGRAL
5,703643652
con n = 8
Aproximacion por la derecha
h= 1
n
f(x)
2 1,44269504
3 0,91023923
4 0,72134752
5 0,62133493
6 0,55811063
7 0,51389834
8 0,48089835
9 0,45511961
10 0,43429448
SUMATORIA
4,695243093
INTEGRAL
4,695243093
con n = 8
Aproximacion por el punto medio
h= 1
n
f(x)
1,09135667
0,7982356
0,6648594
0,58659693
0,53424449
0,4963018
0,46727527
0,44418942
22
10
SUMATORIA
5,083059588
INTEGRAL
5,083059588
5-
METODO DEL TRAPECIO
FUNCION =
f(x)= arctan(1/(1+x2))
Valor de pi =
3,141592654
Limite superior de la Integral (b) =
3,1415927
Limite inferior de la Integral (a) =
0,7853982
Valor de sub-intervalos (n) =
con n = 6
h = 0,392699082
n
f(x)
0,785398163 0,553901703
1,178097245 0,396586824
1,570796327 0,280781318
1,963495408 0,203119604
2,35619449 0,151464252
2,748893572 0,116343571
3,141592654 0,091741419
SUMATORIA
1,148295568
INTEGRAL
0,577706346
23
METODO DE SIMPSON
FUNCION =
f(x)= arctan(1/(1+x2))
Valor de pi =
3,141592654
Limite superior de la Integral (b) =
3,14159265
Limite inferior de la Integral (a) =
0,78539816
Valor de sub-intervalos (n) =
con n = 6
h = 0,392699082
n
f(x)
f(xi ) * K
0,785398163 0,553901703
1,178097245 0,396586824
4 1,58634729
1,570796327 0,280781318
2 0,56156264
1,963495408 0,203119604
4 0,81247841
2,35619449 0,151464252
2,748893572 0,116343571
4 0,46537429
3,141592654 0,091741419
2 0,18348284
SUMATORIA
4,466075674
INTEGRAL
0,584607939
0,5539017
0,3029285
INTEGRACION NUMERICA A PARTIR DE
LA SUMA DE RIEMANN
FUNCION =
f(x)= arctan(1/(1+x2))
Valor de pi =
3,141592654
Limite superior de la Integral (b) =
3,14159265
24
Limite inferior de la Integral (a) =
0,78539816
Valor de sub-intervalos (n) =
con n = 6
Aproximacion por la izquierda
h = 0,392699082
f(x)
0,785398163
0,5539017
1,178097245 0,39658682
1,570796327 0,28078132
1,963495408
0,2031196
2,35619449 0,15146425
2,748893572 0,11634357
3,141592654 0,09174142
SUMATORIA
1,702197271
INTEGRAL
0,668451305
con n = 6
Aproximacion por la derecha
h = 0,392699082
n
f(x)
0,785398163
0,5539017
1,178097245 0,39658682
1,570796327 0,28078132
1,963495408
0,2031196
2,35619449 0,15146425
2,748893572 0,11634357
3,141592654 0,09174142
SUMATORIA
1,240036987
INTEGRAL
0,486961386
25
con n = 6
Aproximacion por el punto medio
h = 0,392699082
n
f(x)
0,785398163
0,47098797
1,178097245
0,33322129
1,570796327
0,23795785
1,963495408
0,17472466
2,35619449
0,13228167
2,748893572
0,10299919
3,141592654
SUMATORIA
1,349173443
INTEGRAL
0,529819172
6-
METODO DEL TRAPECIO
f(x)= e-x^2
FUNCION =
Limite superior de la Integral (b) =
10
Limite inferior de la Integral (a) =
-10
Valor de sub-intervalos (n) =
20
con n = 20
h= 1
n
f(x)
1 0,367879441
26
2 0,018315639
4 1,12535E-07
5 1,38879E-11
6 2,31952E-16
7 5,24289E-22
8 1,60381E-28
9 6,63968E-36
10
10 3,72008E-44
11
11 2,82077E-53
12
12 2,89464E-63
13
13 4,02006E-74
14
14 7,55582E-86
15
15 1,92195E-98
16
16 6,6163E-112
17
17 3,0824E-126
18
18 1,9435E-141
19
19 1,6584E-157
20
20 1,9152E-174
SUMATORIA
0,386318602
INTEGRAL
1,772637205
0,00012341
METODO DE SIMPSON
FUNCION =
f(x)= e-x^2
Limite superior de la Integral (b) =
Limite inferior de la Integral (a) =
Valor de sub-intervalos (n) =
20
con n = 20
h= 1
27
f(x)
1 0,36787944
4 1,47151776
2 0,01831564
2 0,03663128
3 0,00012341
4 0,00049364
4 1,1254E-07
2 2,2507E-07
5 1,3888E-11
4 5,5552E-11
6 2,3195E-16
4,639E-16
7 5,2429E-22
2,8
1,468E-21
8 1,6038E-28
2,8 4,4907E-28
9 6,6397E-36
2,8 1,8591E-35
10
10 3,7201E-44
2,8 1,0416E-43
11
11 2,8208E-53
2,8 7,8982E-53
12
12 2,8946E-63
2,8
13
13 4,0201E-74
2,8 1,1256E-73
14
14 7,5558E-86
2,8 2,1156E-85
15
15 1,9219E-98
2,8 5,3815E-98
16
16 6,616E-112
2,8 1,853E-111
17
17 3,082E-126
2,8 8,631E-126
18
18 1,944E-141
2,8 5,442E-141
19
19 1,658E-157
2,8 4,644E-157
20
20 1,915E-174
2,8 5,362E-174
SUMATORIA
2,508642907
INTEGRAL
1,672428605
f(xi ) * K
1
8,105E-63
INTEGRACION NUMERICA A PARTIR DE
LA SUMA DE RIEMANN
FUNCION =
f(x)= e-x^2
Limite superior de la Integral (b) =
10
Limite inferior de la Integral (a) =
28
Valor de sub-intervalos (n) =
20
con n = 20
Aproximacion por la izquierda
h= 1
n
f(x)
1 0,36787944
2 0,01831564
3 0,00012341
4 1,1254E-07
5 1,3888E-11
6 2,3195E-16
7 5,2429E-22
8 1,6038E-28
9 6,6397E-36
10
10 3,7201E-44
11
11 2,8208E-53
12
12 2,8946E-63
13
13 4,0201E-74
14
14 7,5558E-86
15
15 1,9219E-98
16
16 6,616E-112
17
17 3,082E-126
18
18 1,944E-141
19
19 1,658E-157
20
20 1,915E-174
SUMATORIA
1,386318602
INTEGRAL
1,386318602
con n = 20
Aproximacion por la derecha
29
h= 1
n
f(x)
1 0,36787944
2 0,01831564
3 0,00012341
4 1,1254E-07
5 1,3888E-11
6 2,3195E-16
7 5,2429E-22
8 1,6038E-28
9 6,6397E-36
10
10 3,7201E-44
11
11 2,8208E-53
12
12 2,8946E-63
13
13 4,0201E-74
14
14 7,5558E-86
15
15 1,9219E-98
16
16 6,616E-112
17
17 3,082E-126
18
18 1,944E-141
19
19 1,658E-157
20
20 1,915E-174
SUMATORIA
0,386318602
INTEGRAL
0,386318602
con n = 20
Aproximacion por el punto medio
h= 1
n
f(x)
0,77880078
0,10539922
30
0,00193045
4,7851E-06
1,6052E-09
7,2877E-14
4,4777E-19
3,7234E-25
4,1901E-32
6,3815E-40
10
10
1,3153E-48
11
11
3,6691E-58
12
12
1,3851E-68
13
13
7,0767E-80
14
14
4,8931E-92
15
15
4,579E-105
16
16
5,799E-119
17
17
9,938E-134
18
18
2,305E-149
19
19
7,236E-166
20
20
SUMATORIA
0,886135248
INTEGRAL
0,886135248
7-
METODO DEL TRAPECIO
FUNCION =
f(x)= cos(ln(x))
Limite superior de la Integral (b) =
Limite inferior de la Integral (a) =
Valor de sub-intervalos (n) =
1
8
31
con n = 8
h = 0,125
n
f(x)
1,125 0,99307159
1,25 0,97520661
1,375 0,94972069
1,5 0,91891904
1,625 0,88443809
1,75 0,84745914
1,875 0,80884643
SUMATORIA
6,377661581
INTEGRAL
0,907785129
0,7692389
METODO DE SIMPSON
FUNCION =
f(x)= cos(ln(x))
Limite superior de la Integral (b) =
Limite inferior de la Integral (a) =
Valor de sub-intervalos (n) =
con n = 8
h = 0,125
n
X
0
1
f(x)
1
k
1
1,125 0,99307159
f(xi ) * K
1
4 3,97228637
32
1,25 0,97520661
2 1,95041323
1,375 0,94972069
4 3,79888275
1,5 0,91891904
2 1,83783807
1,625 0,88443809
4 3,53775235
1,75 0,84745914
2 1,69491827
1,875 0,80884643
4 3,23538571
SUMATORIA
21,79671565
INTEGRAL
0,908196486
0,7692389
0,7692389
INTEGRACION NUMERICA A PARTIR DE LA
SUMA DE RIEMANN
FUNCION =
f(x)= cos(ln(x))
Limite superior de la Integral (b) =
Limite inferior de la Integral (a) =
Valor de sub-intervalos (n) =
con n = 8
Aproximacion por la izquierda
h = 0,125
X
0
f(x)
1
1,125 0,993071594
1,25 0,975206613
1,375 0,949720688
1,5 0,918919036
1,625 0,884438087
1,75 0,847459137
1,875 0,808846427
2 0,769238901
33
SUMATORIA
7,377661581
INTEGRAL
0,922207698
con n = 8
Aproximacion por la derecha
h = 0,125
n
X
0
f(x)
1
1,125 0,993071594
1,25 0,975206613
1,375 0,949720688
1,5 0,918919036
1,625 0,884438087
1,75 0,847459137
1,875 0,808846427
2 0,769238901
SUMATORIA
INTEGRAL
7,146900482
0,89336256
con n = 8
Aproximacion por el punto medio
h = 0,125
n
f(x)
0,99816289
1,125
0,98527005
1,25
0,96325331
1,375
0,93486934
1,5
0,90205587
1,625
0,86620067
1,75
0,82831263
1,875
0,78913403
34
SUMATORIA
7,267258795
INTEGRAL
0,908407349
8-
METODO DEL TRAPECIO
f(x)= arcsec(x2)
FUNCION =
Limite superior de la Integral (b) =
10
10
Limite inferior de la Integral (a) =
Valor de sub-intervalos (n) =
con n = 9
h= 1
n
f(x)
1,3181
1,4595
1,5083
1,5308
1,5430
1,5504
1,5552
1,5585
10
1,5608
35
SUMATORIA
12,0238
INTEGRAL
12,8042
METODO DE SIMPSON
FUNCION
=
f(x)= arcsec(x2)
Limite superior de la Integral (b) =
10
10
Limite inferior de la Integral (a) =
Valor de sub-intervalos (n) =
El metodo de aproximacion de Simpsom no acepta valores de subintervalos(n) impares
INTEGRACION NUMERICA A PARTIR DE LA
SUMA DE RIEMANN
f(x)= arcsec(x2)
FUNCION =
Limite superior de la Integral (b) =
10
10
Limite inferior de la Integral (a) =
Valor de sub-intervalos (n) =
con n = 9
Aproximacion por la izquierda
h= 1
n
f(x)
1,3181
1,4595
1,5083
36
1,5308
1,5430
1,5504
1,5552
1,5585
10
1,5608
SUMATORIA
12,0238
INTEGRAL
12,0238
con n = 9
Aproximacion por la derecha
h= 1
n
f(x)
1,3181
1,4595
1,5083
1,5308
1,5430
1,5504
1,5552
1,5585
10
1,5608
SUMATORIA
13,5846
INTEGRAL
13,5846
con n = 9
Aproximacion por el punto medio
h= 1
n
f(x)
1,1102
1,4101
1,4891
1,5214
1,5377
37
1,5471
1,5530
1,5570
1,5597
10
SUMATORIA
11,7256
INTEGRAL
11,7256
9-
METODO DEL TRAPECIO
FUNCION =
f(x)= sin(cos(tan(ln(x))))
Limite superior de la Integral (b) =
5/2
2,50
Limite inferior de la Integral (a) =
0,50
Valor de sub-intervalos (n) =
con n = 4
h = 0,5
n
f(x)
0,5 0,62443089
1 0,84147098
1,5 0,78906044
2 0,62443089
2,5 0,26133803
38
SUMATORIA
2,254962315
INTEGRAL
1,348923388
METODO DE SIMPSON
FUNCION =
f(x)= sin(cos(tan(ln(x))))
Limite superior de la Integral (b) =
5/2
2,50
Limite inferior de la Integral (a) =
0,50
Valor de sub-intervalos (n) =
con n = 4
h = 0,5
n
f(x)
f(xi ) * K
0,5 0,62443089
1 0,62443089
1 0,84147098
4 3,36588394
1,5 0,78906044
2 1,57812089
2 0,62443089
4 2,49772355
2,5 0,26133803
1 0,26133803
SUMATORIA
8,327497295
INTEGRAL
1,387916216
INTEGRACION NUMERICA A PARTIR DE
LA SUMA DE RIEMANN
39
FUNCION =
f(x)= sin(cos(tan(ln(x))))
Limite superior de la Integral (b) =
5/2
2,50
Limite inferior de la Integral (a) =
1/2
0,50
Valor de sub-intervalos (n) =
con n = 4
Aproximacion por la izquierda
h = 0,5
f(x)
0,5 0,62443089
1 0,84147098
1,5 0,78906044
2 0,62443089
2,5 0,26133803
SUMATORIA
2,879393202
INTEGRAL
1,439696601
con n = 4
Aproximacion por la derecha
h = 0,5
n
f(x)
0,5 0,62443089
1 0,84147098
1,5 0,78906044
2 0,62443089
2,5 0,26133803
SUMATORIA
2,516300349
INTEGRAL
1,258150174
40
con n = 4
Aproximacion por el punto medio
h = 0,5
n
SUMATORIA
f(x)
0,5
0,81720399
0,82734447
1,5
0,72438006
0,47545036
2,5
2,844378881
INTEGRAL
1,42218944
10-
METODO DEL TRAPECIO
FUNCION =
f(x)= (COS(|x|))1/2
Valor de pi =
3,141592654
Limite superior de la Integral (b) =
1,57079633
1,57079633
Limite inferior de la Integral (a) =
Valor de sub-intervalos (n) =
con n = 8
h = 0,392699082
41
f(x)
-1,570796327
7,82671E-09
-1,178097245
0,618614122
-0,785398163
0,840896415
-0,392699082
0,961186523
0,392699082
0,961186523
0,785398163
0,840896415
1,178097245
0,618614122
1,570796327
7,82671E-09
SUMATORIA
5,841394122
INTEGRAL
2,49025965
METODO DE
SIMPSON
FUNCION =
f(x)= (COS(|x|))1/2
Valor de pi =
3,141592654
1,5707963
3
1,5707963
3
Limite superior de la Integral (b) =
Limite inferior de la Integral (a) =
Valor de sub-intervalos (n) =
con n = 8
h = 0,392699082
n
f(x)
-1,570796327
7,82671E-09
-1,178097245
0,618614122
f(xi ) * K
1 7,8267E-09
2,4744564
4
9
42
1,6817928
3
3,8447460
9
-0,785398163
0,840896415
-0,392699082
0,961186523
0,392699082
0,961186523
0,785398163
0,840896415
1,178097245
0,618614122
1,570796327
7,82671E-09
1 7,8267E-09
SUMATORI
A
18,00199084
INTEGRAL
2,356455091
2
3,8447460
9
1,6817928
3
2,4744564
9
INTEGRACION NUMERICA A PARTIR DE LA
SUMA DE RIEMANN
FUNCION =
f(x)= (COS(|x|))1/2
Valor de pi =
3,141592654
Limite superior de la Integral (b) =
1,570796327
Limite inferior de la Integral (a) =
-1,570796327
Valor de sub-intervalos (n) =
con n = 8
Aproximacion por la izquierda
h = 0,392699082
f(x)
-1,570796327
7,82671E-09
-1,178097245
0,618614122
-0,785398163
0,840896415
-0,392699082
0,961186523
1
43
0,392699082
0,961186523
0,785398163
0,840896415
1,178097245
0,618614122
8
SUMATORI
A
1,570796327
7,82671E-09
INTEGRAL
2,293910111
5,84139413
con n = 8
Aproximacion por la derecha
h = 0,392699082
n
f(x)
-1,570796327
7,82671E-09
-1,178097245
0,618614122
-0,785398163
0,840896415
-0,392699082
0,961186523
0,392699082
0,961186523
0,785398163
0,840896415
1,178097245
0,618614122
8
SUMATORI
A
1,570796327
7,82671E-09
INTEGRAL
2,293910111
5,84139413
con n = 8
Aproximacion por el punto medio
h = 0,392699082
n
f(x)
-1,570796327
0,441690301
-1,178097245
0,745365838
-0,785398163
0,911849556
-0,392699082
0,990346041
0,990346041
0,392699082
0,911849556
44
0,785398163
0,745365838
1,178097245
0,441690301
1,570796327
SUMATORI
A
6,178503472
INTEGRAL
2,42629264
45
PARTE III
1. y = 4-(x-3)^2;
y = abs(x-2)
a)
Densidad constante
A=
P
=> A=P+Q = 1.2695 + 5.4893 =6.7588
SECCIN P
4-(x-3)^2=-x^2+6x-5
Coordenada X = 1/P
=( 1/1.2695) 2.1933 = 1.7277
Coordenada Y = (1/2P)
=(1/2(1.2695))(3.3234)=
1.3089
Centroide P: (1.7277, 1.3089)
SECCIN Q
Coordenada X = 1/Q
= (1/5.4893)15.9733 = 2.9099
Coordenada Y = (1/2Q)
= (1/2(5.4893))(25.1616)=
2.2919
46
Centroide Q: (2.9099, 2.2919)
EN SUPERFICIE TOTAL
X= [(1.2695(1.3089)+ 5.4893(2.2919)]/A= 2.1073
Y= [(1.2695(1.7277)+ 5.4893(2.9099)]/A= 2.6878
Centro de gravedad del sistema: (2.1073, 2.6878)
Densidad variable
(x^2+sin(x))
SECCIN P
m1 =
= 5.066
M(y)=
= 8.899
Coordenana X = M(y)/m= 1.757
M(x)= 1/2
= (1/2)13.503 = 6.752
Coordenana Y = M(x)/m = 1.333
Centroide P: (1.757, 1.333)
SECCIN Q
m2=
= 49.340
M(y)=
= 152.475
Coordenana X = M(y)/m = 3.090
M(x)= 1/2
=234.662
Coordenana Y = M(x)/m = 4.756
Centroide Q: (3.090,4.756)
47
EN SUPERFICIE TOTAL
X = [(m1(Y1)+m2(Y2)]/(m1+m2)]=[ (5.066(1.333)+ 49.340(4.756)]/(54.406)= 4.441
Y = [(m1(X1)+m2(X2)]/(m1+m2)]=[( (5.066(1.757)+ 49.340(3.090)]/(54.406)= 2.985
Centro de gravedad = (4.441, 2.985)
2. y = sqrt(abs(x-1)); y = -1/10ln(cos(x))
2 a)
Densidad constante
48
3. y = -(4-(x-5)^2); y = 4-(x-5)^2; x = y^2; x = 8-y^2
MASAS
Cada una se multiplica por 2 pues slo se
utiliza el rea sobre el eje X , ya que sta
rea es simtrica.
2P
m1 =
P
P
=11.969
S
S
2Q
m2 =
=
= 25.448
2R
m3 =
= 241.935
2S
m4 =
= 16.832
m=296.186
6.69
En los M(y) se multiplica la densidad por x por 2rea. Cmo = x^2. Se
simplifica para obtener x^3. As: x(x^2)(2rea)= x^3(2rea)
Adems se multiplica por 2 por la simetra antes mencionada.
En los M(x) tenemos como integrando: /2[f(x) +g(x)]2(f(x) -g(x))dx entonces al
simplificarse, sabemos que nuestro g(x) = 0 por lo tanto el integrando ser de la forma:
f(x) dx.
SECCIN P
.
X= 1/m1
= 40.278/11.969 = 3.365
49
Y=
= 7.851
Coordenadas P: (3.365, 7.851)
SECCIN Q
X= 1/m2
= 96.235/25.448 = 3.782
Y=
=24.740
Coordenadas Q: (3.782, 24.740)
SECCIN R
X= 1/m3
= 1318.34/241.935 = 5.449
Y=
= 191.005
Coordenadas R: (5.449, 191.005)
SECCIN S
X= 1/m4
= 114.39/16.832 = 6.796
Y=
= 6.461
Coordenadas R: (6.796, 6.461)
Centro de gravedad:
X = [(m1(Y1)+m2(Y2)+m3(Y3) +m4(Y4)]/(m1+m2+m3+m4)
= 47043.098/296.186
=158.83
X = [(m1(X1)+m2(X2)+m3(X3) +m4(X4)]/(m1+m2+m3+m4)
= 1567.980/296.186
= 5.29
50
CONCLUSIONES
Al desarrollo de todas estas integrales nos podemos dar cuenta de que es posible afirmar
que cualquier integral propia o impropia puede resolverse con estos mtodos, con la nica
condicin de que el resultado deseado sea el rea bajo la curva.
Tambin es posible decir con certeza qu; esto abre expectativas muy alentadoras para
poder evitar el enorme obstculo que pueda llegar a representar una integral demasiado
complicada al momento de querer resolverla de una forma analtica.
RECOMENDACIONES
Es importante recalcar que a la hora de ejecutar alguno de estos mtodos se tome en cuenta
que estos tienen un margen de error que est relacionado con el nmero de particiones.
Los datos generados con estos mtodos de integrales numricas tienden a ser fciles si se
toman en cuenta los pasos correctamente.
El resultado final de una integral evaluada con alguno de estos mtodos puede que este muy
contrastada con la que pueda obtenerse de una integral de modo analtico o generada en
calculadora y/o computadora.
REFERENCIAS.
Leithold, L. (1973). Integracion numerica. En L. Leithold, Calculo con geometra analitica segunda
edicion (pgs. 544-553). California: Harper and Row lationoamerica.
51