OFFSET
1,1
FORMULA
MATHEMATICA
Ceiling[Sqrt[Times@@#]]^2&/@Partition[Prime[Range[50]], 2, 1] (* Harvey P. Dale, Aug 26 2013 *)
PROG
(Python)
from sympy import prime, primerange
def aupton(terms):
primes = list(primerange(3, prime(terms+1)+1))
return [9] + [((p+q)//2)**2 for p, q in zip(primes[:-1], primes[1:])]
print(aupton(43)) # Michael S. Branicky, Sep 16 2021
CROSSREFS
KEYWORD
nonn
AUTHOR
Reinhard Zumkeller, Oct 22 2002
STATUS
approved