PyPlane, un GUI en Python para resolver ecuaciones diferenciales ordinarias. Últimamente estoy viendo cada vez más artículos sobre lo quemados que están muchos desarrolladores en el mundo del software abierto y sobre ciertas actitudes tóxicas que no ayudan. Otro libro: “Research Software Engineering with Python – Building software that makes… Leer más »PyboEnlaces de la semana: 2021-01-11 a 2021-01-17
Hoy vamos a hablar de ecuaciones diferenciales (y de modelos epidemiológicos). Imaginad que existe una pandemia mundial de una nueva enfermedad que llamaremos X. Vamos a establecer una serie de ideas para poder trabajar en un modelo de simulación de estas pandemias. El modelo SIR Voy a empezar hablando un… Leer más »Modelo SIR, modelo epidemiológico, con Python
Introducción En este artículo os voy a presentar el proyecto FEniCS, una colección de bibliotecas escritas con interfaz en Python para la resolución de ecuaciones diferenciales por el método de los elementos finitos. FEniCS proporciona un método muy inteligente para automatizar los procesos más laboriosos de la solución de este… Leer más »FEniCS: Resolución de ecuaciones diferenciales en Python
Introducción Resulta intuitivo pensar que las poblaciones de un animal depredador y su presa están relacionadas de algún modo en el que si una aumenta, la otra lo hace también. Utilizaremos como ejemplo en este artículo un ecosistema aislado y formado por leones y cebras que viven en armonía, es… Leer más »Ecuaciones de Lotka-Volterra: modelo presa depredador
En este rápido apunte vamos a ver cómo construir una matriz tridiagonal en Python utilizando NumPy y SciPy. Una matriz tridiagonal es una matriz cuadrada que solamente tiene elementos distintos de cero en su diagonal principal y en las dos diagonales adyacentes a esta (la superdiagonal y la subdiagonal). Las matrices tridiagonales aparecen mucho en cálculo numérico, por ejemplo en la discretización de ecuaciones diferenciales, y tienen la característica de ser matrices dispersas (en lugar de densas) al ser la mayoría de sus elementos cero.
Sin que sirva de precedente, hoy vamos a escribir código que sea compatible tanto con Python 2 como con Python 3. Es un cambio nimio, pero merece la pena ir acostumbrándose a pensar que tarde o temprano habrá que abandonar Python 2 🙂
En esta entrada se ha usado python 2.7.3, numpy 1.6.1 y scipy 0.10.1 y es compatible con python 3.2.3.
Leer más »Cómo crear una matriz tridiagonal en Python con NumPy y SciPy
Hoy veremos cómo computar con Python la región de estabilidad absoluta de un método numérico para resolver problemas de valores iniciales en ecuaciones diferenciales ordinarias, una herramienta muy importante a la hora de escoger un método numérico adecuado para integrar nuestro problema concreto. Se trata simplemente de otro ejemplo aplicado de lo que publicamos hace unos días sobre cómo pintar curvas de nivel con matplotlib; si quieres ver otro más, puedes leer nuestro ejemplo de uso de Basemap y NetCDF4, donde vimos cómo representar datos climatológicos.
En esta entrada se ha usado python 2.7.3, numpy 1.6.1 y matplotlib 1.1.0.
Leer más »Regiones de estabilidad de métodos numéricos en Python
Introducción En Python tenemos numerosas herramientas listas para que podamos integrar ecuaciones diferenciales ordinarias sin tener que preocuparnos en implementar un esquema numérico. Sin ir más lejos, en el módulo integrate de SciPy existen varias funciones a tal efecto. La función odeint es una interfraz en Python a la biblioteca ODEPACK,… Leer más »Integrando ecuaciones diferenciales: método leapfrog en Python