Returns a matrix containing the 10 000 first primes.
m = number_primes10000 ( )
a 1x10000 matrix of doubles
This function allows to update the table of primes associated to the "-primeslist" option, so that we can use a low discrepancy sequence in higher dimensions.
prarray = number_primes10000 ( ); size(prarray) lds = lowdisc_new("faure"); prarray = number_primes10000 ( ); lds = lowdisc_configure(lds,"-primeslist",prarray); lds = lowdisc_configure(lds,"-dimension",1500); [lds,next] = lowdisc_next ( lds ); next lds = lowdisc_destroy(lds);
http://primes.utm.edu/