Primes in Pi

piprimes.py

#Figure out some primes in consecutive digits of pi
from sympy import *
N = 1000
pidigits = ["3"] + list(str(pi.evalf(N)))[2:]

for end in xrange(1,N-1):
    candidate = int("".join(pidigits[:end]))
    if int(pidigits[end]) >= 5:
        rounded = candidate + 1
        if isprime(rounded):
            print "Found a rounded one with {} digits, is: {}".format(end,rounded)
    if isprime(candidate):
        print "Found one, with {} digits, is: {}".format(end,candidate)