22 lines
447 B
Plaintext
22 lines
447 B
Plaintext
Halton:=proc(n,b)
|
|
local i:=n,k:=1,s:=0,r;
|
|
while i>0 do
|
|
k/=b;
|
|
i:=iquo(i,b,'r');
|
|
s+=k*r
|
|
od;
|
|
s
|
|
end;
|
|
|
|
map(Halton,[$1..10],2);
|
|
# [1/2, 1/4, 3/4, 1/8, 5/8, 3/8, 7/8, 1/16, 9/16, 5/16]
|
|
|
|
map(Halton,[$1..10],3);
|
|
# [1/3, 2/3, 1/9, 4/9, 7/9, 2/9, 5/9, 8/9, 1/27, 10/27]
|
|
|
|
map(Halton,[$1..10],4);
|
|
# [1/4, 1/2, 3/4, 1/16, 5/16, 9/16, 13/16, 1/8, 3/8, 5/8]
|
|
|
|
map(Halton,[$1..10],5);
|
|
[1/5, 2/5, 3/5, 4/5, 1/25, 6/25, 11/25, 16/25, 21/25, 2/25]
|