(* primes.ml - Don Yang (uguu.org)

   Print prime numbers.

   12/30/09
                                                                           *)

let max_prime = 1000000;;

let rec remove_primes number step bits =
   if number < (Array.length bits) then
   (
      bits.(number) <- false;
      remove_primes (number + step) step bits
   );;

let rec list_primes number bits =
   if number < (Array.length bits) then
   (
      if bits.(number) then
      (
         print_int number;
         print_newline();
         remove_primes number number bits
      );
      list_primes (number + 2) bits
   );;

list_primes 3 (Array.make max_prime true);;