plot of Lorenz's strange attractor

Your browser does not support the HTML5 canvas element.

You must enable Javascript to view this page properly.

wx=0.01; wy=0.01; wz=0.01; dt=0.01; p=10; r=30; b=8/3;
x=c(); y=c(); z=c(); 
for (i in 1:50000){
  dx=(-p*wx+p*wy)*dt
  dy=(-wx*wz+r*wx-wy)*dt
  dz=(wx*wy-b*wz)*dt
  wx=wx+dx; wy=wy+dy; wz=wz+dz;
  x=c(x,wx); y=c(y,wy); z=c(z,wz); 
}
plot3d(x,y,z,type=3,col=heat.colors(length(z)))


by Akihiko YAMAZAWA