Programación científica (I): Representando el atractor de Lorenz con C++20, Boost.Numeric.Odeint y Dlib

Artículos de la serie:

1. Representando el atractor de Lorenz con C++20, Boost.Numeric.Odeint y Dlib
2. Análisis dimensional con mp-units, regresiones no-lineales con Dlib.Optimization y visualización de datos con Matplot++
3. Rectas de mejor ajuste por mínimos cuadrados (mp-units, Matplot++)
4. Ceros de una función
5. Comprobando la ley de Zipf
6. Instalación de Matplot++ y mp-units con Mingw-w64

Introducción

Última actualización del artículo: 9 de noviembre de 2021.

El sistema de Lorenz es un sistema dinámico determinista 3-dimensional descrito por las siguientes ecuaciones diferenciales ordinarias (EDOs) no-lineales [1]:

dx/dt = σ(y - x)

dy/dt = x(ρ - z) - y

dz/dt = xy - βz

donde σ, ρ y β son parámetros reales positivos. Dada una posición inicial p₀ = (x₀, y₀, z₀) ∈ ℝ³, el sistema seguirá una órbita en el espacio de fases parametrizada por el tiempo t ∈ [0, +∞).