OFFSET
1,1
MAPLE
f:= proc(n) local p, m;
p:= ithprime(n);
for m from p+1 by 2 do if isprime(m^2-p) then return m fi od:
end proc:
f(1):=2:
map(f, [$1..100]); # Robert Israel, Dec 14 2017
MATHEMATICA
Table[m = p = Prime[n]; While[! PrimeQ[m^2 - p], m++]; m, {n, 60}] (* slightly modified by Robert G. Wilson v, Dec 14 2017 *)
PROG
(PARI) vector(100, n, m=prime(n); while(!isprime(m^2-prime(n)), m++); m)
CROSSREFS
KEYWORD
nonn
AUTHOR
Zak Seidov, Dec 13 2017
STATUS
approved