2Eva_IT2012_T2 Modelo de clima – Métodos numéricos

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