OFFSET
1,1
MATHEMATICA
lst={}; Do[p=Prime[n]; If[PrimeQ[(p-2)*(p+2)-2*p]&&PrimeQ[(p-2)*(p+2)+2*p]&&PrimeQ[(p-3)*(p+3)-3*p]&&PrimeQ[(p-3)*(p+3)+3*p], AppendTo[lst, p]], {n, 9!}]; lst
Select[Prime[Range[200000]], AllTrue[Flatten[{(#-2)(#+2)+{2#, -2#}, (#-3)(#+3)+ {3#, -3#}}], PrimeQ]&] (* The program uses the AllTrue function from Mathematica version 10 *) (* Harvey P. Dale, Apr 26 2015 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Vladimir Joseph Stephan Orlovsky, Jan 18 2009
STATUS
approved