Gaps Between Primes.sagews

for n in [1..50]:
    n,is_prime(n)

(1, False) (2, True) (3, True) (4, False) (5, True) (6, False) (7, True) (8, False) (9, False) (10, False) (11, True) (12, False) (13, True) (14, False) (15, False) (16, False) (17, True) (18, False) (19, True) (20, False) (21, False) (22, False) (23, True) (24, False) (25, False) (26, False) (27, False) (28, False) (29, True) (30, False) (31, True) (32, False) (33, False) (34, False) (35, False) (36, False) (37, True) (38, False) (39, False) (40, False) (41, True) (42, False) (43, True) (44, False) (45, False) (46, False) (47, True) (48, False) (49, False) (50, False)

n=5
for k in [2..(n+1)]:
    factorial(n+1)+k,is_prime(factorial(n+1)+k)

(722, False) (723, False) (724, False) (725, False) (726, False)

n=9
for k in [2..(n+1)]:
    factorial(n+1)+k,is_prime(factorial(n+1)+k)

(39916802, False) (39916803, False) (39916804, False) (39916805, False) (39916806, False) (39916807, False) (39916808, False) (39916809, False) (39916810, False) (39916811, False)

︠497bf843-9d03-48de-9379-0c3796add6d2︠