MÉTODOS NUMÉRICOS

FASE 5: TRABAJO COLABORATIVO 3

APORTES PARTE INDIVIDUAL

LUZ MARY PADILLA GREY

CÓDIGO: 1047428718

WILLIAM SALAZAR

TUTOR

GRUPO:

100401_62

UNIVERSIDAD NACIONAL ABIERTA Y A DISTANCIA – UNAD

FACULTAD DE CIENCIAS E INGENIERÍA

PROGRAMA INGENIERÍA DE ALIMENTOS

2018
 INTRODUCCION

En el siguiente trabajo se resuelve derivadas e integrales de funciones,

usando los métodos iterativos para elegir el intervalo adecuado,
desarrollando el conocimiento de los métodos numéricos; esto con el fin de
identificar, comprender y aplicar los diferentes métodos necesarios para
realizar la solución de los problemas; Además utilizara los diferentes
métodos iterativos buscando soluciones a ecuaciones diferenciales en valor
inicial o de valor en la frontera, comparando los procesos y resultados
arrojados por éstos.

Se utilizan métodos como la Regla Del Trapecio, Regla De Simpson,

Integración De Romberg, Integrales Multiples, Método De Euler, Métodos De
Tylor Y Método De Runge Kutta adquiriendo los conocimientos necesarios
debido a la web conferencia y a los textos presentados en el entorno de
conocimiento en el Curso De Métodos Numéricos En La Universidad Nacional
Abierta Y A Distancia UNAD
 APORTE N° 1

1. Plantee y solucione dos ejercicios sobre Diferenciación Numérica

explicando paso a paso el procedimiento utilizado.

Solución

1.1. Considere la tabla adjunta, en la misma se muestran la posición

para determinado tiempo de una partícula moviéndose en el
espacio, calcule de ser posible la aceleración de la misma. Si la
partícula describe un movimiento parabólico dado por f ( t )=t 2, calcule
el error de la aceleración aproximada.
j ti xi
0 0 0
1 0.1 0.001
2 0.2 0.005
3 0.3 0.008
4 0.4 0.018
5 0.5 0.027

La segunda derivada en su forma aproximada es solo aplicable a puntos

interiores de los dados, es decir que podemos calcular la aceleración en
t=0.1 , 0.2 ,0.3 y 0.4
1
f ' ' ( 0.1 ) ≈ [ 0.005−2 ( 0.001 ) +0 ] ≈ 2
0.12
1
f ' ' ( 0.2 ) ≈ [ 0.008−2 ( 0.005 )+ 0.001 ] ≈ 2
0.12
1
f ' ' ( 0.3 ) ≈ [ 0.018−2 ( 0.008 )+ 0.005 ] ≈2
0.12
1
f ' ' ( 0.4 ) ≈ [ 0.027−2 ( 0.018 ) +0.008 ] ≈ 2
0.12

La segunda derivada de t es igual a 2 en todo el intervalo por lo que este

ejercicio tuvimos un error de cero.
1.2. Use aproximaciones con diferencias finitas hacia adelante y hacia
atrás 0(h) y una aproximación de diferencia centrada de 0(h 2) para
estimar la primera derivada de:
f ( x )=−0.1 x 4 −0.15 x3 −0.5 x2 −0.25 x+1.2
En x=0.5 utilizando un incremento en h=0.5, repetimos el cálculo con
h=0.25 y observamos que la derivada se calcula directamente como:

f ' ( x )=−0.4 x 3 −0.45 x 2−1.0 x−0.25

Y se puede utilizar para calcular el valor verdadero como:

f ' ( 0.5 ) =−0.9125

Para h=0.5 la función se emplea para determinar y evaluamos en:

f ( x )=−0.1 x 4 −0.15 x3 −0.5 x2 −0.25 x+1.2

x i−1=0 f ( x i−1 )=1.2

x i=0.5 f ( x i ) =0.925

x i+1=1.0 f ( xi +1 )=0.2

Con estos valores podemos calcular la diferencia divididas hacia adelante:

f ( x¿¿ i )
f ' (x ¿¿i)=f ( xi +1 )− +O(x i+1−x i )¿ ¿
xi +1−x i

0.2−0.925
f ' ( 0.5 ) ≅ =−1.45
1−0.5

|aproximacionaproximacion
ε t=
real−aproximacion anterior
real |100 %
0.9125−1.45
ε =|
0.9125 |
t 100 %

|ε t|=58.9 %
La diferencia dividida hacia atrás
f (x ¿ ¿ i−1) ∇ f 1
f ' (x ¿¿i)=f ( xi ) − = ¿¿
h h

0.925−1.2
f ' ( 0.5 ) ≅ =−0.55
0.5
 ε t= |aproximacionaproximacion
real−aproximacion anterior
real |100 %

ε t= |0.9125−0.55
0.9125 |
100 %

|ε t|=39.7 %

La diferencia dividida centrada

f ( x¿¿ i−1)
f ' (x ¿¿i)=f ( xi +1 )− −0(h2)¿ ¿
2h

0.2−1.2
f ' ( 0.5 ) ≅ =−1.0
2 ( 0.5 )

ε t= |aproximacionaproximacion
real−aproximacion anterior
real |100 %

ε t= |0.9125−1
0.9125 |
100 %

|ε t|=9.6 %
Para h=0.25, la función se emplea para determinar
x i−1=0.25 f ( x i−1 )=1.10351563

x i=0.5 f ( x i ) =0.925

x i+1=0.75 f ( xi +1 )=0.63632813

Esos valores sirven para calcular las diferencias divididas hacia adelante
 f ( x¿ ¿i)
f ' (x ¿¿i)=f ( xi +1 )− −0 ( x i+1 ) −xi ¿ ¿
( x i+1 )−x i

0.63632813−0.925
f ' ( 0.5 ) ≅ =−1.155
0.25

ε t= |aproximacionaproximacion
real−aproximacion anterior
real |100 %

ε t= |0.9125−1.155
0.9125 |100 %
|ε t|=26.5 %

La diferencia dividida hacia atrás

f (x ¿ ¿ i−1) ∇ f 1
f ' (x ¿¿i)=f ( xi ) − = ¿¿
h h

0.925−1.10351563
f ' ( 0.5 ) ≅ =−0.714
0.25

ε t= |aproximacionaproximacion
real−aproximacion anterior
real |100 %

ε t= |0.9125−0.714
0.9125 |100 %
|ε t|=21.7 %

La diferencia dividida centrada

 f ( x¿¿ i−1)
f ' (x ¿¿i)=f ( xi +1 )− −0(h2)¿ ¿
2h

0.63632813−1.10351563
f ' ( 0.5 ) ≅ =−0.934
0.5

ε t= |aproximacionaproximacion
real−aproximacion anterior
real |100 %

ε t= |0.9125−0.934
0.9125 |100 %
|ε t|=2.4 %

Para ambos tamaños de paso, la aproximación en diferencias centrales es

más exacta que las diferencias hacia adelante y hacia atrás

2. Solucione el siguiente ejercicio utilizando la Regla del Trapecio. (n= 4).

2
x3
a. ∫ 1/ 2
dx
1 3+ x

Solución

Para solucionar el problema se utilizan las siguientes formulas:

b
h
I =∫ f ( x ) dx= ¿ ¿
a 2

a=1 ,b=2 , n=4

x n=x n−1 +h

x 0=a=1

x 1=1+0.25=1.25

x 2=1.25+0.25=1.5
 x 3=1.5+0.25=1.75

x 4 =1.75+ 0.25=2

Luego,
2−1
h= =0.25
4
2 3
∫ 3+xx 1/ 2 dx= h2 ¿
1

Reemplazando en la función, se tiene:

13
f 1=
( ) 1
=0.25
3+1 2

1. 253
f ( 1.25 ) = 1
=0.474285
2
3+ 1.25

1. 53
f ( 1.5 )= 1
=0.798864
2
3+ 1.5

1. 753
f ( 1.75 ) = 1
=1.239770
2
3+ 1.75

23
f ( 2) = 1
=1.812327
2
3+ 2

Reemplazando en la formula, se tiene:

2
x3 0.25
I =∫ 1
dx= [ 0.25+2 ( 0.474285 ) +2 ( 0.798864 ) +2 ( 1.239770 )+ 1.812327 ]
1 2
2
3+ x
I =0.886020
4

b. ∫ √3 x ( e2 x ) dx
2

Solución
 b
h
I =∫ f ( x ) dx= ¿ ¿
a 2

b−a 4−2
h= = =0.5
n 4
a=2 , b=4 , n=4
x n=x n−1 +h

x 0=a=2

x 1=2+0.5=2.5

x 2=2.5+0.5=3

x 3=3+0.5=3.5

x 4 =3.5+0.5=4

Remplazando en la función, se tiene

3
f ' ( 2 ) =√ 2(e 2 2 )=68.789358
( )

3
f ( 2.5 )=√ 2.5 ( e2 ( 2.5) )=201.427646
3
f ' (3)=√3 ( e2 (3 ) )=581.845004
3
f ' ( 3.5)= √ 3.5 ( e 2 (3.5 )) =1665.012077

f ' ( 4 ) =√3 4 (e 2( 4 ))=4731.975844

Reemplazando en la formula, se tiene:

4
3 0.5
I =∫ √ x ( e2 x ) dx= [68.789358+ 2 ( 201.427646 ) +2 ( 581.845004 ) +2 ( 1665.012077 ) +4731.975844 ]
2 2

I =2424.33

3. Soluciones los siguientes ejercicios utilizando la Regla de Simpson 1/3

y 3/8. (n= 4)
4 3
ex
a .∫ dx b . ∫ e 2 x ln ( 2 x ) dx
2 x−1 1
Solución

a. Se utilizan las siguientes fórmulas para Simpson 1/3

b

∫ f ( x ) ≈ h3 ¿ ¿
a

b−a 4−2
h= = =0.5
n 4
b=4 , a=2 , n=4
x n=x n−1 +h

x 0=2

x 1=2+0.5=2.5

x 2=2.5+0.5=3

x 3=3+0.5=3.5

x 4 =3.5+0.5=4

Ahora reemplazamos en la función:

e2
f (x 0)=f ( 2 )= =7.389056
2−1

e 2.5
f ( x 1 ) =f ( 2.5 )= =8.121662
2.5−1

e3
f ( x 2 ) =f ( 3 )= =10.042768
3−1

e 3.5
f ( x 3 ) =f ( 3.5 )= =13.246180
3.5−1

e4
f ( x 4 )=f ( 4 )= =18.199383
4−1

Remplazamos los valores obtenidos en la primera fórmula así:

4
ex
∫ x−1 dx = 0.5
3
[ 7.389056+4 ( 8.121662 )+2 ( 10.042768 ) +4 ( 13.246180 )+ 18.199383 ]
2

4
ex
∫ x−1 dx =21.857557
2
 3

b. ∫ e2 x ln ( 2 x ) dx por Simpson 3/8

1

Es importante notar que en el método de Simpson 3/8, el número de

intervalos “n” solo puede ser un múltiplo de 3, ósea 3, 6, 9, etc. Se puede
verificar por medio de la gráfica siguiente donde se grafican la función y
sus intervalos.
 APORTE N° 2

4. Solucione los siguientes ejercicios utilizando la Integración de

Romberg. Usando segmentos de longitud 1, 1/2 y 1/4.
3 2
3

a .∫ e dx b .∫ e x ln ( x ) dx
x

1 1
Solución

a.
3
3

∫ e x dx
1

Primero calculamos las integrales del nivel 1 usando la regla del trapecio
para las longitudes de segmentos indicados:
h1 =1
1 3

1
h2 =
2 1 2 3

1
h3 =
4 3
1 1.5 2 2.5

Primer nivel: regla del trapecio

h
I= ¿
2
 3−1
( )
1 3 3

I ( h1 )=
2
[ e 1 + e3 ] =5.3204824060
3−1
I ( h )=
( 2 ) 13 23 33
[ e + 2e +e ]=2.6602412328
2
2

3−1
I ( h )=
( 4 )
¿
3
2

Segundo nivel de aproximación

4 1
I ( h 2) − I ( h1 )
3 3

I ( h1 ) → integral menos exacta

I ( h2 )=integral mas exacta

Por medio de diagrama calculamos:

4 1
5.3204824060 ( 2.662412328 )− ( 5.324824060 )=1.7734941750
3 3
4 1
2.6602412328 ( 1.3301511532 )− ( 2.6602412328 ) =0.8867877933
3 3
1.3301511532

Tercer nivel de aproximación

16 1
I − I
15 m 15 L
I m → Integralmas exacta

I L → Integralmenos exacta

16 1
( 0.8867877933 ) − (1.7734941750 )=0.8276740345
15 15
3
3

∫ e x dx ≈ 0.8276740345
1

b.
 2

∫ e x ln ( x ) dx
1

Calculamos las integrales de nivel 1 usando:

h1 =1
1 2

1
h2 =
2 1 1.5 2

1
h3 =
4 2
1 1.25 1.5 1.75

Con esto tenemos:

2−1
I ( h )=
( 1 ) 1 2
[ e ln ( 1 ) +e ln ( 2 ) ]=2.5608517009
1
2

2−1
I ( h )=
( 2 )
¿
2
2

2−1
I ( h )=
( 4 )
¿
3
2

Segundo nivel de aproximación:

4 1
I ( h 2) − I ( h1 )
3 3

I ( h1 ) → integral menos exacta

I ( h2 )=integral mas exacta

Por medio de diagrama calculamos:

4 1
2.5608517009 ( 2.1890101221 )− ( 2.5608517009 )=2.0650629291
3 3
4 1
2.1890101221 ( 1.8995966887 )− ( 2.1890101221 )=1.8031255442
3 3
1.8995966887

Tercer nivel de aproximación

16 1
I m− I L
15 15
I m → Integralmas exacta

I L → Integralmenos exacta

16 1
( 2.650629291 )− ( 1.8031255442 )=2.7071295408
15 15
2

∫ e x ln ( x ) dx ≈ 2.7071295408
1

5. Solucione los siguientes ejercicios de Integrales Múltiples compruebe

que:

a.
2
0.5 x

∫ ∫ e y /x dy dx ≈ 0.0333054
0.1 x 3
b.

1 2x

∫∫ ( x ¿ ¿ 2+ y 3)dy dx ≈ 1.000122 ¿
0 x

Solución

a.
2
0.5 x

∫ ∫ e y /x dy dx ≈ 0.0333054
0.1 x 3

Se calcula primero la integral

2
x y

∫ e x dy
3
x

Aplicamos integración por sustitución:

 y du 1
u= = xdu=dy
x dy x

Calculamos la integral indefinida:

y

∫ eu xdu=x ∫ e u du=x eu=e x +c

Calculamos los límites:
y
( ) 3

lim x e x =x e x / x =xe x
3
2

y →x

2
y x

2
( )
lim x e =x e =xe x x x

y →x

2
x y
2

∫ e x dy=x e x −x e x
3
x

Calculamos la integral
0.5
2

∫ ( x e x−x e x ) dx
0.1

Calculamos la integral indefinida:

2 2

∫ (¿ x e x−x e x )dx=∫ x e x dx−∫ x e x dx ¿

∫ x e x dx → se integra por partes
du
u=x ; =1 ; du=dx ; dv =e x dx ; ∫ dv=∫ e x dx ; v =e x
dx

∫ udv=uv−∫ vdu
∫ x e x dx=x e x −∫ e x dx
∫ x e x dx=x e x−¿ e x ¿
Ahora:
2

∫ x e x dx
Primero aplicamos integración por sustitución:
du du
u=x2 ; =2 x ; =dx
dx 2x
 x eu
∫ 2 x du= 12 ∫ e u= 12 e u= 12 e x
2

2
1 2

∫ (¿ x e x −x e x )dx= x e x −e x − 2 e x + c ¿
Calculamos los límites:
1 2
1 2

lim (¿ x e x −e x − e x )=(0.1)e0.1 −e0.1 − e (0.1) =−1.49967890¿

x→ 0.1 2 2
1 2
1 2

lim (¿ x e x −e x − e x )=(0.5) e 0.5−e 0.5− e(0.5 ) =−1.46637334 ¿

x→ 0.5 2 2
2
0.5 x y
x
∫∫e dy dx=−1.46637334− (−1.4996789 )
0.1 x 3
2
0.5 x

∫ ∫ e y /x dy dx ≈ 0.0333054
0.1 x 3
b.
1 2x

∫∫ ( x ¿ ¿ 2+ y 3)dy dx ≈ 1.000122 ¿
0 x

Calculamos la integral indefinida:

4
∫ (x ¿ ¿ 2+ y3 ) dy=∫ x 2 dy+∫ y3 dy =x2 y + y4 + c ¿
Ahora calculamos los límites:

y4 x4 x4
y →x
(
lim x 2 y +
4 )
=x 2 ( x ) + =x 3+
4 4

y4 2 x4 16 x 4
lim x 2 y+
y →2 x
( 4 )
=x 2 ( 2 x )+
4
=2 x3 +
4

2x 4
x4
∫ x 2 + y 3 dy=2 x 3+ 4 x 4 − x 3 + x4 ( )
=¿ 2 x 3+ 4 x 4 −x3 −
4
¿
x

3 4 x4
¿ x +4 x −
4

3 16 x 4 −x 4
¿x +
4
 4 x 3 +16 x 4−x 4
¿
4

4 x 3 +15 x 4
¿
4

Ahora calculamos la integral:

1
4 x3 +15 x 4
∫ 4 dx
0

1
¿ ¿
4

4∗x 4 15∗x 5
¿
1
4 [( ) ( )]
4
+
5

1
¿ ( x 4 +3 x 5 ) +c
4

Calculamos los límites:

1 1
lim ( (¿ x 4 +3 x 5 ))= ¿ ¿
x →0 4 4
1 1
lim ( (¿ x 4 +3 x 5 ))= ¿ ¿
x →1 4 4
1 2x

∫∫ (x ¿ ¿ 2+ y 3) dy dx=1−0=1 ¿
0 x

6. Demostrar que el valor aproximado en x = 0,2 de la solución del

problema de valor inicial 𝑦’ = 𝑥 + 𝑦 , 𝑦(0) = 0 usando el Método de
Euler con h = 0.05 y Zo = 0 , es 0,016
Solución

y ' =x + y h=0.05

Los modos son:

x 0=0 ; x 1=0.05 ; x 2=0.1 ; x3 =0.15 ; x 4=0.2

El método de Euler es:

y 0= y a ; y j+1= y j hf ( x j , y j ) j=0.1 , … n−1
En este caso:
y 0=0; y j+1= y j +0.05 ( x j , y j ) j=0,1,2,3

Iteraciones

Fase 0
x 0=0 ; y 0= y ( x 0 ) =0

Fase 1
x 0=0 ; y 0=0 ; y 1= y 0 +0.05 ( x 0+ y 0 )=0

y 1=0+ 0.05 ( 0+0 )=0

Fase 2
x 1=0.05 ; y 1=0 ; y 2= y 1 +0.05 ( x 1 + y 1 )

y 2=0+ 0.05 ( 0.05+0 )

y 2=0.0025

Fase 3
x 2=0.1 ; y 2=0.0025 ; y 3= y 2 +0.05 ( x 2+ y 2 )

y 3=0.0025+ 0.05 ( 0.1+ 0.0025 )

y 3=0.007625

Fase 4
x 3=0.15 ; y 3=0.007625 ; y 4 = y 3 +0.05 ( x 3+ y 3 )

y 4 =0.007625+0.05 ( 0.15+0.007625 )

y 4 =0.01550625

Resumimos en una tabla:

j xj yj
0 0 0
1 0.05 0
2 0.1 0.0025
3 0.15 0.007625
4 0.2 0.01550625
 APORTE N°3

7. Aplicar el método de Taylor de orden dos a la ecuación y’= Cos (xy),

con la condición inicial: y(0) = 1. Utilizar h = 0.5
Solución

Aplicamos la siguiente expresión

1
y ( x n+1 ) ≈ y n+1= y n + y ' ( x n ) h+ y '' ( x n) h 2
2

Tenemos que:

y ' ( x n ) =f ( x , y n )=cos ( xn y n )
''
Calculamos y ( x n )

d ( cos ( xy ))
y' '=
dx

Aplicamos la derivada de un producto

y ' ' =−sen ( xy ) y−sen ( xy ) y '

y ' ' =−sen ( xy ) y−sen ( xy ) cos ( xy )=−sen ( xy ) ( y+ cos ( xy ) )

x 1=x 0+ h=0.5

h2
y 1= y 0 +hcos ( x 0 y 0 ) + ( −sen ( x0 y 0 ) )( y 0+ cos ( x 0 y 0 ) )
2

( 0.5 )2
y 1=1+0.5 cos ( 0 )+ (−sen ( 0∗1 ) ) ( 1+cos ( 0∗1 ) )
2
y 1=1.5
x 2=x 1+ h=0.5+ 0.5=1

h2
y 2= y 1 +hcos ( x1 y 1 ) + ¿
2

( 0.5 )2
y 2=1.5+ 0.5cos ( 0.5∗1.5 )+ ¿
2

y 2=1.99581

8. Plantee y solucione paso a paso un ejercicio por el Método de Runge-

Kutta de cuarto orden.

Solución

Con el método de Runge-Kutta de orden 4 obtener una aproximación del valor de

y(1,5) para el siguiente problema de valor inicial, tomando un paso h = 0,1

y ' =2 ty
{ y ( 1 ) =1

El primer paso para resolver este problema es determinar la malla de puntos en

donde se va a obtener la solución.

Como en este caso h está dado, se tiene que N = (1,5 - 1)/0,1 = 5.

Por lo tanto, los puntos en donde se va a determinar la solución, dados por la

fórmula ti = 1 + 0,1 i, para i =1,2,3,4,5, son:
t1 = 1,1

t2 = 1,2

t3 = 1,3

t4 = 1,4

t5 = 1,5

Una vez establecida la malla del problema, tenemos, para i = 0:

t 0=1 , y 0=1

k 1=hf ( t 0 , y 0 )=0.1∗2 t 0 y 0=0.2

1 1 1 1
( ) ( )
k 2=hf t 0 + 0.1; y 0+ k 1 =0.1 f 1+ 0.1 ; 1+ 0.2 =0.1 f (1,2,1,1 ) =0.231
2 2 2 2
 1 1 1 1
( ) ( )
k 3=hf t 0 + 0.1 ; y 0+ k 2 =0.1 f 1+ 0.1 , ; 1+ 0.231 =0.1 f ( 1,2 ; 0.11 155 )=0.234255
2 2 2 2

k 4=hf ( t 0 +0.1 ; y 0 +k 3 ) =0.1 f ( 1+ 0.1; 1+0.2 34255 )=0.1 f (1,2 ; 1.234255 )=¿ 0.715361

Resolviendo tenemos:
1 1
y 1= y 0 + ( k 1 +2 k 2 +2 k 3 +k 4 )=1+ ¿
6 6
y 1=1.23367435

Y aplicando sucesivamente la fórmula de Runge-Kutta de orden 4, para i desde 1

hasta 4, se obtienen los datos que se muestran en la siguiente tabla, donde
además se muestra el valor de la solución exacta para cada punto de la malla.

ti yi Solución exacta Error absoluto

1.00 1.0000 1.0000 0.0000
1.10 1.2337 1.2337 0.0000
1.20 1.5527 1.5527 0.0000
1.30 1.9937 1.9937 0.0000
1.40 2.6116 2.617 0.0001
1.50 3.4902 3.4904 0.0002

Al analizar la tabla anterior y comparar los resultados obtenidos con el método

Runge-Kutta de orden 4 con los valores reales, se ve por qué es tan difundido
este método. En la próxima tabla se comparan los métodos de Euler y Runge
Kutta de orden 4 para el mismo problema.
ti Euler Runge Kutta Solución exacta
1.00 1.0000 1.0000 1.0000
1.10 1.2000 1.2337 1.2337
1.20 1.4640 1.5527 1.5527
1.30 1.8154 1.9937 1.9937
1.40 2.2874 2.6116 2.617
1.50 2.9278 3.4902 3.4904
 CONCLUSIONES

La regla del trapecio es una de las más utilizadas en el cálculo de

aproximaciones de integrales definidas.
La regla de Simpson 1/3 es la mas exacta en 3 orden
Cada ejercicio fue sencillo de resolverlo ya que nos dicen con que
método realizarlo y solo fue leer la teoría con la formula.