Lorentz curve (param3d)

animation 36K

kp=xget("pixmap");xset("pixmap",1);
dr=driver();
if driver()=='Rec' then driver('X11'),end //unset memorisation
//compute lorenz differential eqn solution
y=ode(1.e-8*[1;1;1],0,0:0.003:50,'loren');
[n1,n2]=size(y)	;
wid=xget("white");
q=1:(n2/wid); //select color for each link
ds=xget("dashes");
for i=35:5:80,
  param3d(y(1,:),y(2,:),y(3,:),35,45,'X@Y@Z',[2,4]);
  xset("wwpc");  //erase
  for j=1:wid;
    xset("dashes",j);
    tt= ((j-1)*(n2/wid)+1):(j*(n2/wid));
    param3d(y(1,tt),y(2,tt),y(3,tt),i,45,'X@Y@Z',[0,0]);
  end	
  xset("wshow");  //show
end
xset("dashes",ds(1));xset("pixmap",kp);driver(dr)