Explore BrainMass

# Linearizing Lorenz Equations using the Implicit Euler Method

Not what you're looking for? Search our solutions OR ask your own Custom question.

This content was COPIED from BrainMass.com - View the original, and get the already-completed solution here!

I need help to linearize the Lorenz equations so that I can use Matlab to create the butterfly effect, etc. We were given the linearized equations but a couple of students pointed out that one of them was wrong. I don't know which equation is wrong, so if someone could show me how to at least linearize the first Lorenz equation using implicit Eulers, then I can reproduce the method to check the other two equations. I would like to see step by step how the method is employed. We are to only use implicit Eulers. I need them linearized so that I can use Gauss-Seidel iteration in Matlab to create the butterfly effect (we are not to use ode45 or other solvers, hence the need for GS iteration and linearized equations).

here is what the instructor gave us:

dx
dt = s(y-x)

dy
dt = x(r-z)-y

dz
dt = -bz + xy

s = 10, r = 28, b = 8/3

s = sigma, r = rho, b = beta
dt = 10^-4

initial conditions:

x = 2, y = 2, z = 2

x^n+1- x^n
dt = s*x^n+1 + s*y^n+1 &#8594; x^n+1 = (1+s*y^n+1)/(1+dt * s)

y^n+1 - y^n
dt = -y^n+1 + x^n+1(r-z^n+1) &#8594; y^n+1 = (1 + dt * x^n+1(r-zn+1))/(1 +dt)

z^n+1 - z^n
dt = -b*z^n+1 + x^n+1 * y^n+1 &#8594; zn+1 = (1 + dt * x^n+1 * yn+1)/(1 + b*dt )

Thank you for any help possible.

https://brainmass.com/math/ordinary-differential-equations/linearizing-lorenz-equations-implicit-euler-method-216860

## SOLUTION This solution is FREE courtesy of BrainMass!

The attached pdf file contains the explanations of only the implicit Euler linearization as we discussed in messages.

Referenced web page:
http://en.wikipedia.org/wiki/Numerical_ordinary_differential_equations

This content was COPIED from BrainMass.com - View the original, and get the already-completed solution here!