Hi again,
Here is a small code doing the graphic job:
function X=anglePath(r, th)
cumth = cumsum(th);
X = cumsum([0 r.*cos(cumth);0 r.*sin(cumth)], 2);
end
function N=collatz(n)
N = n;
while n<>1
if modulo(n,2) == 0
n = n/2;
else
n = 3*n+1;
end
N = [n N];
end
end
a=0.08;
e=1.2;
clf
gcf().color_map=graycolormap(128);
for i=1:100
r=collatz(ceil(rand()*2000000000));
X=anglePath(r./(1+r.^1.2),a*%pi*(0.35-modulo(r,2)));
plot(0.01+X(1,:),-0.01+X(2,:));
end collatz
--
Stéphane Mottelet
Ingénieur de recherche
EA 4297 Transformations Intégrées de la Matière Renouvelable
Département Génie des Procédés Industriels
Sorbonne Universités - Université de Technologie de Compiègne
CS 60319, 60203 Compiègne cedex
Tel : +33(0)344234688
http://www.utc.fr/~mottelet
_______________________________________________
users mailing list
users@lists.scilab.org
http://lists.scilab.org/mailman/listinfo/users
