Curso de Métodos Numéricos

Tarea 3- Diferenciación e Integración y Ecuaciones Diferenciales en Valor Inicial

Presentado a:
Jorge Armando Amador

Entregado por:

Sandra Fanny Liliana Picón Díaz

1013692047
María Elena Montoya
49785077
César Augusto Pava Linares
79119166

Grupo: 100401_28

Universidad Nacional Abierta y a Distancia - UNAD

Escuela de Ciencias Básicas, Ingenierías y Tecnologías
Escuela de Ciencias de la Educación
Bogotá D.C. 28 de noviembre 2019
 Tabla de contenido

Introducción ..................................................................................................................................... 3
Objetivos.......................................................................................................................................... 4
General ......................................................................................................................................... 4
Específicos .................................................................................. Error! Bookmark not defined.
Desafio 2 ......................................................................................... Error! Bookmark not defined.
Desafío 3 ........................................................................................................................................ 28
Desafío 4 ........................................................................................................................................ 35
Conclusiones.................................................................................................................................. 42
Referencias Bibliográficas ............................................................................................................. 43

2
Introducción

3
Objetivos

Objetivo General
Resolver problemas aplicando diferenciación e integración numéricas y ecuaciones
diferenciales en valor inicial.

Objetivos Específicos
 Reconocer la diferenciación numérica con base en el desafío planteado, así como la
Regla del Trapecio, Regla de Simpson.
 Aplicar la integración de Romberg e integrales múltiples, según el problema
planteado.
 Diferenciar los Métodos de Euler, Métodos de Taylor y Método de Runge Kutta,
que nos permita su comprensión en el desarrollo de los problemas propuestos.
.

4
Tabla 1: Rol de trabajo y desafío.
Nombre Estudiante Rol Desafío
María Elena Montoya Alertas 2
Evaluador
Cesar Augusto Pava Linares Compilador 4
Revisor
Sandra Fanny Liliana Picón Díaz Entregas 3
Fuente: Tutoría.

Aporte 1. Diferenciación numérica

Individuales
Desafío 2
Evalúe la ganancia en la precisión usando las fórmulas de diferencias finitas hacia adelante,
atrás y central, en la evaluación de la primera derivada de la función 𝑓(𝑥) = 𝑒 𝑥 en 𝑥 = 1,
con ℎ = 0.1, comparada con la solución analítica (exacta).

Desafío 3: Diferencias finitas

Dada la función 𝑓(𝑥) = 𝑡𝑎𝑛(40𝑥), encuentre 𝑓 ′ (0.175) usando las representaciones de

diferencias finitas hacia adelante, atrás y central, con ℎ = 0.075. Compare los resultados obtenidos
con la solución analítica. ¿Tienen sentido las respuestas obtenidas? Haga un análisis basado en la
teoría y el conocimiento de la función acerca de lo que sucede en este caso.

ECUACIÓN RAZÓN O EXPLICACIÓN

5
Desafío 4
Dada la función 𝑓(𝑥) = 𝑠𝑒𝑛(𝑥), encuentre 𝑓 ′ (1) usando los esquemas de diferencias
finitas hacia adelante, atrás y central, con ℎ = 0.2. ¿Es este un tamaño de paso
suficientemente pequeño para este problema? Haga un análisis justificando su respuesta
con la teoría revisada.

6
Aporte 2. Integración numérica

Individuales
Cálculo de la integral correspondiente a su número de desafío con cada uno de los métodos
indicados
Desafío 2:
𝟏
𝟒
∫ 𝑑𝑥
𝟎 𝟏 + 𝒙𝟐

Desafío 3:
𝟏
𝟏
∫ 𝟐
𝑑𝑥
−𝟏 𝟏 + 𝟏𝟎𝟎𝒙

Desafío 3:
a. Regla del trapecio simple y compuesta

Cálculo de la integral:
𝟏
𝟏
∫ 𝟐
𝑑𝑥
−𝟏 𝟏 + 𝟏𝟎𝟎𝒙

Para las reglas compuestas use 𝑛 = 6 en todos los casos

ECUACIÓN RAZÓN O
EXPLICACIÓN

7
Desafío 3:
b. Regla de Simpson1/3 simple y compuesta

Cálculo de la integral:
𝟏
𝟏
∫ 𝟐
𝑑𝑥
−𝟏 𝟏 + 𝟏𝟎𝟎𝒙

Para las reglas compuestas use 𝑛 = 6 en todos los casos

ECUACIÓN RAZÓN O EXPLICACIÓN

Desafío 3:
c. Regla de Simpson3/8 simple y compuesta

Cálculo de la integral:
𝟏
𝟏
∫ 𝟐
𝑑𝑥
−𝟏 𝟏 + 𝟏𝟎𝟎𝒙

Para las reglas compuestas use 𝑛 = 6 en todos los casos

ECUACIÓN RAZÓN O EXPLICACIÓN

Desafío 4:
𝟏
∫ √𝒙𝒍𝒐𝒈(𝒙)𝑑𝑥
𝟎

a. REGLA DEL TRAPECIO

8
Inicialmente se resuelve la integral como se hace en cálculo para luego comparar el
resultado exacto con respecto al resultado numérico y poder calcular el error.


0
x log xdx

1 1
ln x 1
Como log x 
ln 10
entonces 
0
x log xdx 
ln 10 0
x ln xdx

Ahora bien, la integral  x ln xdx se resuelve por partes, recurriendo a la fórmula

3
dx 2
 udv  uv   vdu ; es decir si u  ln x y v  xdx , entonces du  x
y v  x 2 ; por lo
3
3 3 3
2 x 2 ln x 4 x 2 2 x 2  2
3 1
2 ln x 2 2 2
tanto,  x ln xdx 
3
x   x dx 
3 3

9

3   ln x  
3

3
2x  2
2
Finalmente,  x log xdx  
3 ln 10 
ln x    C
3

El problema de esta integral es que la función f ( x)  x log xdx no está definida en x=0;
el valor de la integral se podría calcular haciendo que el límite inferior sea un número
cercano a cero, pero mayor que cero. Es decir, que el valor exacto de la integral no se puede
calcular aun cuando se pueda encontrar la solución de la integral.

De acuerdo a lo anterior utilizando Excel hice un cálculo de la integral haciendo

variaciones del límite inferior acercándome a cero con potencias negativas de 10.

La tabla que a continuación expongo resume estos resultados:

2/(3*LN(10)
) x x^(3/2) LN(x) LN(X)-(2/3)
- -
0,66666666 0,1930197
0,28952965 1 1 0 7 7 -0,1658341
0,28952965 0,1 0,0316227 - -2,96925176 - -
9
 8 2,3025850 0,0271856 0,1914934
9 7 2
- - - -
4,6051701 5,27183685 0,0015263 0,1929504
0,28952965 0,01 0,001 9 3 5 2
- - -
3,1623E- 6,9077552 7,57442194 -6,9349E- 0,1930169
0,28952965 0,001 05 8 6 05 1
- - -
9,2103403 9,87700703 -2,8597E- 0,1930196
0,28952965 0,0001 0,000001 7 9 06 6
- - -
3,1623E- 11,512925 12,1795921 -1,1151E- 0,1930197
0,28952965 0,00001 08 5 3 07 7
- - -
13,815510 14,4821772 0,1930197
0,28952965 0,000001 1E-09 6 2 -4,193E-09 7
- -
3,1623E- 16,118095 16,7847623 -1,5368E-
0,28952965 0,0000001 11 7 2 10
VALORES DE LA INTEGRAL HACIENDO VARIACIONES DEL LIMITE INFERIOR

Utilizando la regla del trapecio se realizaron 4 particiones:

x sqrx Lnx
- -
0,25 0,5 1,38629436 1 0,69314718
- -
0,5 0,70710678 0,69314718 2 0,98025814
- -
0,75 0,8660254 0,28768207 2 0,49827997
1 1 0 1 0
-
2,17168529

0,25
0,125
2,30258509
VALOR DE LA INTEGRAL CON 4 -
TRAPECIOS 0,11789387
Si asumimos que el valor de la integral con el límite inferior igual a 0,0000001 es el 100%
entonces el porcentaje de error con el valor -0,11789387 es 38,92%

10
 x sqrx lnx ln10
0 0
- -
0,1 0,31622777 2,30258509 2,30258509 1 0,316227766
- -
0,2 0,4472136 1,60943791 2,30258509 2 0,625177778
-
0,3 0,54772256 -1,2039728 2,30258509 2 0,572784967
-
0,4 0,63245553 0,91629073 2,30258509 2 -0,50335872
- -
0,5 0,70710678 0,69314718 2,30258509 2 0,425720703
- -
0,6 0,77459667 0,51082562 2,30258509 2 0,343686605
- -
0,7 0,83666003 0,35667494 2,30258509 2 0,259200556
- -
0,8 0,89442719 0,22314355 2,30258509 2 0,173357901
- -
0,9 0,9486833 0,10536052 2,30258509 2 0,086818734
1 1 0 2,30258509 1 0
-3,30633373
-
VALOR INTEGRAL CON DIEZ TRAPECIOS 0,165316686

Si asumimos que el valor de la integral con el límite inferior igual a 0,0000001 es el 100%
entonces el porcentaje de error con el valor -0,165316686 es 14,35%

-
0,01 0,1 4,60517019 2,30258509 1 -0,2
- -
0,02 0,14142136 3,91202301 2,30258509 2 0,480541284
-
0,03 0,17320508 -3,5065579 2,30258509 2 0,527540672
- -
0,04 0,2 3,21887582 2,30258509 2 0,559176003
- -
0,05 0,2236068 2,99573227 2,30258509 2 0,581838302
- -
0,06 0,24494897 2,81341072 2,30258509 2 0,598581196
- -
0,07 0,26457513 2,65926004 2,30258509 2 0,611116675
- -
0,08 0,28284271 2,52572864 2,30258509 2 0,620506007

11
 - -
0,09 0,3 2,40794561 2,30258509 2 0,627454494
- -
0,1 0,31622777 2,30258509 2,30258509 2 0,632455532
- -
0,11 0,33166248 2,20727491 2,30258509 2 0,635868157
- -
0,12 0,34641016 2,12026354 2,30258509 2 0,637961947
- -
0,13 0,36055513 2,04022083 2,30258509 2 0,638944535
- -
0,14 0,37416574 1,96611286 2,30258509 2 0,638979269
- -
0,15 0,38729833 1,89711998 2,30258509 2 0,638196966
- -
0,16 0,4 1,83258146 2,30258509 2 0,636704014
- -
0,17 0,41231056 1,77195684 2,30258509 2 0,634588076
- -
0,18 0,42426407 1,71479843 2,30258509 2 0,631922234
- -
0,19 0,43588989 1,66073121 2,30258509 2 0,628768033
- -
0,2 0,4472136 1,60943791 2,30258509 2 0,625177778
- -
0,21 0,45825757 1,56064775 2,30258509 2 0,621196277
- -
0,22 0,46904158 1,51412773 2,30258509 2 0,616862204
- -
0,23 0,47958315 1,46967597 2,30258509 2 0,612209153
- -
0,24 0,48989795 1,42711636 2,30258509 2 0,607266482
- -
0,25 0,5 1,38629436 2,30258509 2 0,602059991
- -
0,26 0,50990195 1,34707365 2,30258509 2 0,596612463
- -
0,27 0,51961524 1,30933332 2,30258509 2 0,590944111
- -
0,28 0,52915026 1,27296568 2,30258509 2 0,585072945
- -
0,29 0,53851648 1,23787436 2,30258509 2 0,579015076
-
0,3 0,54772256 -1,2039728 2,30258509 2 0,572784967
- -
0,31 0,55677644 1,17118298 2,30258509 2 0,566395647

12
 -
0,32 0,56568542 1,13943428 2,30258509 2 -0,55985889
- -
0,33 0,57445626 1,10866262 2,30258509 2 0,553185367
- -
0,34 0,58309519 1,07880966 2,30258509 2 0,546384779
- -
0,35 0,59160798 1,04982212 2,30258509 2 0,539465965
- -
0,36 0,6 1,02165125 2,30258509 2 0,532436999
- -
0,37 0,60827625 0,99425227 2,30258509 2 0,525305275
- -
0,38 0,6164414 0,96758403 2,30258509 2 0,518077576
- -
0,39 0,6244998 0,94160854 2,30258509 2 0,510760142
-
0,4 0,63245553 0,91629073 2,30258509 2 -0,50335872
- -
0,41 0,64031242 0,89159812 2,30258509 2 0,495878614
- -
0,42 0,64807407 0,86750057 2,30258509 2 0,488324731
- -
0,43 0,65574385 0,84397007 2,30258509 2 0,480701614
- -
0,44 0,66332496 0,82098055 2,30258509 2 0,473013477
-
0,45 0,67082039 -0,7985077 2,30258509 2 0,465264236
- -
0,46 0,678233 0,77652879 2,30258509 2 0,457457534
- -
0,47 0,68556546 0,75502258 2,30258509 2 0,449596766
- -
0,48 0,69282032 0,73396918 2,30258509 2 0,441685098
- -
0,49 0,7 0,71334989 2,30258509 2 0,433725488
- -
0,5 0,70710678 0,69314718 2,30258509 2 0,425720703
- -
0,51 0,71414284 0,67334455 2,30258509 2 0,417673332
- -
0,52 0,72111026 0,65392647 2,30258509 2 0,409585803
- -
0,53 0,72801099 0,63487827 2,30258509 2 0,401460394
- -
0,54 0,73484692 0,61618614 2,30258509 2 0,393299244

13
 -
0,55 0,74161985 -0,597837 2,30258509 2 0,385104366
-
0,56 0,74833148 -0,5798185 2,30258509 2 0,376877652
- -
0,57 0,75498344 0,56211892 2,30258509 2 0,368620884
- -
0,58 0,76157731 0,54472718 2,30258509 2 0,360335745
- -
0,59 0,76811457 0,52763274 2,30258509 2 0,352023819
- -
0,6 0,77459667 0,51082562 2,30258509 2 0,343686605
- -
0,61 0,78102497 0,49429632 2,30258509 2 0,335325517
-
0,62 0,78740079 -0,4780358 2,30258509 2 0,326941894
- -
0,63 0,79372539 0,46203546 2,30258509 2 0,318537003
-
0,64 0,8 -0,4462871 2,30258509 2 0,310112042
- -
0,65 0,80622577 0,43078292 2,30258509 2 0,301668148
-
0,66 0,81240384 0,41551544 2,30258509 2 -0,2932064
- -
0,67 0,81853528 0,40047757 2,30258509 2 0,284727819
- -
0,68 0,82462113 0,38566248 2,30258509 2 0,276233378
- -
0,69 0,83066239 0,37106368 2,30258509 2 0,267723998
- -
0,7 0,83666003 0,35667494 2,30258509 2 0,259200556
- -
0,71 0,84261498 0,34249031 2,30258509 2 0,250663886
- -
0,72 0,84852814 0,32850407 2,30258509 2 0,242114782
- -
0,73 0,85440037 0,31471074 2,30258509 2 0,233553999
- -
0,74 0,86023253 0,30110509 2,30258509 2 0,224982256
-
0,75 0,8660254 0,28768207 2,30258509 2 -0,21640024
- -
0,76 0,87177979 0,27443685 2,30258509 2 0,207808603
- -
0,77 0,87749644 0,26136476 2,30258509 2 0,199207969

14
 - -
0,78 0,88317609 0,24846136 2,30258509 2 0,190598933
- -
0,79 0,88881944 0,23572233 2,30258509 2 0,181982063
- -
0,8 0,89442719 0,22314355 2,30258509 2 0,173357901
- -
0,81 0,9 0,21072103 2,30258509 2 0,164726966
- -
0,82 0,90553851 0,19845094 2,30258509 2 0,156089752
- -
0,83 0,91104336 0,18632958 2,30258509 2 0,147446733
- -
0,84 0,91651514 0,17435339 2,30258509 2 0,138798361
-
0,85 0,92195445 0,16251893 2,30258509 2 -0,13014507
- -
0,86 0,92736185 0,15082289 2,30258509 2 0,121487275
- -
0,87 0,93273791 0,13926207 2,30258509 2 0,112825371
-
0,88 0,93808315 0,12783337 2,30258509 2 -0,10415974
- -
0,89 0,94339811 0,11653382 2,30258509 2 0,095490744
- -
0,9 0,9486833 0,10536052 2,30258509 2 0,086818734
- -
0,91 0,9539392 0,09431068 2,30258509 2 0,078144043
- -
0,92 0,9591663 0,08338161 2,30258509 2 0,069466992
- -
0,93 0,96436508 0,07257069 2,30258509 2 0,060787887
-
0,94 0,96953597 -0,0618754 2,30258509 2 0,052107025
- -
0,95 0,97467943 0,05129329 2,30258509 2 0,043424688
- -
0,96 0,9797959 0,04082199 2,30258509 2 0,034741146
- -
0,97 0,98488578 0,03045921 2,30258509 2 0,026056662
- -
0,98 0,98994949 0,02020271 2,30258509 2 0,017371484
- -
0,99 0,99498744 0,01005034 2,30258509 2 0,008685853
1 1 0 2,30258509 1 0
-
15
 38,28873222
-
VALOR DE LA INTEGRAL CON 100 TRAPECIOS 0,191443661

Si asumimos que el valor de la integral con el límite inferior igual a 0,0000001 es el 100%
entonces el porcentaje de error con el valor -0,191443661 es 0,82%

-
VALOR DE LA INTEGRAL CON 1000 TRAPECIOS 0,192947622

Si asumimos que el valor de la integral con el límite inferior igual a 0,0000001 es el 100%
entonces el porcentaje de error con el valor -0.192947622 es 0,037%

16
Aporte 3. Ecuaciones Diferenciales Ordinarias

Se selecciona el desafío correspondiente al número de desafío seleccionado en las

actividades evaluativas anteriores.

1. Cálculo de la Ecuación Diferencial Ordinaria:

Desafío 2:
𝑥 ′ = 𝑐𝑜𝑠(𝑥), 𝑥(0) = 0

Desafío 3:
𝑡
𝑥′ = , 𝑥(0) = 1
1 + 𝑥2

Desafío 3: Calculo de la Ecuación Diferencial Ordinaria Método de Euler

𝑡
𝑥′ = , 𝑥(0) = 1
1 + 𝑥2
Con ℎ = 0.1. Realice los cálculos hasta 5ℎ

ECUACIÓN RAZÓN O EXPLICACIÓN

Desafío 3: Calculo de la Ecuación Diferencial Ordinaria Método de Runge – Kutta de orden 2

𝑡
𝑥′ = , 𝑥(0) = 1
1 + 𝑥2
Con ℎ = 0.1. Realice los cálculos hasta 5ℎ

17
 ECUACIÓN RAZÓN O EXPLICACIÓN

Desafío 3: Calculo de la Ecuación Diferencial Ordinaria Método de Runge – Kutta de orden 4

𝑡
𝑥′ = , 𝑥(0) = 1
1 + 𝑥2
Con ℎ = 0.1. Realice los cálculos hasta 5ℎ

ECUACIÓN RAZÓN O EXPLICACIÓN

2. Realice una gráfica comparando cada uno de los resultados encontrados en el

numeral anterior con la solución exacta y realice un análisis basándose en la teoría
revisada. Qué sucede si aumenta el tamaño de ℎ? qué sucede si disminuye el tamaño
de ℎ?

Desafío 4:
1
𝑥′ = , 𝑥(1) = 1
𝑥

18
Conclusiones

19
Referencias Bibliográficas

Alfa, T. (2012). Método de Euler [Video] Recuperado

de https://www.youtube.com/watch?v=F-dxL9Dnxc0&feature=youtu.be

Asesores. (2011). Método de Runge Kutta [Video] Recuperado

de https://www.youtube.com/watch?v=6bP6gcuLQoM&feature=youtu.be

Cetremo. (2014). Regla del Trapecio [Video] Recuperado de

https://www.youtube.com/watch?v=v0iIhdP9oxE&feature=youtu.be

Centremo. (2014). Regla de Simpson [Video] Recuperado

de https://www.youtube.com/watch?v=jJdp1n4vaGg

Chapra, S. C., & Canale, R. P. (2007). Métodos numéricos para ingenieros (5a. ed.). Pág.
719 – 756. Recuperado de https://ebookcentral-proquest-
com.bibliotecavirtual.unad.edu.co/lib/unadsp/reader.action?docID=4508648&ppg=745

Gómez, M. (27,03,2017). Métodos de Integración Numérica. [Archivo de video].

Recuperado de http://repository.unad.edu.co/handle/10596/11809

Nieves, H. A. (2014). Métodos numéricos: aplicados a la ingeniería: aplicados a la

ingeniería. México, D.F., MX: Larousse - Grupo Editorial Patria. Pág. 454 – 505.
Recuperado de https://ebookcentral-proquest-
com.bibliotecavirtual.unad.edu.co/lib/unadsp/reader.action?docID=3227640&ppg=471

Ríos, J.A. (2012). Integrales Múltiples [Video] Recuperado

de https://www.youtube.com/watch?v=hHJzmjkVVs8&feature=youtu.be

Vera, I. (2013). Integración de Romberg [Video] Recuperado

de https://www.youtube.com/watch?v=l44hSGuXYK8&feature=youtu.be

20