attractor/16 Lorenz Ⅲ (stop_motion)

> a = 10, b = 28, c = 6, ``` def lorenz_III(a, b, c): h = lambda x, y: sqrt(x*x + y*y) return lambda x, y, z:( ((-a + 1)*x + (a - b + z)*y)/3 + ((1 - a)*(x*x - y*y) + 2*(b + a - z)*x*y)/(3*h(x, y)), ((b - a - z)*x - (a + 1)*y)/3 + (2*(a - 1)*x*y + (b + a - z)*(x*x - y*y))/(3*h(x, y)), (3*x*x*y - y*y*y)/2 - c*z ) ``` ## References > Lorenz-Attraktor > http://page.math.tu-berlin.de/~gunn/Files/mvws13html/html2420/lorenz2420.html > > R. Miranda & E. Stone, The proto-Lorenz system, Physics Letters A 178 (1993) 105-113.