The Sieve of Eratosthenes (see http://en.wikipedia.org/wiki/Sieve_of_Eratosthenes) is an algorithm that filters prime numbers. It iteratively identifies multiples of found primes. The multiples are, by definition, not primes and can be eliminated. This sieve is efficient for primes less than 10 million. Let's now try to find the 10001st prime number.
The first mandatory step is to create a list of natural numbers:
arange()function for that:
a = np.arange(i, i + LIM, 2)
We are not sure if this is what Eratosthenes wanted us to do, but it works. In the following code, we are passing a NumPy array and getting ...