7 lines
293 B
Mathematica
7 lines
293 B
Mathematica
|
function x = corput (n)
|
||
|
|
b = dec2bin(1:n)-'0'; % generate sequence of binary numbers from 1 to n
|
||
|
|
l = size(b,2); % get number of binary digits
|
||
|
|
w = (1:l)-l-1; % 2.^w are the weights
|
||
|
|
x = b * ( 2.^w'); % matrix times vector multiplication for
|
||
|
|
end;
|