    matlabpool ( 'open', 'local', 4 )

    n = 50;

    while ( n <= 500000 )
      primes = prime_number_parallel ( n );
      fprintf ( 1, '  %8d  %8d\n', n, primes );
      n = n * 10;
    end

    matlabpool ( 'close' )
