Twin Primes.sagews

p=3;
q=5;
for n in [1000000..10000000]:
    p = 6*n - 1;
    q = 6*n + 1;
    if (is_prime(p) and is_prime(q)):
        n,p,q