it5Q[n_]:=AllTrue[NestList[4#+3&, n, 5], PrimeQ]; Select[Prime[ Range[ 60000]], it5Q] (* The program uses the AllTrue function from Mathematica version 10 *) (* Harvey P. Dale, Jul 05 2015 *)
PROG
(Magma) [n: n in [1..10000000] | IsPrime(n) and IsPrime(4*n+3) and IsPrime(16*n+15) and IsPrime(64*n+63) and IsPrime(256*n+255) and IsPrime(1024*n+1023)]; // Vincenzo Librandi, Aug 05 2010