function x=halton(n,s)
%x has dimension n by s and is the first n terms of the halton sequence of 
%dimension s.
p=primes(s*6);    p=p(1:s);  x=[];   
for i=1:s
   x=[x (corput(n,p(i)))'];
end