UNIVERSIDAD NACIONAL

AGRARIA LA MOLINA

MÉTODO DE EULER
Ecuaciones diferenciales ordinarias

PROFESOR: Edwin Baldeón.

INTEGRANTES:
 Ñahui Gutiérrez, Alejandra Carolina
 Vega León, Esther Bertha
 Velásquez Quispe, Laura Stefany
 Villalobos Picón, Fabricio Miguel Ángel
 Vilcapoma Torres, Rodrigo Amador
 Zúñiga Pariona, Alejandra María
GRUPO: 10

I. INTRODUCCION

2018 - I
El método de Euler es un procedimiento de integración numérica que se utiliza para
poder resolver ecuaciones diferenciales ordinarias (EDO) a partir de un valor inicial
dado, es decir, con el método de Euler se logra obtener una solución aproximada en un
conjunto finito de puntos partiendo de la ecuación de la recta, pero con el método de
Euler Mejorado es más exacto al momento de converger, y de las formulas establecidas
por dichos métodos teniendo en cuenta que su solución general es compleja[ CITATION
Ste07 \l 10250 ].

Es importante hacer mención, de lo útil que resulta esta aplicación a la hora de moldear
situaciones físicas en la ingeniería donde se encuentran funciones desconocidas con
respecto a una o varias variables independientes[ CITATION Ste07 \l 10250 ].

Los objetivos son los siguientes:

1. Conocer y desarrollar el Método de Euler tanto su aplicación en los problemas

matemáticos y el desarrollo en la programación (en este caso en el lenguaje C++
y en su estructura en código para razonar su procedimiento).
2. Observar y analizar las ventajas y desventajas que tiene el Método de Euler para
encontrar las aproximaciones de una función con un porcentaje de error
aceptable.
3. Obtener una visión clara de cuál es el comportamiento de este Método tanto en
las matemáticas como en la programación.

II. REVISION LITERARIA

II.1 Ecuaciones diferenciales ordinarias

La gran mayoría de los fenómenos físicos y químicos aplicados a la ingeniería, pueden

ser modelados mediante ecuaciones diferenciales ordinarias, por lo cual se hace
necesario estudiar las técnicas más apropiadas que nos permiten resolver estos modelos
que por lo general no contienen soluciones analíticas. Dentro de estas técnicas de
resolución destacan aquellas de solución inmediata y otras que involucran ecuaciones
predictores y concretas las cuales van a ser analizadas aquí.

II.1.1 Concepto

Es aquella ecuación que presenta una sola variable independiente, por lo tanto sus

dy
derivadas son totales. =f ´ ( x ; y)
dx

X: variable independiente

Y: variables dependientes

Donde X puede estar en forma explícita o implícita:

dy dy
= y+ x = y+ seny
dx dx

Explicita, en variable independiente Implícita en variable dependiente

II.2 Ecuación diferencial parcial

Es aquella en la que existen 2 o más variables independientes, y por ello sus derivadas
son parciales.

Por ejemplo, la ecuación de transferencia de calor en estado no estacionario en una

dimensión:

dT d ² T
= T= f (t; x) t y x: variables independientes.
dt dx ²

II.3 Orden de Ecuación diferencial:

El orden de una ecuación diferencial es la derivada de mayor orden que aparece en una
ecuación
 2
yᴹ + x yⁿ + ( 1−x ) y ´ + y=0 3° orden

II.4 Grado de Ecuación diferencial

Es el grado algebraico de la derivada de mayor orden que s encuentra en la ecuación

yᴹ + x 2 ( yⁿ ) √ y ´ + y=0 Es de 1° grado, ya que el grado de la derivada mayor es 1.

II.5 Métodos para resolver una ecuación diferencial ordinaria:

II.5.1 Sin computadora:

Las ecuaciones diferenciales ordinarias podrían resolverse usando técnicas de

integración analíticas:

Por ejemplo

Integración indefinida debido a que no se

especifica límites de integración

Una solución analítica a la integral indefinida es evaluándola de forma exacta como una
ecuación.

No están disponibles las soluciones exactas para muchas ecuaciones diferenciales

ordinarias de importancia práctica. Los métodos numéricos ofrecen la única alternativa
viable para estos casos. Como estos métodos numéricos requieren de computadoras, los
ingenieros antes de las computadoras se veían muy limitados en el alcance de sus
investigaciones.

El método que utilizaron y el más importante fue la linealizacion, en el cual consistía

convertir una ecuación diferencial ordinaria no lineal, en una ecuación diferencial
ordinaria lineal. La importancia de una ecuación diferencial ordinaria lineal es que se
resuelven analíticamente, en cambio la mayoría de las ecuaciones no lineales no pueden
resolverse de manera exacta. Así una táctica para resolver ecuaciones no lineales era
linealizarlas.

Un ejemplo:
 ᶿ= Angulo de desplazamiento

g =constante gravitatorio

l =longitud de péndulo

Predecir el momento de un péndulo oscilante usando la segunda ley de newton para

obtener una ecuación diferencial.

La ecuación diferencial no es lineal debido a la función de sen.

Una forma de darle solución es que para pequeños desplazamientos del Angulo, se

considere la siguiente deducción: Senᶿ=ᶿ

Y reemplazando, queda una ecuación diferencial ordinaria lineal:

El cual se puede resolver en forma analítica.

La linealizacion es una buena herramienta pero no se puede usar en todos los casos, por
ejemplo en grandes desplazamientos del Angulo.

Las computadoras usando métodos numéricos, coloca esta opción al alcance de todos
los ingenieros.

II.5.2 Con computadora:

Luego de una aplicación para la solución de ecuaciones diferenciales ordinarias no

lineales y lineales a través de métodos numéricos, se establecieron distintos métodos
para su resolución.

Ejem: El programa C++

Resultados del programa en C++

Figura 1. Resultados del programa C++

Al principio del ejecutable del código en C++ del método de Euler nos piden ingresar
los valores correspondientes para realizar el programa y la solución aproximada para la
solución de ecuaciones diferenciales. El condigo en si realiza con ciclos por
las operaciones y cálculos indicados por el usuario donde este mismo introdujo las
formulas del método y la organización para que se fuera desarrollando. También en otro
ciclo por se colocaron las fórmulas para los valores verdaderos de la ecuación
diferencial y los errores que se fueron calculando entre el método de Euler y el del valor
verdadero. La pantalla muestra cómo se imprime en forma de lista las operaciones
realizadas con el desarrollo del método en el programa.

ANÁLISIS DEL PROGRAMA EN C++ DEL “MÉTODO DE EULER”

El desarrollo del programa en c++ se basó en la creación de dos ciclos por para ir
desarrollando las formulas y ecuaciones del problema a partir de los ciclos con ciertas
condiciones de inicio que fueron ingresadas por el usuario. El primer por desarrolla
la parte de las formulas las cuales sustituyen los valores iniciales ingresados por el
usuario y después de sacar el limite al cual debe de detenerse el ciclo como también
para el segundo ciclo por. Ya resueltos la sustitución y desarrollo de las fórmulas para la
aproximación de los resultados el segundo ciclo por realiza las operaciones
correspondientes a las funciones verdaderas y los errores para también imprimirlos en
ese mismo por y así mostrarlos en el ejecutable
II.5.2.1 Método de Euler:

Es uno de los métodos más sencillos de implementar para resolver ecuaciones

diferenciales ordinarias, el cual consiste en hacer una integración de una ecuación
diferencial.

En este método se resuelve la ecuación diferencial de manera numérica. Es decir que

usamos la ecuación diferencial como base de un algoritmo para aproximar la solución
desconocida. (Chapra, 2007)

Se llama método de Euler al método numérico consistente en ir incrementando paso a

paso la variable independiente y hallando la siguiente imagen con la derivada.
(Seminario, 2004)

El método de Euler se presenta mediante le problema modelo. El primer paso es dividir

el intervalo [a,b] en m subintervalos de longitud h = b-a/m mediante puntos xi, con i=o,
…...,m. (Huerta, 1998)

Figura 2. Modelo de grafica de Euler

Debe verificarse para todo los puntos [a,b]. Si se particulariza para un punto xi, de la
discretizacion resulta:

dy
( xi )=f ( xi , yi)
dx
Donde se emplea la notación yi = y(xi). La idea básica del método de Euler es
aproximar la derivada en el punto xi, mediante un cociente incremental hacia adelante,
utilizando los puntos xi y xi+1.

dy Yi+1−Yi yi+1− yi
( xi )= =
dx xi+1−xi h

Esta ecuación puede interpretarse geométricamente como aproximar la recta tangente a

la función y en xi por la recta secante que pasa por los puntos (xi,yi) y (xi+1,yi+1).

Fundamento de la ecuación de Euler

Figura 3. Fundamento del metido de Euler

Desarrollo:

Sea el problema de valor inicial

Ec. 1

Supongamos que tiene una solución única φ(x) en un algún intervalo con centro en x o.
Sea h > 0 y consideremos puntos igualmente espaciados
Los valores de la solución φ (xn) se pueden aproximar con y n, donde los valores de yn se
obtienen como sigue.

En el punto (x0,y0) la pendiente de la solución de es dy/ dx = f(x 0,y0). Por lo tanto, la

recta tangente a la curva solución en el punto (x0,y0) es:

Ec. 2

Si se usa Ec. 2, como una aproximación a φ(x), en el punto x1 = x0 + h

Ec. 3

En seguida, empezando en el punto (x1,y1) con pendiente f(x1,y1), se tiene la recta y = y 1

+ (x – x1)f(x1,y1) al pasar de x1 a x2 = x1 + h nos da la aproximación:

Y al repetir el procedimiento se obtiene:

Así, sucesivamente hasta obtener un número “n” constante. Este sencillo procedimiento
se llama Método de Euler y se resume mediante las siguientes fórmulas sucesivas.

x (n+1 )=x n+ h

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

x (n +1)−x n
h=
n
II.5.2.2 Método de Taylor:

Una manera para reducir el error con el método de Euler, seria incluir términos de orden
superior en la expansión de la serie de Taylor para la solución.

II.5.2.3 Método de Rounge Kutta:

Logran la exactitud del procedimiento de la serie de Taylor sin necesitar el cálculo de

derivadas de orden superior.

II.5.2.4 Métodos de mejora de Euler

Un motivo fundamental de error para Euler era suponer que la derivada al inicio del
intervalo es la misma durante todo el intervalo.

 Método de Heun
 Método de punto medio

ANALISIS DEL ERROR PARA EL METODO DE EULER

La solución numérica de las EDO implica dos tipos de error

Errores de truncamiento o de discretizacion

Son originados por la naturaleza de las técnicas empleadas para aproximar los valores
de la función de y”

Errores de redondeo

Causados por el número limitado de cifras significativas que una computadora puede
retener.

Los errores de truncamiento se componen de dos partes: la primera, es un error de

truncamiento local que resulta de una aplicación del método considerado, en un solo
paso. La segunda es un error de truncamiento propagado que resulta de las
aproximaciones producidas durante los pasos previos. La suma de los dos es el error de
truncamiento global o total.

Al adquirir cierta comprensión de la magnitud y de las propiedades del error de

truncamiento puede desarrollarse el método de Euler directamente de la expansión de la

Y”= f (x, y)
serie de Taylor. Para ello, observe que la ecuación diferencial que se va a integrar será
de la forma general:

Donde y” es igual a dy/dx, x y son las variables independiente y dependiente,

respectivamente. Si la solución (es decir, la función que describe el comportamiento de
y) tiene derivadas continuas, se representa por una expansión de la serie de Taylor
respecto a un valor inicial. (Chapra, 2011)

Figura4. Error de truncamiento

La serie de Taylor permite solo una estimación de error de truncamiento local; es decir,
el error generado durante un solo pasó del método. No ofrece una medida de error
propagado, por lo tanto, ni del error de truncamiento global.

Aunque estas limitaciones impiden el análisis exacto del error, en la mayoría de los
problemas prácticos, la serie de Taylor brinda una valiosa ayuda en la comprensión en el
comportamiento del método de Euler.

Según Carnahan, las observaciones permiten establecer las siguientes conclusiones

útiles:
a) se puede reducir el error disminuyendo el tamaño del paso

b) el método dará como resultado predicciones sin error si la función que se analiza (es
decir, la solución de la ecuación diferencial) es lineal, debido a que en una línea recta la
segunda derivada es 0.

Esta última conclusión tiene un sentido intuitivo, puesto que el método del Euler usa
segmentos de línea recta para aproximar la solución de ahí el método de Euler se le
conoce como un método de primer orden. (Carnahan y sus colaboradores, 1969)

ALGORITMO PARA EL METODO DE EULER

Los algoritmos para las técnicas como de un paso como el método de Euler son muy
simples de programar. Como se especificó al inicio de este capítulo, todos los métodos
de un paso tienen la formula general

Nuevo valor= valor anterior + pendiente* tamaño de paso

Es lo único que difieren los métodos es en el cálculo de la pendiente. Suponga que usted
quiere realizar el cálculo simple. Es decir, a usted le gustaría utilizar el método de Euler
para integrar.

Con la condición inicial de qué y= 1, en x=0.

Usted quiere integrarla hasta X=4 usando un tamaño de aso de 0.5, y desplegar todos los
resultados. (Davis, 1990)

MEJORAS DEL MÉTODO DE EULER

Un motivo fundamental de error en el método de Euler, es suponer que la derivada al

inicio del intervalo es la misma durante todo el intervalo. Has dos modificaciones
simples para evitar esta consideración.
METODO DE HEUN

Es un método para mejorar la estimación de la pendiente, emplea la determinación de

dos derivadas en el intervalo (una en el punto inicial y otro en el final). Las dos
derivadas s e promedian después con la finalidad de obtener una mejor estimación de la
pendiente en todo el intervalo. Este procedimiento conocido como el método de He un,
se presenta en forma gráfica. (Riggs, 1994)

Figura 5. A) Predictor B) corrector

Ejemplo 1:

Con el método de Euler integre numéricamente la ecuación

dy
=−2 x 3+ 12 x 2−20 x+ 8.5
dx

Desde x=0 hasta x=4 con un tamaño de paso 0,5. La condición inicial en x=0 es y=1.

1. Planteamiento del problema:

Nuevo valor= Valor anterior

+ Pendiente x Tamaño de
paso
Recordar el método general

En forma matemática y i+1 = yi + ø xh

ec.1
x
y i+1 = y i + f (¿ ¿ i, y i ) xh
¿
La solución exacta está dada por
4 3 2
ec.2 y=−0.5 x +4 x −10 x +8.5 x+ 1

2. Solución del problema:

Sabiendo que X= [0,4] y h=0, 5

Hallamos el valor de todas las x en ese intervalo:

x i=x o +h

x o=0

x 1=0+0,5=0,5

x 2=0,5+0,5=1

x 3=1+0,5=1,5

x 4=1,5+ 0,5=2

x 5=2+0,5=2,5

x 6=2,5+0,5=3

x 7=3+0,5=3,5

x 8=3,5+0,5=4

Luego, estos valores ingresaran a la ec.1

x
y i+1 = y i + f (¿ ¿ i, y i ) xh
¿
Sabiendo que y i = 1 y el valor de la pendiente es igual a:

x
 f (¿ ¿ i, y i ) = −2(0)3+12( 0)2 −20 ( 0 ) +8,5=8,5
¿
y1 = 1 + 8,5 x 0,5 = 5,25

Ahora, sabiendo que y 1 = 5,25 y el valor de la pendiente es igual a:

x
 f (¿ ¿ 1 , y 1 ) = −2( 0.5)3 +12(0.5)2−20 ( 0.5 ) +8,5=1,25
¿
y 2 = 5,25 + 1,25 x 0,5 = 5,875

x
 f (¿ ¿ 2 , y 2 ) = −2(1)3 +12(1)2−20 ( 1 )+8,5=−1,5
¿
y 3 = 5,875 + (-1,5) x 0,5 = 5,125

x
 f (¿ ¿ 3 , y3 ) = −2( 1,5)3 +12(1,5)2−20 ( 1,5 ) +8,5=−1,25
¿
y 4 = 5,125 + (-1,25) x 0,5 = 4,5

x
 f (¿ 4 , y 4)
¿ = −2( 2)3 +12(2)2−20 ( 2 ) +8,5=0,5
¿
y 5 = 4,5 + 0,5 x 0,5 = 4,75

x
 f (¿ ¿ 5 , y5 ) = −2( 2,5)3 +12(2,5)2−20 ( 2,5 ) + 8,5=2,25
¿
y 6 = 4,75 + 2,25 x 0,5 = 5,875

x
 f (¿ ¿ 6 , y 6) = −2( 3)3 +12(3)2−20 ( 3 )+ 8,5=2,5
¿
y 7 = 5,875 + 2,5 x 0,5 = 7,125

x
 f (¿ ¿ 7 , y 7) = −2( 3.5)3 +12(3.5)2−20 ( 3.5 ) +8,5=−0,25
¿
y 8 = 7,125 + (-0,25) x 0,5 = 7

Para hallar el y verdadero , se utiliza la ec. 2

 y=−0.5 x 4 +4 x3 −10 x 2 +8.5 x+ 1

4 3 2
y i=−0.5 ( 0 ) + 4 ( 0 ) −10 ( 0 ) +8.5 ( 0 ) +1=1

y 1=−0.5 ( 0,5 )4 + 4 ( 0,5 )3 −10 ( 0,5 )2+ 8.5 ( 0,5 ) +1=3,21875

4 3 2
y 2=−0.5(1) + 4(1) −10 ( 1 ) +8.5 ( 1 )+ 1=3

y 3=−0.5 ( 1,5 )4 + 4 ( 1,5 )3−10 ( 1,5 )2 +8.5 ( 1,5 ) +1=2,21875

4 3 2
y 4 =−0.5 ( 2 ) + 4 ( 2 ) −10 ( 2 ) +8.5 ( 2 ) +1=2

y 5=−0.5 ( 2,5 )4 + 4 ( 2,5 )3−10 ( 2,5 )2+8.5 ( 2,5 )+1=2,71875

y 6=−0.5 ( 3 )4 + 4 ( 3 )3 −10 (3 )2+ 8.5 (3 )+ 1=4

y 7=−0.5 ( 3,5 )4 + 4 ( 3,5 )3−10 ( 3,5 )2+8.5 ( 3,5 )+1=4,71875

y 8=−0.5( 4)4 + 4 (4)3−10 ( 4 )2 +8.5 ( 4 ) +1=3

Comparación de los valores verdadero y aproximado de la integral

y ' =−2 x 3 +12 x 2−20 x +8.5 , con la condición inicial de que y=1 en x=0.

x y verdadero y Euler Error relativo porcentual

Global Local
0,0 1 1
0,5 3,21875 5,25 -63,1 -63,1
1 3 5.875 -95,8 -28
1,5 2,21875 5,125 131,0 -1,41
2 2 4,5 -125.0 20,5
2,5 2,71875 4,75 -74.7 17,3
3 4 5,875 46,9 4,0
3,5 4,71875 7,125 -51,0 -11,3
4 3 7 -133.3 -53,0
 Comparación de la solución verdadera con una solucion numérica usando el método de Euler
8

4
y

0
0 0.5 1 1.5 2 2.5 3 3.5 4 4.5

Figura 6. Dependencia de tamaño de paso en la exactitud

Aplicación en la ingeniería

Ejemplo 2:

Si se drena el agua desde un tanque cilíndrico vertical por medio de abrir una válvula en
la base, el líquido fluirá rápido cuando el tanque esté lleno y despacio conforme se
drene. Como se ve, la tasa a la que el nivel del agua disminuye es:

dy
=−k √ y
dt

Donde k es una constante que depende de la forma del agujero y del área de la sección
transversal del tanque y agujero de drenaje. La profundidad del agua y se mide en
metros y el tiempo t en minutos. Si k = 0,06, determine cuándo tiempo se requiere para
vaciar el tanque si el nivel del fluido se encuentra en un inicio a 3 m. Resuelva con la
aplicación de la ecuación de Euler y escriba un programa de computadora en Excel.
Utilice un paso de 5 minutos.
 k 0.06

tamaño de paso (h) 5 min

dy
Y,m dt
=−k √ y

Y(0) = 3m
Y(35) = ?

Datos

X=Xi + n*h
n = (X-Xi)/h
n= (35-0)/5
n= 7

a) Resolviendo con el Método de Euler:

Xi+1 = Xi + h Yi+1 = Yi + pendiente*h

X0 = -5 + 5 = 0 Yi+1 = Yi +(−�√Yi) *h

X1 = 0 + 5 = 5
Y1 = 3 + (−0.06√3)*5= 2.48
X2 = 5 + 5 =10
Y2 = 2.48 + (−0.06√3)*5= 2.008
X3 = 10 + 5 =15
Y3 = 2.008 + (−0.06√3)*5= 1.583
X4 = 15 + 5 = 20
Y4 = 1.583 + (−0.06√3)*5= 1.206
X5 = 20 + 5 = 25
Y5 = 1.206 + (−0.06√3)*5= 0.877
X6 = 25 + 5 = 30
Y6 = 0.877 + (−0.06√3)*5= 0.596
X7 = 30 + 5 =35
Y7 = 0.596 + (−0.06√3)*5= 0.364

Intervalos (i) Xi Yi Yi+1

0 0 3 2.48

1 5 2.48 2.008
 2 10 2.008 1.583

3 15 1.583 1.206

4 20 1.206 0.877

5 25 0.877 0.596

6 30 0.596 0.364

7 35 0.364 0.183

b) Tablas y gráficos
3.5

2.5

1.5

0.5

0
0 5 10 15 20 25 30 35 40

Gráfico de f(xi; yi)

Respuesta: El nivel del tanque

III. CONCLUSIONES

En lo general sobre el Método pues es muy favorable su procedimiento ya que

cuenta con menos cálculos a comparación de otros métodos, pero una desventaja es
que su margen de error que es muy grande dependiendo al valor de intervalo que se
encuentre.

En conclusión es un método muy sencillo y recomendable para la práctica pero para

la exactitud estaría mejor La ventaja de este método seria que su desarrollo es muy
fácil y sencillo de realizar, pero su desventaja sería el de hacer muchas operaciones
si queremos llegar a la exactitud.

IV. BIBLIOGRAFIA
 CHAPRA, S. 2007. Métodos numéricos para ingenieros. México. Quinta
Edición.
 DAVIS M.E; 1990. Métodos numéricos para ingenieros químicos. Compañía
Editorial Continental de C.V. México. Cap. 1
 HUERTA, 1998. Métodos Numéricos: introducción, aplicaciones y
programación. Ediciones de la Universidad Politécnica de Catalunya. Catalunya,
España. Cap 5
 RIGGS J.B; 1994. An Introduction to Numerical Methodos for Chemical
Engeneers. Texas Tech University Press, Lubbock. Texas. Cap 4

 SEMINARIO, R. 2004. Métodos Numéricos para Ingeniería. En Línea.