OFFSET
1,1
MATHEMATICA
fQ[n_]:=PrimeQ[n^2-48]&&PrimeQ[n^2+48]; lst={}; Do[If[fQ@Prime[n], AppendTo[lst, Prime[n]]], {n, 7!}]; lst
Select[Prime[Range[PrimePi[4500]]], AllTrue[#^2+{48, -48}, PrimeQ]&] (* The program uses the AllTrue function from Mathematica version 10 *) (* Harvey P. Dale, Dec 08 2015 *)
PROG
(Magma) [p: p in PrimesUpTo(5000)|IsPrime(p^2-48) and IsPrime(p^2+48)] // Vincenzo Librandi, Jan 30 2011
CROSSREFS
KEYWORD
nonn
AUTHOR
Vladimir Joseph Stephan Orlovsky, Dec 23 2008
STATUS
approved