Scribd
Cargar
300 vistas2 páginas

Resolución de EDOs en Python

Para resolver ecuaciones diferenciales de segundo orden o de orden superior en Python, primero se debe reducir la ecuación a un sistema de ecuaciones de primer orden. Luego, la librería SciPy permite resolver numéricamente el sistema de ecuaciones de primer orden mediante el método odeint, el cual requiere definir la función que describe el sistema y las condiciones iniciales.

Cargado por

ISAEL Jaimes
Derechos de autor
© © All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como DOCX, PDF, TXT o lee en línea desde Scribd
300 vistas2 páginas

Resolución de EDOs en Python

Para resolver ecuaciones diferenciales de segundo orden o de orden superior en Python, primero se debe reducir la ecuación a un sistema de ecuaciones de primer orden. Luego, la librería SciPy permite resolver numéricamente el sistema de ecuaciones de primer orden mediante el método odeint, el cual requiere definir la función que describe el sistema y las condiciones iniciales.

Cargado por

ISAEL Jaimes
Derechos de autor
© © All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como DOCX, PDF, TXT o lee en línea desde Scribd

¿Por que se cambian las ecuaciones de 2do. orden a 1er. orden en Python?

Python solo es capas de resolver ecuaciones de primer orden, por lo cual para resolver
ecuaciones de 2do. Orden o de orden superior tenemos que hacer la reducción la ecuación a
una de primer orden.

Solución de un PVI usando Python:

Para solucionar un PVI en Python se hará uso de la instrucción odeint, que es un método
que esta implementado en la librería scipy. odeint soluciona la ecuación diferencial de
forma numérica en una determinada región (Se debe garantizar solución endicha región) La
sintaxis que usaremos será la siguiente:

Solucionar la ecuación diferencial:


Para solucionar la ecuación diferencial usando “odeint” lo primero es dividir el problema
de forma que se tenga el siguiente sistema de ecuaciones.

import numpy as np
import [Link] as plt
from [Link] import odeint
def df(y, t):

y1 = y[0]
y2 = y[1]
# Se le asignan una posición en el vector solución a y1, y2

dy1 =
dy2 =

return (dy1, dy2)

y0 = (0, 0) , # Condiciones iniciales

t = [Link](0, 0, 0) , # Definición del rango

solucion = odeint(sistema, y0, t) ,


# La función odeint se encarga de dar los argumentos a la función df(el rango “x” y la
variable “y”)

y1 = solucion[:, 0] , # Toma el vector correspondiente a la solución de y1


y2 = solucion[:, 1] , # Toma el vector correspondiente a la solución de y2

[Link](t, y1)
[Link](t, y2)
[Link](True)

También podría gustarte