a b r dt

10 iterations per frame

Drag mouse to rotate model. Hold shift key or use mouse wheel to zoom it (canvas is matched to the browser window). α determines point brightness.
Dynamics of N = 4 194 304 independent points is plotted. GPU based simulation of the Lorenz strange attractor. Equations of the particle motion are
    x' = a(y - x),
    y' = x(r - z) - y,
    z' = xy - bz.

The script uses the second order Runge-Kutta integration algorithm (x is vector (x,y,z) below)
    xi+1 = xi + [F(xi + F(xi )dt/2)]dt.
It makes 10 iterations per frame with time step dt.
Simulations on GPU
updated 28 Aug 2013