Flows with Strange Attractors

## Select Attractor

### Lorenz Strange Attractor

dx/dt = σ(y - x)
dy/dt = x(ρ - z) - y
dz/dt = xy - βz

σ = 10
ρ = 28
β = 8/3

### Rössler Attractor

dx/dt = - y - z
dy/dt = x + ay
dz/dt = b + z(x - c)

a = 0.1
b = 0.1
c = 14

### Hénon Map

xn+1 = yn + 1 - axn2
yn+1 = bxn
(This is a 2-D attractor, so there is no z coordinate, technically)

a = 1.4
b = 0.3

### Chua's Attractor (Double Scroll Attractor)

dx/dt = α(y - x - f(x))
dy/dt = x - y + z
dz/dt = - βy
f(x) = m1x + (m0-m1)/2 * |x+1| - |x-1|)

The equation for f(x) was obtained from matlab code provided here

α = 15.6
β = 28
m0 = -1.143
m1 = -0.714

### Liénard-Van der Pol Oscillator

dx/dt = y
dy/dt = μ(1 - x2)y - x

This is a 2-D oscillator. I showed flow for various random values of μ. It also looks nice if you set z = μ

### Rabinovich-Fabrikant System

dx/dt = y(z - 1 + x2) + γx
dy/dt = x(3z + 1 - x2) + γy
dz/dt = - 2z(α + xy)

α = 0.98
γ = 0.1

### Tinkerbell Map

xn+1 = xn2 - yn2 + axn + byn
yn+1 = 2xnyn + cxn + dyn
(This is a 2-D dynamical system, so there is no z coordinate, technically)

a = 0.9
b = -0.6013
c = 2
d = 0.5

### Ikeda Map

xn+1 = 1 + u(xncos(tn) - ynsin(tn))
yn+1 = u(xnsin(tn) + yncos(tn))
tn = 0.4 - 6/(1 + xn2 + yn2)

(This is a 2-D dynamical system, so there is no z coordinate, technically)

u = 0.85

### Rainey Dynamical System

dx/dt = yz
dy/dt = xz
dz/dt = xy

Chris Rainey is a former REU/illiMATH student who made up his own set of ODEs for a dynamical system

### Fu Attractor

dx/dt = α(y + x)/z
dy/dt = x(β - x - y - z)
dz/dt = xyz/γ

α = 18
β = 42
γ = 15

My own creation (a reward for hard work)