UNIVERSIDAD ESTATAL

PENINSULA DE SANTA ELENA

FACULTAD DE CIENCIAS DE LA
INGENIERÍA

CARRERA DE PETRÓLEOS

ASIGNATURA:
ANALISIS NUMERICO

DOCENTE:
ING. JORGE PAUL ORDOÑEZ Ph.D.

ESTUDIANTE:
PAOLA DOILE GONZALEZ
THALIA RODRIGUEZ MUÑIZ

PERIODO ACADÉMICO:
2022-1
 METODO DE EULER
Fórmula de Euler
Las siguientes fórmulas constituyen los métodos clásicos para resolver numéricamente
ecuaciones diferenciales ordinarias. Son equivalentes a usar varios términos de la serie
de Taylor, pero sustituyen las derivadas por aproximaciones, de tal manera que no se
requiere especificarlas o usar el método de derivación implícita anterior. Las fórmulas y
los algoritmos resultantes son independientes de una EDO particular. Sea una ecuación
diferencial ordinaria explícita de primer orden con una condición en el inicio:

F (x, y) = y´(x), y (𝒙𝟎 ) = 𝒚𝟎

ℎ2
𝑦𝑖+1 = 𝑦𝑖 + ℎ𝑦´𝑖 + 𝑦′′(𝑧)
2!
ℎ2
= 𝑦𝑖 + ℎ𝑓(𝑥𝑖, 𝑦𝑖 ) + 𝑦 ′′ (𝑧), 𝑥𝑖 ≤ 𝑧 ≤ 𝑥𝑖+1
2!
Definición: Fórmula de Euler

Algoritmo para calcular puntos de la solución de una EDO de primer orden con la fórmula de
Euler.

1) Defina f(x,y) y la condición inicial (x0, y0)

Ejemplo 1. Obtenga dos puntos de la solución de la siguiente ecuación diferencial con

i=0: y1 = y0 + h f(x0, y0) = 1 + 0.1 f(0, 1) = 1 + 0.1(1 + 0 - 02 + 1) = 1.2000;

x1 = x0 + h = 0 + 0.1 = 0.1

i=1: y2 = y1 + h f(x1, y1) = 1.2 + 0.1 f(0.1, 1.2) = 1.2 + 0.1 (1.2 + 0.1 - 0.12 + 1] = 1.4290 x2 = x1 +
h = 0.1 + 0.1 = 0.2

Para comprobar comparamos con la solución exacta: y(x) = x + x2 + ex

y(0.1) = 1.2152

y(0.2) = 1.4614

El error es significativo. Para reducirlo se pudiera reducir h. Esto haría que el error de
truncamiento se reduzca, pero si la cantidad de cálculos es muy grande, pudiera acumular
error de redondeo. Una mejor estrategia es usar métodos más precisos que no requieran que
h sea muy pequeño

El método de Euler es la técnica de aproximación más básica para resolver problemas

𝒅𝒚
= 𝒇(𝒙, 𝒚), 𝒂 ≤ 𝒕 ≤ 𝒃, 𝒚(𝒂) = 𝜶.
𝒅𝒙
No se obtendrá una aproximación continua a la solución 𝒚(𝒕); en su lugar, las
aproximaciones para 𝒚 se generan en varios valores, llamados puntos de malla, en el
intervalo [𝒂,𝒃].
Ejemplo 2.
Dada la siguiente ecuación diferencial con la condición inicial:

Solución Numérica
Aplicamos el método de Euler y para ello, observamos que la distancia entre y
no es lo suficientemente pequeña. Si dividimos esta distancia entre cinco obtenemos
un valor de y por lo tanto, obtendremos la aproximación deseada en cinco pasos.
Sustituyendo estos datos en la fórmula de Euler, tenemos, en un primer paso:

De esta forma, tenemos los siguientes datos:

Y así sucesivamente hasta obtener Resumimos los resultados en la siguiente tabla

N 𝑥𝑛 𝑦𝑛
0 0 1
1 0.1 1
2 0.2 1.02
3 0.3 1.0608
4 0.4 1.12445
5 0.5 1.2144
Concluimos que el valor aproximado, usando el método de Euler es:

Puesto que, en este caso, conocemos el valor verdadero, podemos usarlo para calcular
el error relativo porcentual que se cometió al aplicar la fórmula de Euler .
Ejemplo 3. MÉTODO DE EULER MEJORADO

Aplicar el metodo de euler mejorado, para aproximar y(0.5) si:

𝒚′ = 𝟐𝒙𝒚
y(0) =1
Solucion
Vemos que este es el mismo ejemplo 1 delmetodo anterior, asi que definimos h = 0.1 y
encontraremos la aproximacion despues de cinco iteraciones. A diferencia del metodo
de Euler 1, en cada iteraciones requerimos de dos calculos en vez de uno solo: el de
𝑦′𝑛 primero y posteriormente el de 𝑦𝑛 para aclarar el metodo veamos ccon detalle las
primeras dos iteraciones. Primero que nada, aclaramos que tenemos los siguientes
datos iniciales

En nuestra primera iteración tenemos: Segunda iteración:

Resultado de la tabla: N 𝒙𝒏 𝒚𝒏
0 0 1
1 0.1 1.01
2 0.2 1.040704
3 0.3 1.093988
4 0.4 1.173192
5 0.5 1.28336
Concluimos entonces que la aproximación obtenida con el método de Euler mejorado
es:

Con fines de comparación, calculamos el error relativo verdadero:

Fórmula de Heun
Ejemplo 4.

Algoritmo para resolver una EDO de primer orden con la fórmula de Heun
1) Defina f(x, y) y la condición incial (x0, y0)
2) Defina h y la cantidad de puntos a calcular m
3) Para i = 1, 2, ..., m
4) K1 = hf(xi, yi)
5) K2 = hf(xi + h, yi + K1))

6) yi+1 = yi +1/2 (K1 + K2) .

7) xi+1 = xi + h
8) fin
Ejemplo. Obtener dos puntos de la solución de la siguiente ecuación diferencial con
la fórmula de Heun. Use h = 0.1
y’ - y - x + x2 - 1 = 0, y(0) = 1
Solución
y’ = f(x, y) = x - x2 + y + 1, x0 = 0, y0 = 1, h = 0.1
Cálculos
i=0: K1 = hf(x0, y0) = 0.1 f(0, 1) = 0.1 (0 - 02 +1 + 1) = 0.2000;
K2 = hf(x0 + h, y0 + K1) = 0.1 f(0.1, 1.2) = 0.1 [ 0.1 – 0.12 + 1.2 +1] = 0.2290
y1 = y0 +1/2 (K1 + K2) = 1 + 0.5(0.2000 + 0.2290) = 1.2145
x1 = x0 + h = 0 + 0.1 = 0.1

i=1: K1 = hf(x1, y1) = 0.1 f(0.1, 1.2145) = 0.1 (0.1 – 0.12 + 1.2145 + 1) =0.2305;
K2 = hf(x1 + h, y1 + K1) = 0.1 f(0.2, 1.4450) = 0.1 [ 0.2 – 0.22 + 1.4450 + 1] = 0.2605
y2 = y1 + 1/2 (K1 + K2) = 1.2145 + 0.5(0.2305 + 0.2605) = 1.4600
x2 = x1 + h = 0.1 + 0.1 = 0.2
Para comprobar comparamos con la solución exacta: y(x) = x + x2 + ex
y(0.1) = 1.2152
y(0.2) = 1.4614
El error de truncamiento en cada paso está en el orden de los milésimos, coincidiendo
aproximadamente con E=O(h3).

Ejemplo 5.

Y′=2x–3y+1, y(1)=5, y(1.2)

Primer caso h=0.1

Pasos:

Escribimos la ED en la forma: dy/dx = f(x, y), para extraer su segundo miembro

dy/dx =2x–3y+1

ii. Definimos x0, y0 y h de acuerdo a los datos del problema

x0=1,

y0=5,

Para este primer caso: h=0.1

iii. Planteamos la ecuación de Euler utilizando los datos iniciales:

𝒚𝟏 = 𝒚𝟎 + h∗(2×0–3 𝒚𝟎 + 1)

𝒚𝟏 = 5 + (0.1) (2(1)–3 (5) +1)

Desarrollamos hasta el valor buscado en xx, en este caso: x=1.2 x=1.2. (Ver datos del
problema)

y1= 5 + (0.1) (2–15+1)

=5+(0.1) (−12)

=5–1.2y1

= y (1.1) = 3.8000

𝒚𝟏+𝟏 = 𝒚𝟏 + 𝒉𝒇(𝒙𝟏 , 𝒚𝟏 )

𝒚𝟐 = 𝒚𝟏 + 𝒉 ∗ (𝟐𝒙𝟏 − 𝟑𝒚𝟏 + 𝟏)

𝒚𝟐 = 𝟑. 𝟖 + (𝟎. 𝟏) (𝟐(𝟏. 𝟏) − 𝟑(𝟑. 𝟖) + 𝟏)

𝒚𝟐 = 𝟑. 𝟖 + (𝟎. 𝟏)(𝟐. 𝟐 − 𝟏𝟏. 𝟒 + 𝟏)

= 3.8 + (0.1)(-8.2)

= 3.8 – 0.82

𝒚𝟐 = 𝒚(𝟏. 𝟐) = 𝟐. 𝟗𝟖𝟎𝟎