2Eva_IT2012_T2 Modelo de clima

2da Evaluación I Término 2012-2013. 28/Agosto/2012. ICM00158

Tema 2. (20 puntos) El meteorólogo Edward Lorenz propuso inicialmente el siguiente sistema para predecir el comportamiento del clima:

\begin{cases} x'(t) = \alpha (y(t) - x(t)) \\ y'(t) = \rho x(t) - y(t) - x(t) z(t) \\ z'(t) = -\beta z(t) + x(t) y(t) \end{cases}

Para su estudio eligió los parámetros α = 10, β = 8/3 , ρ=28 con las condiciones iniciales x(0) = 10, y(0) = 7, z(0) = 7

Use el método de Runge-Kutta clásico de cuarto orden con h = 0.25 para calcular la solución cuando t=1


Referencia: Chapra 28.2 p833 pdf857

 

https://en.wikipedia.org/wiki/Lorenz_system