TAREA 3 – DIFERENCIACIÓN E INTEGRACIÓN NUMÉRICA Y EDO

PRESENTADO POR

PRESENTADO A:

UNIVERSIDAD NACIONAL ABIERTA Y A DISTANCIA - UNAD

METODOS NUMERICOS

MAYO DE 2020
 Introducción

En el presenta trabajo se realizará la solución de una serie de problemas asociados

a los temas de integración numérica, ecuaciones diferenciales y diferenciación
numérica, con el objetivo cumplir con las actividades estipuladas para la actividad
número 3 en la asignatura métodos numéricos, además de realizar reforzar los
conocimientos que se presentan en el área.
 Tema 1. Diferenciación numérica
Desafío 4:

Dada la función f ( x )=sen ( πx /3 ), encuentre f ' ( 1 ) usando los esquemas de

diferencias finitas hacia adelante, atrás y central, con h=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.
La integración numérica es una técnica de análisis numérico para calcular una
aproximación a la derivada de una función en un punto utilizando los valores y
propiedades de la misma.
Por definición la derivada de una función es:

f ( x +h )−f ( x)
f ' ( x )=lim fffff
h→ 0 h

Las aproximaciones numéricas que podamos hacer (para h > 0) serán:

Diferencias hacia adelante:

Diferencias hacia atrás:

Diferencias centrales:

La solución analítica es:

f ( x )=sen ( πx /3 ), encuentre f ' ( 1 )

f ' ( 1 )=cos ( π3. x ) . π3 = 16 π=0,52399

f(x) 0,7431 0,8660 0,9511
X 0,8 1 1,2

a) Diferencias finitas progresivas:

Se desea conocer el valor aproximado de f ' (1), por ende, se requiere determinar los

valores de f ( x 0 ) y f ( x 0 +h ), sabiendo que x 0=1 y h=0.2

f ( x 0 ) f ( 1 ) =0,8660

f ( x 0 +h ) f ( 1,2 )=0,9511

'
f ( x 0 +h )−f ( x0 ) 0,9511−0,866
f ( x 0 )= = =0,4252
h 0,2

Error Porcentual= |0,5236−0,4252

0,5236 |∗100=19 %
b) Diferencias finitas centradas:

Se desea conocer el valor aproximado de f ' (1), por ende, por ende, se requiere

determinar los valores de f ( x 0 −h ) y f ( x 0+ h ), sabiendo que x 0=1 y h=0.2

f ( x 0 −h ) f ( 0,8 ) =0,7431

f ( x 0 +h ) f ( 1,2 )=0,9511

'
f ( x 0 +h )−f ( x0 ) 0,9511−(0,7431)
f ( x 0 )= = =0,5198
h 2∗0,2

Error Porcentual= |0,5236−0,5198

0,5236 |∗100=1 %

C) Diferencias finitas regresivas:

Se desea conocer el valor aproximado de f ' (1), por ende, se requiere determinar los

valores de f ( x 0 ) y f ( x 0−h ), sabiendo que x 0=1 y h=0.2

f ( x 0 ) f ( 1 ) =0,866

f ( x 0 −h ) f ( 0,8 ) =0,7431

'
f ( x 0 ) −f ( x 0−h ) 0,866−(0,7431)
f ( x 0 )= = =0,6144
h 0,2

Por medio de los 3 métodos, se obtuvo que los resultados derivados a través de las
diferencias hacia arriba y hacia abajo son consideradamente lejanas al valor exacto,
presentando un error del 19% y 17% respectivamente. El método que resulta ser
más exacto es el de las diferencias centradas dando un valor que genera un error
del 1%.

La selección del valor de h no es la mejor opción, para este caso se recomienda que
sea menor que h<0,2 para que la solución sea adecuada para cada uno de los
métodos realizados. De realizar lo contrario, el error incrementará.

Tema 2. Integración numérica

Problema 4:
π /3

∫ √ sen(x )dx
0

Valor exacto = 0,6865911344756129

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

Tabla 1. Componentes X vs f(x)

0,174 0,52359 0,698131 0,872 1,047

x 0 5 0,34906585 9 7 7 2
0,416 0,58482488 0,70710 0,801740 0,875 0,930
f(x) 0 7 3 7 4 2 6

a) Regla del trapecio simple y compuesta

La primera aproximación para la integración numérica es la regla trapezoidal. La

curva a integrar se aproxima con trapezoides. Cada trapezoide posee un área
proporcional al promedio de la altura de los lados multiplicado por la base
(espesor). Como se muestra a continuación:

Fig. 1 Regla Trapezoidal.

Se expresa de la siguiente manera:

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

Intervalo [ 0 , π /3 ] a=0 ; b=π /3

n=6

π
−0
b−a 3
h= = =π /18
n 6

A través de los valores de la tabla 1, se procede a sustituir en la ecuación

obteniendo que:

π /3

∫ f ( x ) dx ≅ I = π /18
2
[ 0+2∗( 0,4167+ 0,5848+0,7071+0,8017+0,8752 ) +0,9306 ]
0

I =0.6721

Error Porcentual

Error Porcentual %= |0,6866−0,6721

0,6866 |∗100=2 %
b) Regla de Simpson1/3 simple y compuesta

Este método se aproxima a un polinomio de grado 2, el cual se expresa de la

siguiente manera:

n−2 n−2
b

a
h
∫ f ( x ) dx ≅ I = 3 ∗¿ ( 2 2

)
f ( a)+ 4 ∑ f ( x2 k +1)+2 ∑ f ( x 2 k )+ f (b) ¿
k=0 k=1
 Fig. 2 Regla de Simpson 1/3.

Intervalo [ 0 , π /3 ] a=0 ; b=π /3

n=6

b−a π /3−(0)
h= = =π /18
n 6

A través de los valores de la tabla 1, se procede a sustituir en la ecuación

obteniendo que:

π /3

∫ f ( x ) dx ≅ I = π /18
3
[ 0+4∗( 0,4167+ 0,7071+ 0,8752 )+ 2(0,5848+0,8017)+0,9306 ]
0

I =0 .6807

|0,6866−0,6807
Error Porcentual %=
0,6866 |∗100=1 %
c) Regla de Simpson3/8 simple y compuesta

f se aproxima por un polinomio de grado 3, el cual se expresa de la siguiente

manera:

∫ f ( x ) dx ≅ I = 38h ∗¿ ¿ ¿
a
 Fig. 3 Regla de Simpson 3/8.

Intervalo [ 0 , π /3 ] a=0 ; b=π /3

n=6

b−a π /3−(0)
h= = =π /18
n 6

A través de los valores de la tabla 1, se procede a sustituir en la ecuación

obteniendo que:

π /3

∫ f ( x ) dx ≅ I = 3∗π8/ 18 [ 0+3∗( 0,4167+ 0,5848+0,7071+0,8017+0,8752 ) +0,9306 ]

0

I =0,7257

Error Porcentual %= |0,6866−0,7257

0,6866 |∗100=6 %
Por medio de los métodos realizados, se puede observar que el método más exacto
para la función en estudio es el de la Regla de Simpsons 1/3 el cual genera un error
del 1% con respecto al valor exacto de la integral.

Tema 3. Ecuaciones Diferenciales Ordinarias

Calcule la Ecuación Diferencial Ordinaria:

Problema 4:
x 2 ( )
()
3
x ' ( t )=( tx ) −
t
, x 1 =1

1. Empleando cada uno de los siguientes métodos:

a) Método de Euler

El Método de Euler o de las Tangentes constituye el primer y más sencillo ejemplo

de método numérico para la resolución de un problema de valor inicial.

y ' =f ( x , y ) , y ( x 0 )= y 0

donde suponemos además que se verifican las hipótesis del Teorema de Picard , y
en consecuencia existe solución única para el problema.

Interpretando la E.D.O. y’ = f(x, y) como un campo de direcciones en el plano x − y

y la condición inicial y(x0) = y0 como un punto (x0, y0) de dicho plano, podemos
aproximar la función solución y(x) por medio de la recta tangente a la misma que
pasa por ese punto:

y ' ≅ y 0+ f ( x 0 , y 0)(x−x 0 )

donde se ha utilizado que la pendiente de dicha tangente es: m = y’ (x0) y, en

consecuencia: m = f(x0, y0).

Calculamos así de manera aproximada el valor de la solución y en el punto de

abscisa x1 como:

y ( x 1) ≅ y 1= y 0 +f ( x0 , y 0)( x 1−x 0 )

y con este punto (aproximado) ya calculado, podemos repetir el método para

obtener otro punto aproximado (x2, y2) de la forma:

y ( x 2) ≅ y 2= y 1+ f (x 1 , y1 )( x 2−x 1)

y así sucesivamente. Es habitual en este método tomar abscisas equiespaciadas, es

decir, calcular la solución aproximada en puntos de la forma: xn = xn−1 + h = x0 +
nh, siendo h el paso del método. De esta forma se obtienen las fórmulas que nos
determinan la solución aproximada en la forma:

x n=x n−1 +h y n= y n−1 + f ( x n−1 , y n−1 ) h

Desde el punto de vista geométrico, tenemos en definitiva que el Método de Euler
aproxima a la función solución por medio de una línea poligonal, la aproximación
será tanto peor cuanto mayor sea en número de pasos, es decir, cuanto más “lejos”
nos encontremos del punto inicial (x0, y0). Por otro lado, el error será
evidentemente tanto mayor cuanto más grande sea el “paso” del método, h.

Para el ejercicio seleccionado se utilizará h=0,1. Se conoce que x (1)=1. Entonces

x i=x i−1+ h

t i= y i−1+ hf ( xi −1 , y i−1)

x t f(x,t)
1 1 0
1,1 1,0000 0,121
1,2 1,0121 0,38571401
1,3 1,0507 1,01726267
1,4 1,1524 2,72355749
1,5 1,4248 8,65255142

b) Método de Runge – Kutta de orden 2

La primera opción que podemos aplicar es integrar mediante el método de los

trapecios, es decir tomando:
xn+ 1
1
∫ f ( x , y ( x ) ) dx ≅ 2 h ( f ( x n , y n ) +f ( x n+ 1 , y n+1 ) )
xn

Al ser desconocida yn+1 en la expresión anterior, lo aproximaremos por ¯yn+1,

donde y¯n+1 es la estimación de yn+1 que resultaría aplicando el método de Euler.
Tendremos así:
xn+ 1
1
∫ f ( x , y ( x ) ) dx ≅ 2 h ( f ( x n , y n ) +f ( x n+ 1 , ý n+1 ) )
xn
Con

ý n+ 1= y n +hf ( x n , y n )

y llegaremos a la expresión del método:

h
y n+ 1= y n + f ( x , y ) + f ( xn +1 , ý n +1) )
2( n n
Lo normal es presentar el método con las expresiones siguientes:

k 1=hf ( x n , y n ) ; k 2=hf ( xn +1 , y n +k 1)

1
y n+ 1= y n + ( k 1 +k 2)
2

x t f(xn,yn) k1 f(xn+1;yn+k1) K2
1 1 0 0 0,121 0,0121
1,1 1,00605 0,1598133 0,01598133 0,46615843 0,04661584
1,2 1,03734859 0,59076091 0,05907609 1,489967981 0,1489968
1,3 1,14138503 1,96958424 0,19695842 5,483632629 0,54836326
1,4 1,51404587 8,66857797 0,8668578 45,15422386 4,51542239
1,5 4,20518597 250,847147

C) Método de Runge – Kutta de orden 4

Ahora se introduce un nuevo paso intermedio en la evaluación de la derivada. Una

vez más se presentan varias opciones en la evaluación y es posible ajustar de tal
manera que se garantice el error local de manera proporcional a h 5 (es decir
garantizando exactitud en el cuarto orden en el polinomio de Taylor), lo cual lleva a
un error global proporcional a h4. El Método de cuarto orden más habitual es el
determinado por las fórmulas siguientes:

x i+1=x i +h
 h
t i+1=t i + ( k 1+2 k 2+2 k 3 + k 4 )
6
donde
k 1=f (x i , t i )
h h
k 2=f ( xi + ,t i+ k 1 )
2 2
h h
k 3=f (xi + ,t i + k 2 )
2 2
k 4=f ( x i +h , t i +k 3 h)

x+h/ x+h/
x t k1 2 t+k1h/2 k2 2 t+k2h/2 k3 x+h t+k3h k4
0,0551 0,0707 1,0070 0,1664
1 1 0 1,05 1 3 1,051,00275625 8 1,1 8 1
1,0069701 1,01525 0,3601 1,0429 0,6368
1,1 8 0,16572 1,15 6 0,3085 1,15 1,022395085 6 1,2 9 1
1,0426344 0,63393 1,07433 1,0680 1,2709 1,1697 2,2811
1,2 3 2 1,25 1 7 1,25 1,096037808 4 1,3 3 6
1,1691862 1,28294 4,0881 5,4104 1,7102 13,056
1,3 9 2,27512 1,35 2 7 1,35 1,373595024 9 1,4 4 1
1,7413292 13,8422 43,575 183,51 20,092 27377,
1,4 1 2 1,45 2,43344 4 1,45 3,920101247 5 1,5 9 8
465,83853 3,41E+0 1705932 1,8E+2 2,9E+6 2,9E+6
1,5 9 8 1,55 7 2 1,55 9,2438E+20 3 1,6 2 1E+188

D) 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 h ? qué sucede si
disminuye el tamaño de h ?
 Euler Runge-Kutta Orden 2 Runge-Kutta Orden 4
500
400
300

t
200
100
0
0.9 1 1.1 1.2 1.3 1.4 1.5 1.6
x

Como se puede observar por medio de los métodos realizados, se obtienen los
puntos que satisfacen la ecuación diferencial bajo la condición inicial x(1)=1.

La diferencia tan grande que existe en cada método radica principalmente en el

valor h exigido para determinar los valores. Si se utiliza un h<0,1 se podrá observar
que los resultados son muy cercanos.

Para visualizar lo anterior mencionado, se presenta una gráfica utilizando un valor

de h=0,01

Euler Runge-Kutta Orden 2 Runge-Kutta Orden 4

1
1
1
1
t

1
1
1
0.99 1 1.01 1.02 1.03 1.04 1.05 1.06
x

Se observa que los métodos presentan una solución más cercana.

Conclusión

A través de los problemas modelados, se pudieron concluir los siguientes aspectos:

 En la diferenciación numérica, para el caso en estudio. El método que
presenta el valor más cercano con la solución exacta, es el método de
diferencias finitas centradas. El cual presenta un error del 1%.
 Para la integración numérica se tiene que el valor que mejor se aproxima al
valor exacto de la integral, es el método de la regla de Simpson 1/3,
generando un error del 1% con respecto a la solución exacta.
 Para la E.D.O. se tiene que entre los 3 métodos una diferencia considerable.
Para este caso, se recomienda utilizar un valor de h<1 para que los métodos
presenten una convergencia adecuada.

Referencias Bibliográficas
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 – 467.

Chapra, S. C., & Canale, R. P. (2007). Métodos numéricos para ingenieros (5a. ed.).
Pág. 719 – 730.