SIMULACIÓN DE BIOSISTEMAS, ABRIL 2020 1

Metodo de Runge-Kutta de orden 4, Heun e

Interpolacion polinomica de Lagrange
Facultad de Ingenierı́a, Ingenierı́a Biomédica
Universidad Autonóma de Manizales
Simulación de Biosistemas
Esteban Mercado Ruiz, Jacobo Quintero Alzate, Angie Rosseth Valencia.

Resumen—La solución de ecuaciones diferenciales se puede cartesianas xy (o simplemente plano xy), donde se muestra
hacer por medio de diferentes métodos. Los métodos numéricos el comportamiento de la pendiente (derivada) que le corres-
como el campo de direcciones, la interpolacion polinomica de ponde a la curva solución. En matemáticas, en el ámbito de
Lagrange, el metodo de Heun y el Metodo de Runge-Kutta de
orden 4 sirven para solucionar ecuaciones diferenciales. Para ecuaciones diferenciales ordinarias (EDO), se emplea este pro-
poder probarlo, se utiliza una ecuación diferencial propuesta: f= cedimiento para evaluar el comportamiento de las soluciones
x*y - x2 , a la cual se le debe hallar la solución mediante los (que corresponden a funciones). Este elemento es parte de lo
métodos mencionados anteriormente. El objetivo es implementar que se conoce como teorı́a cualitativa para EDO[1].
los códigos de los métodos en Matlab Para obtener las gráficas
que describen y se aproximan a la solución a la ecuación
diferencial. Se evidencia que el método de Runge kutta resulta
ser más efectivo al acercarse a la solución real y normal de la II-B. Interpolacion polinomica de Lagrange
Ecuacion diferencial.
La interpolacion es una tecnica matematica que permite
I. INTRODUCCI ÓN describir un conjunto de puntos mediante alguna funcion,
adicionalemnete tien utilidad cuando se desea aproximar una
L os métodos numéricos se presentan como una opción
factible para la soluciones de problemas que tienen implica-
funcion por otra funcion matematica mas simple.

El polinomio de interpolación de Lagrange, simplemente

ción o una relación con las ecuaciones diferenciales, muchos
es una reformulación del polinomio de Newton que evita los
de los modelos que rigen diferenres problemas estan dados
cálculos de las diferencias divididas.
por ecuaciones diferenciales que describen estos de manera
que se toman en cuenta factores y variables que la modelan, Demostracion: La funcion que se busca es una funcion
para esto se necesitan metodos que permitan darle forma a polinomica L(x) de grado k con el problema de interpolacion
estos, dentro de ellos se pueden encontrar formas como el que puede tener tan solo una solucion, pues la diferencia entre
campo de direcciones, Interpolacion polinomica de Lagrange, dos tales soluciones, serıa otro polinomio de grado k a lo
metodo de Heun y Runge-Kutta 4, con los cuales se ve sumo, con k+1 ceros. Por lo tanto, L(x) es el unico polinomio
una alternativa diferente al los vistos anteriormente, de esta interpolador[2]
manera se percibe una gran diferencia respecto al porcentaje
de error y la exactitud de los valores encontrados, pues los
metodos nombrados entregan una cantidad de puntos que
aproximadamente describen la pendiente de la curva próxima
a la función por la cual se rige el sistema, siempre teniendo
en cuenta los valores iniciales y las condiciones que dan un
punto de partida de este. El objetivo del presente documento
es poder realizar un programa que permita usar los metodos ya Figura 1. Formula general del polinomio de Lagrange
nombrados para darle una solución aproximada o semejante a
la ecuacion diferencial presentada.

II. M ARCO TE ÓRICO II-C. Metodo de Heun

II-A. Campos de direcciones El metodo de Heun es un procedimiento numérico para

resolver ecuaciones diferenciales ordinarias (EDO) con un
El campo de direcciones, es un bosquejo con pequeños valor inicial dado. Es un método para mejorar la estimación
segmentos de recta trazados en un sistema de coordenadas de la pendiente que involucra la determinación y promediado
de dos derivadas para el intervalo (una en el punto inicial y
otra en el punto final), la pendiente calculada en la estimación
previa, no es para la respuesta final, sino para una predicción
SHELL et al.: BARE DEMO OF IEEETRAN.CLS FOR JOURNALS SIMULACIÓN DE BIOSISTEMAS, ABRIL 2020

método Runge-Kutta)). Se puede ver como resultado de la

aplicacion de la regla de Simpson como formula de cuadratura
numerica[?]. Para determinar cada elemento de la sucesion
xk se realizan cuatro estimaciones previas, de acuerdo con el
siguiente esquema:

Figura 5. Metodo de Runge-Kutta

Figura 2. Ejemplo de grafica representativa de la interpolacion polinomica

de Lagrange
Ası́, el siguiente valor (yn+1) es determinado por el presente
valor (yn) más el producto del tamaño del intervalo (h) por una
pendiente estimada. La pendiente es un promedio ponderado
intermedia. Esta ecuación es llamada predictor. Mejora una de pendientes, donde k1 es la pendiente al principio del
estimación de yi+1 que permite el cálculo de una estimación intervalo, k2 es la pendiente en el punto medio del intervalo,
de la pendiente al final del intervalo[5]. usando k1 para determinar el valor de y en el punto (xn+(h/2))
usando el método de Euler. K3 es otra vez la pendiente del
punto medio, pero ahora usando k2 para determinar el valor
de y; k4 es la pendiente al final del intervalo, con el valor de y
determinado por k3. Esta forma del método de Runge-Kutta,
es un método de cuarto orden lo cual significa que el error
Figura 3. Ecuacion predictora del metodo de Heun por paso es del orden de O(h5 ) mientras que el error total
acumulado tiene el orden O(h4 ).

Y con lo anterior de la figura 5. se define:

Figura 6. Metodo de Runge-Kutta ecuacion

III. M ETODOLOG ÍA

Para la implementacion del metodo de Runge-Kutta de

orden cuatro en conjunto con el campo de direcciones, la
Figura 4. Grafica representativa del metodo de Heun interpolacion polinomica de Lagrange y el metodo de Heun,
se realizo inicialmente una declaracion de variable, vectores y
elementos a utilizar para el correcto desarollo, empezando asi
II-D. Metodo de Runge-Kutta por nuestra ecuacion diferencial que la alberga f que puede
ser descrita como: f= x*y - x2 .
Los métodos de Runge-Kutta son un conjunto de métodos
genéricos iterativos, explı́citos e implı́citos, de resolución Pasando asi por (xi) y (xf) que corresponden al punto
numérica de ecuaciones diferenciales. Se define como el inicial que es igual a -1.5 y punto final que es igual a 1.5
Método de iteraciones que extiende la idea geométrica y utiliza respectivamente, es decir, desde que valor de partida quiero la
varias derivadas o tangentes intermedias para poder aproximar funcion y en que valor quiero que me finalice, siguiendo en
la función desconocida[4]. un orden consecuente la variable (n) me es igual al numero
de pasos que tiene un valor de 100, la variable (h) que
Un miembro de la familia de los métodos Runge-Kutta me representa el numero de tramos que hay desde el punto
usado ampliamente es el de cuarto orden. Es usado tanto inicial hasta el punto final divido en el numero de pasos, es
que a menudo es referenciado como ((RK4)) o como ((el decir repartiendo y organizando el vector en 100 posiciones;

2
SHELL et al.: BARE DEMO OF IEEETRAN.CLS FOR JOURNALS SIMULACIÓN DE BIOSISTEMAS, ABRIL 2020

pasando a el vector (x) que evalua desde el inicio hasta el final

incrementandose en un paso igual a h.

Posteriormente nos encontramos con (y0) que es la condi-

cion inicial con un valor igual a cero de la ecuacion diferencial,
Luego esta (y) que es una matriz que procede como la solucion
a la ecuacion diferencial y que en su primera posicion tendra
como valor la condicion inicial mencionada anteriormente
denotandose de la siguiente manera: y(1)= y0. Transcurriendo
asi hasta el ciclo for, en donde vemos la aplicacion del metodo
de Runge-Kutta de orden cuatro, que es el aumento de paso
por la ecuacion diferencial evaluada en un punto de (i,i)
demostrando asi el avance de (k) con el orden dado por el
metodo aplicado k1, k2, k3 y k4, calculando como por ultima
instancia del ciclo (y) que como se menciono es la solucion
a la ecuacion diferencial que es y(i) en una posicion (i), mas Figura 8. Grafica de Runge-Kutta obtenida aumentada
la suma de las pendientes como corresponde.
V. C ONCLUSI ÓN

los métodos numéricos son útiles para la resolución de

IV. R ESULTADOS
ecuaciones diferenciales. La implementación de métodos co-
mo la interpolacion polinomica de Lagrange, el metodo de
Al finalizar el codigo del metodo de Runge-Kutta en con- Heun y el metodo de Runge-Kutta 4, ofrecen soluciones muy
junto con el campo de direcciones, la interpolacion polinomica aproximadas a la solución real y normal de la ecuación dife-
de Lagrange y el metodo de Heun, se obtuvo un grafica rencial, siendo entre estos unos más efectivos que otros pero
que va desde el punto inicial (xi) hasta el punto final (xf) igualmente útiles como es el caso del método RK 4, el cual
que corresponde a -1.5 y 1.5 respectivamente, la grafica es más utilizado para resolver numéricamente problemas de
que se expresa en color amarillo corresponde al metodo de ecuaciones diferenciales ordinarias, con condiciones iniciales
Runge-Kutta, la de color azul a la interpolacion polinomica e y que proporciona un pequeño margen de error con respecto
Lagrange, la de color marron al metodo de Heun y por ultimo a la solución real del problema y es fácilmente programable
en el fonde de color azul se encuentra el campo de direcciones. en un software para realizar todas las iteraciones necesarias.

R EFERENCIAS

[1] Campo de direcciones - Wikipedia, la enciclopedia libre. Retrieved April

02, 2020, from https://es.wikipedia.org/wiki/Campodedirecciones
[2] Interpolación de Lagrange - Teorı́a - Interpolación. Retrieved April 17,
2020, from http://interpolacion.wikidot.com/lag-teoria
[3] Matemática Departamento de Aplicada. Grado de ingerniria. Universidad
de Granada. from http:// http://www.ugr.es/ arobles/Telecos/Tema11.pdf
[4] Rodriguez Barajas, M. A. (2016). Euler y runge kutta p. Retrieved from
https://es.slideshare.net/MarcoAntonioRodrguez11/euler-y-runge-kutta
[5] Weisstein, E. W., Weisstein, E. W. Runge-Kutta Method. MathWorld.
Retrieved from http://mathworld.wolfram.com/Runge-KuttaMethod.html

Figura 7. Grafica de Runge-Kutta obtenida

Para expresar una comparacion entre los metodos empleados

para hallar la solucion de la ecuacion diferencial, se puede
decir que los metodos utilizados hacen una aproximacion a la
solucion de la ED, en este caso en particular el mas usado a
estas aplicaciones es el RK4”por lo que se puede decir que
es el metodo que mas se asemeja a hallar la solucion, ya que
estos utilizados poseen un error que para poder comparar se
debe utilizar la solucion real de la ecuacion diferencial, con
este metodo se obtienen unas mejores aproximaciones a la
solucion exacta.