Todometodos PDF

MÉTODOS NUMÉRICOS
Fabián Ordóñez M.1

1 Universidad de las Fuerzas Armadas
ESPE
NRC: 6035-6036-6037-6045
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
Fabián Ordóñez M. (ESPE) Métodos Numéricos Mayo-Septiembre 2020 1 / 16
III Parcial:ECUACIONES DIFERENCIALES ORDINARIAS Método de Euler
ECUACIONES DIFERENCIALES
ORDINARIAS
MÉTODO DE EULER
Sea la EDO siguiente:

y ′ = f (x, y)
f (x) =
y(a) = y0
Para determinar una solución numérica aproximada a
esta EDO, en un intervalo [a, b], se debe construir un
conjunto de puntos xk , yk que son aproximaciones
numéricas a la solución. Para ello, se discretiza el
intervalo [a, b].
b−a
h=
n .
.
.
.
.
. . . . .
. . . .
. . . .
. . . .
. . . .
. . . . .
.
.
.
.
.
.
.
.
.

Las abscisas de los nodos equiespaciados que se

forman, se determinan con la relación: xk = a + kh con
k = 0, 1, · · · n. Aplicando el Teorema de Taylor con un
orden de aproximación de 2, a la función y(x) alrededor
de x0 , tenemos:
y ′′ (c)(x − x0 )2
′
y(x) = y(x0 ) + y (x − x0 ) +
2
La ecuación diferencial planteada, evaluada en x0 , es:
y ′ (x0 ) = f(x0 , y(x0 )) y si evaluamos el desarrollo de
Taylor para algún valor, por ejemplo x1 , tenemos:
y ′′ (c)(x1 − x0 )2
y(x1 ) = y(x0 ) + (x1 − x0 )f(x0 , y(x0 )) +
2
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
Además, se cumple que: x1 − x0 = h

0
z }| {
y ′′ (c)h2
y1 = y0 + h.f(x0 , y0 ) +
2
Al realizar el proceso iterativamente, se obtiene una
sucesión de valores, cuyos nodos formados e interpolados
representan la solución numérica aproximada de la EDO.
Al generalizar el método de Euler, éste viene dado
por:
y′
z }| {
yk+1 = yk + h. f(xk , yk )
con k = 0, 1, · · · (n − 1) y xk+1 = xk + h . . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
EJERCICIO Sea la EDO siguiente:

′ x−y
y = y
f (x) =
y(1) = 1
Determinar la solución numérica aproximada en [1; 1.5]
con h = 0,1, utilizando el Método de Euler.
yk+1 = yk + h.f(xk , yk )

′ xk −yk
Reemplazando y k = f(xk , yk ) por yk , tenemos:

xk − yk
yk+1 = yk + h.
yk
y2k + h(xk − yk )
yk+1 =
yk . . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
12 + (0,1)(1 − 1)
Si k = 0, y1 = =1
1
12 + (0,1)(1,1 − 1)
Si k = 1, y2 = = 1,01
1
1,012 + (0,1)(1,2 − 1,01)
Si k = 2, y3 = = 1,02881188
1,01
1,028811882 + (0,1)(1,3 − 1,02881188)
Si k = 3, y4 =
1,02881188
= 1,05517123
1,0551712322 + (0,1)(1,4 − 1,05517123)
Si k = 4, y5 =
1,05517123
= 1,087851117 .
.
.
.
.
. . . . .
. . . .
. . . .
. . . .
. . . .
. . . . .
.
.
.
.
.
.
.
.
.

k xk yk
0 1 1
1 1.1 1
2 1.2 1.01
3 1.3 1.02881188
4 1.4 1.05517123
5 1.5 1.087851117
Interpolando los nodos formados, obtenemos la siguiente
solución:
P5 (x) = 0,094814x5 −0,60067x4 +1,3088x3 −0,71931x2

− 0,61268x1 + 1,5291
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
1
La solución de la EDO, es: y(x) = 2x + 2x , que
contrastada con la ecuación numérica aproximada, se
puede apreciar en las siguientes gráficas.
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
III Parcial:ECUACIONES DIFERENCIALES ORDINARIAS Método de Heun
MÉTODO DE HEUN
Este método de mayor precisión, aplica primero el
teorema fundamental del cálculo para determinar el
primer punto P1 (x1 , y1 ), en el intervalo [x0 , x1 ].
Z x1 Z x1
f(x, y(x))dx = y ′ (x)dx = y(x1 ) − y(x0 )
x0 x0
Z x1
y(x1 ) = y(x0 ) + f(x, y(x))dx
x0
Si integramos utilizando la regla del trapecio simple, en
donde h = x1 − x0 , tenemos:
h
y(x1 ) = y(x0 ) + [f(x0 , y(x0 )) + f(x1 , y(x1 ))] dx
2 . . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
Ahora, se ha de reemplazar y(x1 ) por la aproximación

de Euler.
h
y(x1 ) = y(x0 ) + [f(x0 , y(x0 )) + f(x1 , hf(x0 , y(x0 ))]
2
Al realizar el proceso iterativamente, obtenemos una
sucesión de puntos, que al ser interpolados permiten
encontrar una solución aproximada a la EDO.
Generalizando, tenemos:
h
y(xk+1 ) = y(xk ) + [f(xk , y(xk )) + f(xk+1 , yk + hf(xk , y(xk ))]
2
Como f(xk ) = fk , el método Heun se presenta, así:
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
h
yk+1 = yk + [f(xk , yk ) + f(xk+1 , yk + hf(xk , yk )]
2
xk+1 = xk + h
Este método realiza predicción y corrección, con la
aproximación de Euler y luego con la regla del trapecio
EJERCICIO 2 Determinar la solución numérica
aproximada en [1; 1.5] con h = 0,1, utilizando el Método
de Heun. ′ x−y
y = y
f (x) =
y(1) = 1
Encontremos primero la aproximación de
Euler,
−y
reemplazando y ′k = f(xk , yk ) por kyk k .
x
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .

xk − yk y2k + h(xk − yk )
yk+1 = yk + h. =
yk yk
Ahora, reemplazando en la aproximación de Heun y
aplicado la funcion f a estos puntos, tenemos:
 y 
z }|k {
h y2 + h(xk − yk ) 
xk
 xk − yk z}|{ 
yk+1 = yk +  + f( xk+1 , k 
2 yk yk
 
y2 +h(x −y )
h  xk − yk xk+1 − k ykk k
= yk + + 
2 yk y2k +h(xk −yk )
yk
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .

h xk − yk xk+1 yk − y2k − h(xk − yk )
yk+1 = yk + +
2 yk y2k + h(xk − yk )
Con esta relación, podemos calcular los valores de yk .
Si k = 1,
[ ]
0,1 1−1 (1,1(1) − 12 − 0,1(1 − 1)
y1 = 1 + + = 1,005
2 1 12 + 0,1(1 − 1)
Si k = 2,
[ ]
0,1 1,1 − 1,005 1,2(1,005) − 1,0052 − 0,1(1,1 − 1,005)
y2 = 1,005 + +
2 1,005 1,0052 + 0,1(1,1 − 1,005)
= 1,018871558
Si k = 2,
[ ]
0,1 1,1 − 1,005 1,2(1,005) − 1,0052 − 0,1(1,1 − 1,005)
y2 = 1,005 + +
2 1,005 1,0052 + 0,1(1,1 − 1,005)
= 1,018871558
Si k = 3,
[ ]
0,1 1,1 − 1,018871558 1,3(1,018871558) − 1,0188715582 − 0,1(1,1 − 1,018871558)
y3 = 1,018871558 + +
2 1,018871558 1,0188715582 + 0,1(1,1 − 1,018871558)
. . . . . . . . . . . . . . . . . . . .
= 1,040462275
. . . . . . . . . . . . . . . . . . . .
Si k = 4,
[ ]
0,1 1,3 − 1,040462275 1,4(1,040462275) − 1,0404622752 − 0,1(1,3 − 1,040462275)
y4 = 1,040462275 + +
2 1,040462275 1,0404622752 + 0,1(1,3 − 1,040462275)
= 1,069257013
Si k = 5,
[ ]
0,1 1,4 − 1,069257013 1,5(1,069257013) − 1,0692570132 − 0,1(1,4 − 1,069257013)
y5 = 1,069257013 + +
2 1,069257013 1,0692570132 + 0,1(1,4 − 1,069257013)
= 1,101469574
Interpolando los nodos formados, obtenemos la

siguiente solución:
P5 (x) = 0,055446x5 − 0,32544x4 + 0,5703x3 + 0,1865x2

− 1,0575x1 + 1,5707
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
k xk yk
0 1 1
1 1.1 1.005
2 1.2 1.018957575
3 1.3 1.040760123
4 1.4 1.071669669
5 1.5 1.105543128
1
La solución de la EDO, es: y(x) = 2x + 2x , que
contrastada con la ecuación numérica aproximada
utilizando el método de Heun, se puede apreciar en las
siguientes gráficas.
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
MÉTODOS NUMÉRICOS

ESPE
NRC: 6035-6036-6037-6045
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
III Parcial:ECUACIONES DIFERENCIALES ORDINARIAS Método de diferencias finitas
ORDINARIAS
Consideremos la siguiente ecuación diferencial de
orden 2, resolver utilizando el método de diferencias
finitas.  ′′
 y − y = f(x)
f(x) = y(a) = y0

y(b) = yn
Para ello, primero se discretiza al intervalo [a, b] en n
subintervalos, así:
b−a
h=
n
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
Con este valor (salto), se construye las abscisas

equiespaciadas de los nodos a utilizar, con
k = 1, 2, · · · n − 1
xk = a + hk
La fórmula de diferencias finitas centrada para la segunda
derivada de orden 2, viene dada por la siguiente relación.
yk+1 − 2yk + yk−1
y′′k =
h2
A la relación planteada se presenta con una notación
adecuada y se reemplaza en la EDO, obteniendo la
siguiente expresión:
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
yk+1 − 2yk + yk−1

− yk = f(xk )
h2
Ordenando y operando, se tiene:
yk−1 − 2yk + yk+1 − h2 yk = h2 f(xk )
yk−1 + yk (−2 − h2 ) + yk+1 = h2 f(xk )
Luego, se ha de construir un sistema de (n-2) ecuaciones,
para cada valor posible que puede tomar k, así:
Si k = 1,
y0 − 2y1 + y2 − h2 y1 = h2 f(x1 )
y1 (−2 − h2 ) + y2 = h2 f(x1 ) − y0 (1)
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
Si k = 2,
y1 − 2y2 + y3 − h2 y2 = h2 f(x2 )
y1 + y2 (−2 − h2 ) + y3 = h2 f(x2 ) (2)
Si k = 3,
y2 − 2y3 + y4 − h2 y3 = h2 f(x3 )
y2 + y3 (−2 − h2 ) + y4 = h2 f(x3 ) (3)
... ··· ··· ...
Si k = n − 1,
yn−2 − 2yn−1 + yn − h2 yn−1 = h2 f(xn−1 )
yn−2 + yn−1 (−2 − h2 ) = h2 f(x2 ) − yn (n − 1)
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
El sistema planteado tiene (n-1) ecuaciones y va a

determinar (n-1) incógnitas, y1 , y2 , · · · yn−1 , puesto que
y0 y yn son conocidos, constituyen las condiciones de
frontera.
Estos datos, forman el vector columna que
representan las ordenadas de los nodos. Las abscisas, se
determinan con la ecuación de abscisas equiespaciadas.
Conocidos los nodos, que deben ser (n + 1)
intervalos, se puede formar los vectores x e y, con los que
se puede interpolar por cualquiera de los métodos
estudiados (Lagrange, Newton, mínimos cuadrados) y así
encontrar una solución numérica aproximada de la
ecuación diferencial planteada.
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
 
x0 = a
 x1 = x0 + h 
 
 x2 = x0 + 2h 
 
x=
 x. 3 = x0 + 3h.


 .. · · · .. 
 
 xn−1 = x0 + (n − 1)h 
xn = b
 
y0 = y(x0 )
 y1 
 
 y2 
 
y=  .y3 

 .. · · · ... 
 
 yn−1 
yn = y(xn ) .
.
.
.
.
. . . . .
. . . .
. . . .
. . . .
. . . .
. . . . .
.
.
.
.
.
.
.
.
.

EJERCICIO
Sea la EDO siguiente, encontrar una solución
numérica aproximada utilizando el método de las
diferencias finitas, para 8 subintervalos.


 y′′ − y = f(x)

y(0) = 0
f(x) =

 y(0,8) = 0,5 √

f(x) = 1 − cos( x)
0,8 − 0
h=
= 0,1
8
y0 = y(0) = 0, y8 = y(0,8) = 0,5 son los valores de
frontera. La ecuación que permite determinar el sistema
de ecuaciones para este ejercicio, viene dado por: .
.
.
.
.
. . . . .
. . . .
. . . .
. . . .
. . . .
. . . . .
.
.
.
.
.
.
.
.
.

yk−1 + yk (−2 − h2 ) + yk+1 = h2 f(xk )

Si k = 1,
y0 + y1 (−2 − h2 ) + y2 = h2 f(x1 )
y1 (−2 − (0,1)2 ) + y2 = (0,1)2 f(0,1) − 0 (1)
Si k = 2,
y1 + y2 (−2 − (0,1)2 ) + y3 = (0,1)2 f(0,2) (2)
Si k = 3,
y2 + y3 (−2 − (0,1)2 ) + y4 = (0,1)2 f(0,3) (3)
Si k = 4,
y3 + y4 (−2 − (0,1)2 ) + y5 = (0,1)2 f(0,4) (4)
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
Si k = 5,
y4 + y5 (−2 − (0,1)2 ) + y6 = (0,1)2 f(0,5) (5)
Si k = 6,
y5 + y6 (−2 − (0,1)2 ) + y7 = (0,1)2 f(0,6) (6)
Si k = 7,
y6 + y7 (−2 − (0,1)2 ) = (0,1)2 f(0,7) − 0,5 (7)
Al resolver este sistema de ecuaciones, obtendremos los

valores de las ordenadas yk , k = 1, · · · n − 1. Utilizando
el método de Gauss, para resolver este sistema de
ecuaciones, generan las siguientes matrices: . . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
  
−2,01 1 0 0 0 0 0 y1
 1 −2,01 1 0 0 0 0   y2 
   
 0 1 −2,01 1 0 0 0   y3 
   
 0 0 1 −2,01 1 0 0   y4 =
   
 0 0 0 1 −2,01 1 0   y5 
   
 0 0 0 0 1 −2,01 1   y6 
0 0 0 0 0 1 −2,01 y7
 
0, 04958472
 0, 09834440 
 
 0, 14628730 
 
 0, 19342159 
 
 0, 23975540 
 
 0, 28529681 
0, 33005381
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
Los vectores formados x e y, son:

   
x0 = 0,0 y0 = 0
 x1 = 0,1   y1 = −0,472285277 
   
 x2 = 0,2   y2 = −0,899708687 
   
 x3 = 0,3   y3 = −1,237784784 
   
x=  x 4 = 0,4  y =  y4 = −1,441951429
 


 x5 = 0,5   y5 = −1,467115998 
   
 x6 = 0,6   y6 = −1,267196327 
   
 x7 = 0,7   y7 = −0,794651801 
x8 = 0,8 y8 = 0, 5
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
Al interpolar estos puntos, obtenemos la solución

numérica aproximada a la EDO.
y(x) = 1240,08x8 −3472,214x7 +3993,044x6 −2430,173x5

+ 839,063x4 − 155,32x3 + 16,21x2 − 5,423x
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
Bibliografía
Bibliografía
[1] Nieves,A., Domínguez, F., Métodos Numéricos

aplicados a la ingeniería, Grupo Editorial Patria,
México, 2014
[2] Mathews, J. Kurtis, F., Métodos Numéricos con
Matlab, 3era ed., Prentice Hall, 2000
[3] Shoichiro, N., Métodos Numéricos aplicados con
software, 1era ed., Prentice-Hall Hispanoamericana
S.A. , 1992
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
MÉTODOS NUMÉRICOS

ESPE
NRC: 6035-6036-6037-6045
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
III Parcial:INTEGRACIÓN NUMÉRICA Métodos simples: Fórmulas de Newton - Codes
INTEGRACIÓN NUMÉRICA
Fórmulas de cuadratura cerrada
A los nodos equiespaciados los podemos representar de
la siguiente manera: xk = x0 + kh, k = 0, 1, 2, . . . , M, en
donde f(xk ) = fk .
Sea una función f continua en [a, b] y derivable por
lo menos hasta el orden indicado, las fórmulas de
cuadratura cerrada de Newton -Codes, vienen dadas por:
REGLA DEL TRAPECIO
∫ x1
h h3 f ′′ (c)
f(x)dx = [f(x + h) + f(x − h)] −
x0 2 12 . . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
REGLA DE SIMPSON (1/3)

∫ x2
h h5 f (4) (c)
f(x)dx = [f(x) + 4f(x + h) + f(x + 2h)] −
x0 3 90
Ordeon de
REGLA DE SIMPSON (3/8) aproximacion
∫ x3
3h
f(x)dx = [f(x) + 3f(x + h) + 3f(x + 2h) + f(x + 3h)]
x0 8
3h5 f (4) (c)
−
80
REGLA DE BOOLE
∫ x4
2h
f(x)dx = [7f(x) + 32f(x + h) + 12f(x + 2h)
x0 45 . . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
8h7 f (6) (c)

+32f(x + 3h) + 7f(x + 4h)] −
945
EJERCICIO 1
Demostrar la fórmula de cuadratura cerrada de
Simpson (1/3)
∫ x2
h
f(x)dx = [f0 + 4f1 + f2 ]
x0 3
DEMOSTRACIÓN
Con los puntos coordenados P0 (x0 , f(x0 )),
P1 (x1 , f(x1 )), y P2 (x2 , f(x2 )), construimos el polinomio
interpolador de Lagrange, así: .
.
.
.
.
. . . . .
. . . .
. . . .
. . . .
. . . .
. . . . .
.
.
.
.
.
.
.
.
.

(x − x1 )(x − x2 ) (x − x0 )(x − x2 )
P2 (x) = f(x0 ) + f(x1 )
(x0 − x1 )(x0 − x2 ) (x1 − x0 )(x1 − x2 )
(x − x0 )(x − x1 )
+ f(x2 )
(x2 − x0 )(x2 − x1 )
Si f(x) ≈ P2 (x) y f(xk ) = fk e integramos a ambos lados,
tenemos:
∫ x2 constant
∫ x2
(x − x1 )(x − x2 )
f(x)dx ≈ f0 dx+
x0 x0 (x0 − x1 )(x0 − x2 )
∫ x2 ∫ x2
(x − x0 )(x − x2 ) (x − x0 )(x − x1 )
f1 dx + f2 dx
x0 (x1 − x0 )(x1 − x2 ) x0 (x2 − x0 )(x2 − x1 ) . . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
Si se reemplaza x = x0 + ht de manera que se realiza

un cambio de variable y calculamos su derivada,
tenemos: dx = hdt. Como los nodos son equiespaciados, N
podemos reemplazar xi − xjD= (i − j)h y x − xi = (t − i)h
y los nuevos límites de integración, serán: t = 0 y t = 2.
∫ 2 ∫ 2
h(t − 1)h(t − 2)
f(x)dx ≈ f0 hdt+
0 0 (−h)(−2h)
∫ 2 ∫ 2
h(t)h(t − 2) h(t)h(t − 1)
f1 hdt + f2 hdt
0 (h)(−h) 0 (2h)(h)
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
∫ 2 ∫ 2
h
f(x)dx ≈ f0 (t2 − 3t + 2)dt−
0 2 0
∫ 2 ∫ 2
2 h
hf1 (t − 2t)dt + f2 (t2 − t)dt
0 2 0
∫ 2 ( ) ( ) ( )
h 2 −4 h 2
f(x)dx ≈ f0 − hf1 + f2
0 2 3 3 2 3
∫ 2
h
f(x)dx ≈ (f0 + 4f1 + f2 )
0 3
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
EJERCICIO 2 ( )
Dada la función f(x) = e−x sinh πx4 , encontrar la
integral numérica de cada fórmula de cuadratura
cerrada, en el intervalo de [1; 2].
h, se determina con la siguiente expresión:
b−a
h=
n−1
donde n es el número de nodos (n − 1) es el número de
subintervalos.
∫ TRAPECIO
x1
h
f(x)dx = [f0 + f1 ]
x0 2
2−1
El valor de h, es h = 2−1 = 1; entonces x0 = 1 y x1 = 2
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
∫ 2 ( πx )
−x h
e sinh dx = [f(1) + f(2)]
1 4 2
1
= [0,319566188 + 0,311446939] = 0,315506564
2
∫ x2
h
f(x)dx = [f0 + 4f1 + f2 )]
x0 3
2−1
El valor de h, es h = 3−1 = 12 ; entonces x0 = 1, x1 = 1
2 y
x2 = 2.
∫ 2 ( πx )
−x h
e sinh dx = [f(1) + 4f(3/2) + f(2)]
1 4 3
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
1
= [0,319566188 + 4(0,328037465) + 0,311446939]
6
= 0,323860498
∫ x3
3h
f(x)dx = [f0 + 3f1 + 3f2 + f3 ]
x0 8
El valor de h, es h = 2−1 1 4
4−1 = 3 ; entonces x0 = 1, x1 = 3 ,
x2 = 35 y x3 = 2.
∫ 2 ( πx )
−x 1
e sinh dx = [f(1) + 3f(4/3) + 3f(5/3) + f(2)]
1 4 8
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
1
= [0,319566188+3(0,329329579+3(0,324145214)+0,311446939]
8
= 0,323929688
REGLA DE BOOLE
∫ x4
2h
f(x)dx = [7f0 + 32f1 + 12f2 + 32f3 + 7f4 ]
x0 45
El valor de h, es h = 2−1 1 5
5−1 = 4 ; entonces x0 = 1, x1 = 4 ,
x2 = 23 , x3 = 74 y x4 = 2.
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
∫ 2 ( πx )
−x 1
e sinh dx = [7f(1) + 32f(5/4) + 12f(3/2)
1 4 90
+32f(7/4) + 7f(2)]
1
= [7(0,319566188)+32(0,328687078]+12(0,28037465)
90
+32(0,321473959) + 7(0,311446939)] = 0,323985496
El valor que podemos considerar exacto, es:
∫ 2 ( πx )
−x
e sinh dx = 0,323986088
1 4
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
III Parcial:INTEGRACIÓN NUMÉRICA Integración numérica compuesta
METODO COMPUESTO DEL

TRAPECIO
La integral simple de una función dada representa el área

del trapecio formado, lo que nos permite observar la
diferencia existente con el área real. Para encontrar una
área con mejor aproximación, se utiliza la regla .
.
.
.
.
. . . . .
. . . .
. . . .
. . . .
. . . .
. . . . .
.
.
.
.
.
.
.
.
.

b−a
h=
M (n-1)
El área total de la curva, es la suma de las áreas
parciales, así:
∫ b ∑ M
Atotal = f(x)dx = A1 + A2 + · · · + AM = Ak
a k=1
con K = 1, · · · , M
Ahora, cada área Ak viene dado por la integral simple:
∫ xk
h
Ak = f(x)dx = [fk−1 + fk ]
xk−1 2
Y el área total, será la suma de las áreas parciales. .
.
.
.
.
. . . . .
. . . .
. . . .
. . . .
. . . .
. . . . .
.
.
.
.
.
.
.
.
.

∫ ∑ h∑
b M M
h
f(x)dx = [fx(k−1) + fx(k) ] = [fk−1 + fk ]
a 2 2
k=1 k=1
Si desarrollamos la sumatoria, obtenemos:

∫ b
h
f(x)dx = [f0 +f1 +f−1+f2 +f2 +f3 +· · ·+fM−1 +fM ]
a 2
∫ b
h
f(x)dx = [f0 + 2f1 + 2f2 + · · · + 2fM−1 + fM ]
a 2
∫ b ∑
M−1
h
f(x)dx = [fa + fb ] + h fk
a 2
.
k=1
. . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
METODO COMPUESTO DE SIMPSON (1/3)
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
En este caso, se divide el intervalo [a, b] en 2M

subintervalos iguales, cuya dimensión, es:
b−a
h=
2M
De manera similar a la regla anterior, se deben sumar las
áreas parciales formadas y substituir cada área por la
fórmula de Simpson. Cada área Ak viene dado por la
integral simple de Simpson:
∫ xk
h
Ak = f(x)dx = [fk−2 + 4fk−1 + fk ]
xk−2 3
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
∫ b ∑
M
h
f(x)dx = [fx(2k−2) + 4fx(2k−1) + fx(2k) ]
a 3
k=1
∫
h∑
b M
f(x)dx = [f2k−2 + 4f2k−1 + f2 k ]
a 3
k=1
Si desarrollamos la sumatoria, obtenemos:
∫ b
f(x)dx = [f0 + 4f1 + f2 + f2 + 4f3 + f4
a
+ · · · + f2M−2 + 4f2M−1 + f2M ]
∫
2h ∑ 4h ∑
b M−1 M
h
f(x)dx = [f(a) + f(b)] + f2k + f2k−1
3 3 3
. . . . . . . . . . . . . . . . . . . .
a
Fabián Ordóñez M. (ESPE) Métodos Numéricos k=1
. . . . . . . . . . .
k=1
Mayo-Septiembre 2020
. . . . . . . .
18 / 28
.
METODO SIMPLE Y COMPUESTO DEL

RECTÁNGULO (PUNTO MEDIO)
En este método la integral en el intervalo [a, b] se

calcula como el área del rectángulo formado, es decir:
∫ b ( )
a+b
f(x)dx = (b − a)f
a 2 .
.
.
.
.
. . . . .
. . . .
. . . .
. . . .
. . . .
. . . . .
.
.
.
.
.
.
.
.
.

El método compuesto, divide al intervalo [a, b] en M

subintervalos iguales.
b−a
h= = xk − xk−1
M
Cada área Ak viene dado por la integral simple:
∫ xk ( )
xk−1 + xk
Ak = f(x)dx = (xk − xk−1 ) f
xk−1 2
Sumando todas las áreas, tenemos la integral buscada.
∫ b ∑
M ( )
xk−1 + xk
f(x)dx = (xk − xk−1 ) f
a 2
k=1
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
∫ b ∑
M ( )
xk−1 + xk
f(x)dx = h f
a 2
k=1
EJERCICIO 3 ∫2
Determinar la integral I = 1 (x2 .e−x + 1) en [1, 2],
para 7 nodos por los métodos compuestos del trapecio,
Simpson (1/3) y punto medio.
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
METODO TRAPECIO COMPUESTO

∫ b ∑ M−1
h
f(x)dx = [fa + fb ] + h fk
a 2
k=1
∫ 2
h
(x2 .e−x + 1)dx = [f0 + f6 ] + h[f1 + f2 + f3 + f4 + f5 ]
1 2
2−1 1
h= =
6 6
∫ 2
1
(x2 .e−x + 1)dx = [f(1) + f(2)]
1 12
1
+ [f(7/6) + f(8/6) + f(9/6) + f(10/6) + f(11/6)]
6 . . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
1 1
= [1,367879441 + 1,541341133] + [1,423854388+
12 6
1,468617134+1,50204286+1,524654452+1,537373591]
∫ 2
(x2 .e−x + 1)dx = 1,485192119
1
METODO SIMPSON COMPUESTO (1/3)
∫
2h ∑ 4h ∑
b M−1 M
h
f(x)dx = [f(a) + f(b)] + f2k + f2k−1
a 3 3 3
k=1 k=1
∫ 2
h 2h 4h
(x2 .e−x +1)dx = [f0 +f6 ]+ [f2 +f4 ]+ [f1 +f3 +f5 ]
1 2 3 3 . . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
∫ 2
1
(x2 .e−x + 1)dx = [f(1) + f(2)]
1 18
1 2
+ [f(8/6) + f(10/6)] + [f(7/6) + f(9/6) + f(11/6)]
9 9
∫ 2
1
(x2 .e−x + 1)dx = [1,367879441 + 1,541341133]
1 18
1
+ [1,468617134 + 1,524654452]
9
2
+ [1,423854388 + 1,50204286 + 1,537373591]
9
∫ 2
(x2 .e−x + 1)dx = 1,486653448 .
.
.
.
.
. . . . .
. . . .
. . . .
. . . .
. . . .
. . . . .
.
.
.
.
.
.
.
.
.
1
METODO COMPUESTO DEL PUNTO MEDIO

∫ b ∑
M ( )
xk−1 + xk
f(x)dx = h f
a 2
k=1
∫ 2 ∑
M ( )
2 −x xk−1 + xk
(x .e + 1)dx = h f
1 2
k=1
∫ 2 [ ( ) ( )
x 0 + x1 x1 + x2
(x2 .e−x + 1)dx = h f +f
2 2
(
1
) ( ) ( ) ( )
x2 + x3 x3 + x4 x4 + x5 x5 + x6
+f +f +f +f
2 2 2 2
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
∫ 2 [ ( ) ( )
2 −x1 1 + 7/6 7/6 + 8/6
(x .e +1)dx = f +f
6 2 2
(
1
) ( ) ( )
8/6 + 9/6 9/6 + 10/6 10/6 + 11/6
+f +f +f
2 2 2
( )
11/6 + 2
+f
2
∫ 2
1
(x2 .e−x + 1)dx = [1,397226784 + 1,447663745
1 6
+1,486726323+1,514649767+1,532182702+1,540375217]
∫ 2
(x2 .e−x + 1)dx = 1,486470756 .
.
.
.
.
. . . . .
. . . .
. . . .
. . . .
. . . .
. . . . .
.
.
.
.
.
.
.
.
.
1
ERROR ABSOLUTO
Ie − Itrap. = |1,486044373 − 1,485192119| = 0,000852254

Ie − IS = |1,486044373 − 1,486653448| = 0,000609075
Ie − IPM = |1,486044373 − 1,486470756| = 0,00042638
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
Bibliografía
Bibliografía

México, 2014
S.A. , 1992
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
MÉTODOS NUMÉRICOS

ESPE
NRC: 6035-6036-6037-6045
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
III Parcial: Derivación Numérica Definición, interpretación geométrica y análisis
DERIVACIÓN NUMÉRICA
La derivada de una función por incrementos, viene dado por:
f(x + h) − f(x)
f ′ (x) = lı́m
h→0 h
Gráficamente, tenemos:
La derivada en un determinado valor x0 , se realiza reemplazando x por x0 ,

así:
f(x0 + h) − f(x0 )
f ′ (x0 ) = lı́m
h→0 h . . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
Para determinar la derivada en cálculo numérico, se utiliza una sucesión

{an } = {Dn }, tal que hn → 0, en donde se realiza cálculos sucesivos para
valores de hn cada vez más cercanos a 0,así:
f(x + hn ) − f(x)
Dn =
hn
con n ∈ N. La precisión que se consiga al valor buscado depende entonces
del valor de hn , pues mientras más cercano a 0, la secante se acerca a la
tangente.
EJERCICIO 1
Sea la función f(x) = ln(x),se desea determinar la derivada en el punto
x = 1. En la tabla presentada a continuación, podemos ver el
comportamiento de la derivada a medida que h se hace más pequeño. La
función que se va a reemplazar, es:
ln(1 + hn ) − ln(1)
f ′ (1) =
hn
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
hn Dn
10−1 0.953102
10−2 0.9950331
10−3 0.99950033
10−4 0.9999500033
10−5 0.999995
10−10 1
Cuadro: Valores de la derivada de la función ln(x) con hn decreciente
EJERCICIO 2
√
Sea g(x) = x, encontrar la mejor aproximación para g ′ (2)
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
√ √
2+hn − 2
Los valores de hn , se reemplazan en: g ′ (2) = hn
hn Dn
10−1 0.349241122458
10−2 0.353112550269
10−3 0.353509207464
10−4 0.353548971286
10−5 0.353552948651
10−6 0.353553346377
10−7 0.353553384347
10−8 0.353553386567
10−9 0.353553408772
10−10 0.353552742638
√
Cuadro: Valores de la derivada de la función x con hn decreciente
Si se conoce que el valor de f ′ (2) = 0,3535533906, entonces la mejor

aproximación será cuando hn = 10−5 . Para determinar el mejor valor, se
puede aplicar la siguiente comparación: .
.
.
.
.
. . . . .
. . . .
. . . .
. . . .
. . . .
. . . . .
.
.
.
.
.
.
.
.
.

|Dn+1 − Dn | ≥ |Dn − Dn−1 |

Para el ejemplo, analicemos a D5 sería:
|D6 − D5 | ?? |D5 − D4 |
|0,353553346377−0,353552948651| |0,353552948651−0,353548971286|
3,97726x10−7 ≤ 3,977365x10−6
No cumple, analicemos para otro valor
|D10 − D9 | ?? |D9 − D8 |
|0,353552742638−0,353553408772| |0,353553408772−0,353553386567|
6,66134x10−7 ≥ 2,2205x10−8
Entonces, D9 es la mejor aproximación a la derivada de la función
analizada en dicho punto. . . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
III Parcial: Derivación Numérica Derivadas centradas, progresivas, regresivas
La derivada numérica de una función dada se puede determinar en forma

centrada, que se puede calcular cuando se conoce o se puede conocer
datos antes y después del valor buscado, progresiva cuando se conoce o se
puede conocer datos posteriores al buscado y regresiva si se conoce o se
puede conocer datos anteriores al buscado.
TEOREMA: Sea f una función continua en [a,b] y derivable hasta orden
3, en donde x-h y x + h pertenecen a [a, b]. Demostrar que la fórmula
centrada de la primera derivada de O(h2 ), viene dada por:
f(x + h) − f(x − h)
f ′ (x) = + 0(h2 )
2h
DEMOSTRACIÓN
f ′′ (x) ∗ h2 f (3) (c) ∗ h.32

f(x + h) = f(x) + f ′ (x) ∗ h + + (1)
2! 3!
f ′′ (x) ∗ h2 f (3) (c) ∗ h.3

2
f(x − h) = f(x) − f ′ (x) ∗ h + − (2)
2! .3! . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
f (3) (c) ∗ h2 f(3) (c) ∗ h2

f(x + h) − f(x − h) = 2f ′ (x) ∗ h + + (1) − (2)
3! 3!
f(x + h) − f(x − h) f(3) (c) ∗ h2
f ′ (x) = −
2h 6
TEOREMA:Sea f una función continua en [a,b] y derivable hasta orden 5,
en donde x-h, x + h, x + 2h y x-2h pertenecen a [a, b]. Demostrar que la
fórmula centrada de la primera derivada de O(h4 ), viene dada por:
f(x + 2h) + 8f(x + h) − 8f(x − h) + f(x − 2h)
f ′ (x) = − + 0(h4 )
12h
DEMOSTRACIÓN
f ′′ (x) ∗ h2 f (3) (x) ∗ h2 f (4) (x) ∗ h2 f (5) (c) ∗ h5

f(x+h) = f(x)+f ′ (x)∗h+ + + +
2! 3! 4! 5!
(3)
′′
f (x) ∗ h f (x) ∗ h f (x) ∗ h f (c) ∗ h5
2 (3) 2 (4) 2 (5)
f(x−h) = f(x)−f ′ (x)∗h+ − + −
2! 3! 4! . . 5!
. . . . . . . . . . . . . . . . . .
Fabián Ordóñez M. (ESPE) Métodos Numéricos

. . .
Mayo-Septiembre 2020
(4)
. . . .
8 / 20
. . . . . . . . . . . . .
4f ′′ (x) ∗ h2 8f (3) (x) ∗ h2 16f (4) (x) ∗ h2 32f (5

f(x+2h) = f(x)+2f ′ (x)∗h+ + + +
2! 3! 4!
(5)
′′ 2 2
4f (x) ∗ h 8f (x) ∗ h 16f (x) ∗ h 32f (5
(3) (4) 2
f(x−2h) = f(x)−2f ′ (x)∗h+ − + −
2! 3! 4!
(6)
Ahora, si observamos la fórmula que debemos demostrar, es necesario que
la ecuación (3) este multiplicada por el valor (8), la ecuación (4) por (-8)
y la ecuación (5) por (-1), tenemos:
8f ′′ (x) ∗ h2 8f (3) (x) ∗ h2 8f (4) (x) ∗ h2 8f (5)

8f(x+h) = 8f(x)+8f ′ (x)∗h+ + + +
2! 3! 4!
(7)
′′ 2
8f (x) ∗ h 8f (x) ∗ h 8f (x) ∗ h2 f
(3) 2 (4)
−8f(x−h) = −8f(x)+8f ′ (x)∗h− + − −
2! 3! 4!
(8)
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
4f ′′ (x) ∗ h2 8f(3) (x) ∗ h 2 16f (4) (x) ∗ h2 32

−f(x+2h) = −f(x)+2f ′ (x)∗h− + − +
2! 3! 4!
(9)
Sumando las ecuaciones (6), (7), (8) y (9), obtenemos:
−48f (5) (c)h5
−f(x + 2h) + 8(f(x + h) − 8f(x − h) + f(x − 2h) = 12f ′ (x)h +
120
f(x + 2h) + 8f(x + h) − 8f(x − h) + f(x − 2h) f (5) (c)h4
f ′ (x) = − +
12h 30
NOTACIÓN
f(x + h) = fk+1 = f1
f(x + 2h) = fk+2 = f2
f(x − h) = fk−1 = f−1
f(x − 2h) = fk−2 = f−2 . . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
En cálculos que utilizan la derivada, es muy común utilizar la primera,

segunda y de mayor orden, para lo cual se presenta las siguientes
relaciones:
f1 − f−1
f ′ (x0 ) = + O(h2 )
2h
f1 − 2f0 + f1
f ′′ (x0 ) = + O(h2 )
h2
f2 − 2f1 + 2f−1 − f−2
f (3) (x0 ) = + O(h2 )
2h3
f2 − 4f1 + 6f0 − 4f−1 + f−2
f (4) (x0 ) = + O(h2 )
2h3
−f2 + 8f1 − 8f−1 + f−2
f ′ (x0 ) = + O(h4 )
12h
−f2 + 16f1 − 30f0 + 16f−1 − f−2
f ′′ (x0 ) = + O(h4 )
12h2
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
Además, vamos a necesitar la primera, segunda o de orden superior

derivadas tanto progresivas como regresivas, según la información que se
tiene. Para ello, se indican las siguientes relaciones.
Progresivas
−3f0 + 4f1 − f2
f ′ (x0 ) = + O(h2 )
2h
2f0 − 5f1 + 4f2 − f3
f′′ (x0 ) = + O(h2 )
h2
Regresivas
3f0 − 4f−1 + f−2
f ′ (x0 ) = + O(h2 )
2h
2f0 − 5f−1 + 4f−2 − f−3
f ′′ (x0 ) = + O(h2 )
h2
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
TEOREMA Demostrar la segunda derivada progresiva de orden 2

2f0 − 5f1 + 4f2 − f3
f ′′ (x0 ) = + O(h2 )
h2
Esta demostración, vamos a realizar utilizando el polinomio interpolador
de Lagrange.
Al observar la fórmula que se va a demostrar, esta depende de 4 valores
f0 , f1 , f2 , f3 , lo que nos indica que debemos utilizar 4
nodos,P0 (x0 , y0 ); P1 (x1 , y1 ); P0 (x2 , y2 ); P0 (x3 , y3 ) y planteamos el
polinomio interpolador P3 (x).
(x − x1 )(x − x2 )(x − x3 ) (x − x0 )(x − x2 )(x − x3 )
P3 (x) = y0 + y1 +
(x0 − x1 )(x0 − x2 )(x0 − x3 ) (x1 − x0 )(x1 − x2 )(x1 − x3 )
(x − x0 )(x − x1 )(x − x3 ) (x − x0 )(x − x1 )(x − x2 )
+ y2 + y3
(x2 − x0 )(x2 − x1 )(x2 − x3 ) (x3 − x0 )(x3 − x1 )(x3 − x2 )
Como se desea determinar la 2 derivada, debemos entonces derivar 2
veces, además, sabemos que el polinomio aproxima a la función, lo que nos
permite reemplazar,así: .
.
.
.
.
. . . . .
. . . .
. . . .
. . . .
. . . .
. . . . .
.
.
.
.
.
.
.
.
.

(1)(x − x2 )(x − x3 ) + (x − x1 )(1)(x − x3 ) + (x − x1 )(x − x2 )(1)

f ′ (x) = y0 +
(x0 − x1 )(x0 − x2 )(x0 − x3 )
(1)(x − x2 )(x − x3 ) + (x − x0 )(1)(x − x3 ) + (x − x0 )(x − x2 )(1)

y1 +
(x1 − x0 )(x1 − x2 )(x1 − x3 )
(1)(x − x1 )(x − x3 ) + (x − x0 )(1)(x − x3 ) + (x − x0 )(x − x1 )(1)
y2 +
(x2 − x0 )(x2 − x1 )(x2 − x3 )
(1)(x − x1 )(x − x2 ) + (x − x0 )(1)(x − x2 ) + (x − x0 )(x − x1 )(1)
y3
(x3 − x0 )(x3 − x1 )(x3 − x2 )
(x − x2 ) + (x − x3 ) + (x − x1 ) + (x − x3 ) + (x − x1 ) + (x − x2 )
f ′′ (x) = y0 +
(x0 − x1 )(x0 − x2 )(x0 − x3 )
(x − x2 ) + (x − x3 ) + (x − x0 ) + (x − x3 ) + (x − x0 ) + (x − x2 )
y1 +
(x1 − x0 )(x1 − x2 )(x1 − x3 )
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
(x − x1 ) + (x − x3 ) + (x − x0 ) + (x − x3 ) + (x − x0 ) + (x − x1 )
y2 +
(x2 − x0 )(x2 − x1 )(x2 − x3 )
(x − x1 ) + (x − x2 ) + (x − x0 ) + (x − x2 ) + (x − x0 ) + (x − x1 )
y3
(x3 − x0 )(x3 − x1 )(x3 − x2 )
2[(x − x1 ) + (x − x2 ) + (x − x3 )]
f′′ (x) = y0 +
(x0 − x1 )(x0 − x2 )(x0 − x3 )
2[(x − x0 ) + (x − x2 ) + (x − x3 )]
y1 +
(x1 − x0 )(x1 − x2 )(x1 − x3 )
2[(x − x0 ) + (x − x2 ) + (x − x3 )]
y2 +
(x2 − x0 )(x2 − x1 )(x2 − x3 )
2[(x − x0 ) + (x − x1 ) + (x − x2 )]
y3
(x3 − x0 )(x3 − x1 )(x3 − x2 )
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
Reemplazando x = x0 y xi − xj = (i − j)h,yi = fi tenemos:
2[(−h) + (−2h) + (−3h)] 2[0 + (−2h) + (−3h]

f ′′ (x) = f0 + f1 +
(−h)(−2h)(−3h)] h(−h)(−2h)
2[0 + (−h) + (−3h)] 2[0 + (−h) + (−2h)]

f2 + f3
(2h)(h)(−h) (3h)(2h)(h)
f0 − 5f1 + 4f2 − f3
f′′ (x) =
h2
El orden de aproximación por este método no determina como lo realiza el
Teorema de Taylor, en el que se incluye en su desarrollo.
EJERCICIO 3
Utilizando la siguiente información, determinar la primera y segunda
derivada centrada, progresiva y regresiva de orden 2, en el punto x = 0,12
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
x sen(x)
0.09 0.08987855
0.10 0.09983342
0.11 0.10977830
0.12 0.11971221
0.13 0.12963414
0.14 0.13954311
0.15 0.14943813
CENTRADAS
f1 − f−1
f ′ (x0 ) =
2h
0,12963414 − 0,10977830
f ′ (0,12) = = 0,992792
2(0,01
f1 − 2f0 + f1
f ′′ (x0 ) =
h2
0,12963414 − 2(0,11971221) + 0,10977830
f ′′ (0,12) = = −0,1198
0,012 . . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
La primera y segunda derivada de la función dada, sus valoes considerados

exactos, son:
f ′ (0, 12) = cos(0, 12) = 0,992808664
f ′′ (0, 12) = −sin(0, 12) = −0,119712207
PROGRESIVAS
−3f0 + 4f1 − f2
f ′ (x0 ) = + O(h2 )
2h
−3(0,11971221) + 4(0,12963414) − 0,13954311
f ′ (0,12) = 2(0,01)
2(0,01)
f ′ (0,12) = 0,992841
2f0 − 5f1 + 4f2 − f3
f ′′ (x0 ) = + O(h2 )
h2
2(0,11971221) − 5(0,12963414) + 4(0,13954311) − 0,14943813
f ′′ (0,12) =
0,012
f ′′ (0,12) = −0,1197 . . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
REGRESIVAS
3f0 − 4f−1 + f−2
f ′ (x0 ) =
2h
3(0,11971221) − 4(0,10977830) + 0,09983342
f ′ (0,12) =
2(0,01)
f ′ (0,12) = 0,9928425
2f0 − 5f−1 + 4f−2 − f−3
f ′′ (x0 ) =
h2
2(0,11971221) − 5(0,10977830) + 4(0,09983342) − 0,08987855
f ′′ (0,12) =
0,012
f ′′ (0,12) = −0,1195
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
Bibliografía
Bibliografía
[1] Nieves,A., Domínguez, F., Métodos Numéricos aplicados a la

ingeniería, Grupo Editorial Patria, México, 2014
[2] Mathews, J. Kurtis, F., Métodos Numéricos con Matlab, 3era ed.,
Prentice Hall, 2000
[3] Shoichiro, N., Métodos Numéricos aplicados con software, 1era ed.,
Prentice-Hall Hispanoamericana S.A. , 1992
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
MÉTODOS NUMÉRICOS

ESPE
NRC: 6035-6036-6037-6045
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
III Parcial:ECUACIONES DIFERENCIALES ORDINARIAS Método de Taylor
ORDINARIAS
Consideremos la siguiente ecuación diferencial de orden
2, resolver utilizando el método de diferencias finitas.
{ ′
y = f(x, y)
f(x) =
y(a) = y0
Para ello, primero se discretiza al intervalo en el que
existe la solución [a, b] en n subintervalos, así:
b−a
h=
n
y se construye las abscisas equiespaciadas de los nodos a
utilizar, con k = 0, 1, 2, · · · n . . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
La solución aproximada y(x), se puede encontrar

utilizando el Teorema de Taylor, alrededor de x0 ,
tenemos:
(x − x0 )2 y ′′ (x0 )
y(x) = y(x0 ) + (x − x0 )y ′ (x0 ) + +
2!
(x − x0 )3 y (3) (x0 ) (x − x0 )n y (n) (x0 ) (x − x0 )n+1 y (n+1) (c)
+ +· · ·+ +
3! n! (n + 1)!
Si h = x − x0 y generalizamos en forma similar a los

métodos anteriores, tenemos:
h2 y ′′k h3 y k
(3) (n)
hn y k
yk+1 = yk + hy ′k + + + ··· + (1)
2! 3! n!
Ahora, se conoce que y ′ = f(x, y), entonces: . . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
d d ′ ∂f ∂f dy ∂f ∂f
y ′′ = f(x, y) = y = + = + f(x, y)
dx dx ∂x ∂y dx ∂x ∂y
d2 d ′′ ∂y ′′ ∂y ′′ dy ∂y ′′ ∂y ′′
y ′′′ = f(x, y) = y = + = + f(x, y)
dx2 dx ∂x ∂y dx ∂x ∂y
.. .. .. .. ..
. . . . .
d (n−1) d (n−1) ∂y (n−1) ∂y (n−1) dy

y (n) = f(x, y) = y = +
dx(n−1) dx ∂x ∂y dx
(n−1)
∂y ∂y (n−1)
= + f(x, y)
∂x ∂y
Reemplazando y ′′ , y ′′′ , · · · y(n−1) en (1), obtenemos los

yk , k = 0, 1, · · · , n − 1 que permiten construir los
respectivos nodos. .
.
.
.
.
. . . . .
. . . .
. . . .
. . . .
. . . .
. . . . .
.
.
.
.
.
.
.
.
.

Con los nodos formados, se puede interpolar para

determinar la solución aproximada a la EDO planteada.
El número de derivadas que utiliza el método nos da la
precisión de la aproximación, se conoce como orden (n =
2,3,4, ...)
EJERCICIO Dada la EDO siguiente:
{ ′ x−y
y = y
f (x) =
y(1) = 1
Determinar la solución numérica aproximada en [1; 1.5]

con h = 0,1, utilizando el Método de Taylor, de orden 2.
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
1,5 − 1
0,1 =
n
n=5
h2 y ′′k
yk+1 = yk + +hy ′k
2!
xk − yk
y ′k = f(xk , yk ) =
yk
′′ ∂f ∂f 1 xk xk − yk y2k − x2k + xk yk
y = + f(x, y) = − 2 . =
∂x ∂y yk yk yk y3k
Luego, se calcula los valores que de yK para cada valor
de xk , k = 0, 1, · · · , n − 1
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
( ) ( 2 )
xk − yk h2 yk − x2k + xk yk
yk+1 = yk + h +
yk 2 y3k
Si k = 0,
( ) ( 2 )
1−1 0,12 1 − 12 + 1(1)
y1 = 1 + 0,1 + = 1,005
1 2 13
Si k = 1
( ) ( )
1,1 − 1,005 0,12 1,0052 − 1,12 + 1,1(1,005)
y2 = 1,005 + 0,1 +
1,005 2 1,0053
= 1,01891312
Si k = 2,
( ) ( )
1,2 − 1,01891312 0,12 1,018913122 − 1,22 + 1,2(1,01891312)
y3 = 1,01891312 + 0,1 +
1,01891312 2 1,018913123
= 1,02665621
Si k = 3,
( ) ( )
1,3 − 1,02665621 0,12 1,026656212 − 1,32 + 1,3(1,02665621)
y4 = 1,02665621 + 0,1 +
1,02665621 2 1,026656213
= 1,056509153
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
Si k = 4,
( ) ( )
1,4 − 1,056509153 0,12 1,0565091532 − 1,42 + 1,4(1,056509153)
y5 = 1,056509153 + 0,1 +
1,056509153 2 1,0565091533
= 1,092027735
Los nodos formados, se resumen en la siguiente tabla.
k xk yk
0 1 1
1 1.1 1.005
2 1.2 1.01891312
3 1.3 1.02665621
4 1.4 1.056509153
5 1.5 1.092027735
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
Al interpolar estos puntos, obtenemos la solución

numérica aproximada a la EDO.
y(x) = −73,406 ∗ x5 + 458,5 ∗ x4 − 1139 ∗ x3

+ 1406,9 ∗ x2 − 864,11 ∗ x1 + 212,11
La solución exacta de la EDO, es:

x 1
y(x) = +
2 2x
Al contrastar con la ecuación numérica aproximada
utilizando el método de Taylor, se puede apreciar en las
siguientes gráficas.
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
Bibliografía
Bibliografía

México, 2014
S.A. , 1992
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .

Todometodos PDF

Cargado por

Copyright:

Formatos disponibles

Todometodos PDF

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Todometodos PDF

Cargado por

Copyright:

Formatos disponibles

MÉTODOS NUMÉRICOS

Fabián Ordóñez M.1

Fabián Ordóñez M. (ESPE) Métodos Numéricos Mayo-Septiembre 2020 2 / 16

Las abscisas de los nodos equiespaciados que se

Además, se cumple que: x1 − x0 = h

EJERCICIO Sea la EDO siguiente:

Fabián Ordóñez M. (ESPE) Métodos Numéricos Mayo-Septiembre 2020 6 / 16

P5 (x) = 0,094814x5 −0,60067x4 +1,3088x3 −0,71931x2

Ahora, se ha de reemplazar y(x1 ) por la aproximación

Interpolando los nodos formados, obtenemos la

P5 (x) = 0,055446x5 − 0,32544x4 + 0,5703x3 + 0,1865x2

Fabián Ordóñez M.1

Con este valor (salto), se construye las abscisas

yk+1 − 2yk + yk−1

El sistema planteado tiene (n-1) ecuaciones y va a

Fabián Ordóñez M. (ESPE) Métodos Numéricos Mayo-Septiembre 2020 7 / 15

Fabián Ordóñez M. (ESPE) Métodos Numéricos Mayo-Septiembre 2020 8 / 15

yk−1 + yk (−2 − h2 ) + yk+1 = h2 f(xk )

Al resolver este sistema de ecuaciones, obtendremos los

Los vectores formados x e y, son:

Al interpolar estos puntos, obtenemos la solución

y(x) = 1240,08x8 −3472,214x7 +3993,044x6 −2430,173x5

[1] Nieves,A., Domínguez, F., Métodos Numéricos

Fabián Ordóñez M.1

REGLA DE SIMPSON (1/3)

8h7 f (6) (c)

Fabián Ordóñez M. (ESPE) Métodos Numéricos Mayo-Septiembre 2020 4 / 28

Si se reemplaza x = x0 + ht de manera que se realiza

METODO COMPUESTO DEL

La integral simple de una función dada representa el área

Fabián Ordóñez M. (ESPE) Métodos Numéricos Mayo-Septiembre 2020 13 / 28

Fabián Ordóñez M. (ESPE) Métodos Numéricos Mayo-Septiembre 2020 14 / 28

Si desarrollamos la sumatoria, obtenemos:

METODO COMPUESTO DE SIMPSON (1/3)

En este caso, se divide el intervalo [a, b] en 2M

METODO SIMPLE Y COMPUESTO DEL

En este método la integral en el intervalo [a, b] se

Fabián Ordóñez M. (ESPE) Métodos Numéricos Mayo-Septiembre 2020 19 / 28

El método compuesto, divide al intervalo [a, b] en M

METODO TRAPECIO COMPUESTO

METODO COMPUESTO DEL PUNTO MEDIO

Ie − Itrap. = |1,486044373 − 1,485192119| = 0,000852254

[1] Nieves,A., Domínguez, F., Métodos Numéricos

Fabián Ordóñez M.1

La derivada en un determinado valor x0 , se realiza reemplazando x por x0 ,

Para determinar la derivada en cálculo numérico, se utiliza una sucesión

Si se conoce que el valor de f ′ (2) = 0,3535533906, entonces la mejor

Fabián Ordóñez M. (ESPE) Métodos Numéricos Mayo-Septiembre 2020 5 / 20

|Dn+1 − Dn | ≥ |Dn − Dn−1 |

La derivada numérica de una función dada se puede determinar en forma

f ′′ (x) ∗ h2 f (3) (c) ∗ h.32

f ′′ (x) ∗ h2 f (3) (c) ∗ h.3

f (3) (c) ∗ h2 f(3) (c) ∗ h2

f ′′ (x) ∗ h2 f (3) (x) ∗ h2 f (4) (x) ∗ h2 f (5) (c) ∗ h5

Fabián Ordóñez M. (ESPE) Métodos Numéricos

4f ′′ (x) ∗ h2 8f (3) (x) ∗ h2 16f (4) (x) ∗ h2 32f (5

8f ′′ (x) ∗ h2 8f (3) (x) ∗ h2 8f (4) (x) ∗ h2 8f (5)

4f ′′ (x) ∗ h2 8f(3) (x) ∗ h 2 16f (4) (x) ∗ h2 32

En cálculos que utilizan la derivada, es muy común utilizar la primera,