OFFSET
1,1
MATHEMATICA
lst={}; Do[p=Prime[n]; If[PrimeQ[p^3-p^2-p-4], AppendTo[lst, p]], {n, 6!}]; lst
PROG
(Magma) [ n: n in PrimesUpTo(3000) | IsPrime(a) where a is n^3-n^2-n-4 ]; // Vincenzo Librandi, Dec 03 2010; corrected by Klaus Brockhaus, Dec 03 2010
CROSSREFS
KEYWORD
nonn
AUTHOR
Vladimir Joseph Stephan Orlovsky, Apr 04 2009
STATUS
approved