function d=lorlambda(trans,iter,delta);
dt=0.01;
% E.g. d=lorlambda(160,4000,0.00000000000001);
x(1)=2.9; y(1)=-1.3; z(1)=25.9;
for i=2:iter
    [x(i) y(i) z(i)]=nextit(x(i-1),y(i-1),z(i-1),dt);
end
xx(1)=x(1)+delta*(rand-1/2); yy(1)=y(1)+delta*(rand-1/2); zz(1)=z(1)+delta*(rand-1/2);
for i=2:iter
    [xx(i) yy(i) zz(i)]=nextit(xx(i-1),yy(i-1),zz(i-1),dt);
end
for i=1:iter;
   % d(i)=((x(i)-xx(i))*(x(i)-xx(i))+(y(i)-yy(i))*(y(i)-yy(i))+(z(i)-zz(i))*(z(i)-zz(i))^(1/2));
   % d(i)=((x(i)-xx(i))^2+(y(i)-yy(i))^2+(z(i)-zz(i))^2);
   d(i)=(x(i)-xx(i))*(x(i)-xx(i))+(y(i)-yy(i))*(y(i)-yy(i))+(z(i)-zz(i))*(z(i)-zz(i));
end
plot(linspace(trans,iter,iter-trans+1),log(d(trans:iter)));
%polyfit(log(ratio(1,floor(q/3):ceil(2*q/3))),log(ratio(2,floor(q/3):ceil(2*q/3))),1)


function [XO YO ZO]=nextit(X,Y,Z,dt)
    r=28; s=8/3; t=10;
    x1=X+t*(Y-X)*dt/2;
	y1=Y+(X*(r-Z)-Y)*dt/2;
	z1=Z+(X*Y-s*Z)*dt/2;
	XO=X+t*(y1-x1)*dt;
	YO=Y+(x1*(r-z1)-y1)*dt;
	ZO=Z+(x1*y1-s*z1)*dt;