Implement euler integration scheme and implement

Question

In 1963 Edward Lorenz derived a simple set of equations describing convection in the atmosphere:

dx/dt = σ(y-x),

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

dz/dt = xy - βz

Even though these equations are simple and deterministic, long-term behavior of solutions for some particular values of parameters (e.g. σ = 10, ρ = 28, β = 8/3) could be highly unpredictable. Small variations of initial conditions could result in drastic difference of the corresponding solutions of the system. The latter phenomenon is known as the butterfly effect: small perturbations of the atmosphere caused by the butterfly wings at one location on Earth can result (according to the model) in substantial changes in the atmosphere at another location.

Consider (1) and let the initial conditions and the interval of interest be as follows

x(0) = 0, y(0) =1, z(0) = 0, t  ∈ [0,50].

Problem 1:

Please implement Euler integration scheme for (1) with integration step Δ = 0.0001.

Problem 2:

Please implement improved Euler integration scheme for (1) with Δ = 0.0001.

Problem 3:

Please implement Runge-Kutta integration scheme for (1) with Δ = 0.0001. Change x(0) = 0 to x(0) = 0.0001.

Observe the butterflieffect.

Problem 4:

Please quantify how different are these numerical solutions from each other (take Runge-Kutta scheme with x(0) = 0 as a reference and plot their differences). Investigate what happens with these differences when the value of Δ is made 5 times smaller/larger.

Answer 1

By implementing the Euler integration scheme solve the following question showing the whole calculation part. In the year 1963 Edward Lorenz derived a simple set of equations explaining convection in the atmosphere: dx/dt = s(y-x), dy/dt = x(?-z)-y dz/dt = xy - ßz Even though such equations are simple and deterministic, long-term behavior of solutions for some specific values of parameters (example: s = 10, ? = 28, ß = 8/3) could be highly unpredictable. Little variations of initial conditions could outcome in drastic difference of the corresponding solutions of the system. The latter phenomenon is termed as the butterfly effect: small perturbations of the atmosphere caused due to butterfly wings at one place on Earth can result in substantial changes in the atmosphere at the other location. Let the initial conditions and the interval of interest be as : x(0) = 0, y(0) =1, z(0) = 0, t ? [0,50]. Please implement the Euler integration scheme for by integration step ? = 0.0001