This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.
%I A059677 #10 Oct 26 2019 17:53:57 %S A059677 20,53,341,536,2312,2477,3380,3665,3686,4547,5009,5105,6458,6488,6731, %T A059677 6845,7499,7508,7562,7835,8411,8831,9032,9386,9764,9839,11027,11885, %U A059677 14990,19589,20498,21080,22844,24821,25220,26393,27593,29864,29921 %N A059677 Numbers n such that 1n1, 3n3, 7n7 and 9n9 are all primes. %e A059677 2312 is a term because 123121, 323123, 723127 and 923129 are all primes. %t A059677 Select[ Range[ 30000 ], PrimeQ[ ToExpression[ StringInsert[ ToString[ # ], "1", {1, -1} ] ] ] && PrimeQ[ ToExpression[ StringInsert[ ToString[ # ], "3", {1, -1} ] ] ] && PrimeQ[ ToExpression[ StringInsert[ ToString[ # ], "7", {1, -1} ] ] ] && PrimeQ[ ToExpression[ StringInsert[ ToString[ # ], "9", {1, -1} ] ] ] & ] %t A059677 enclose[n_]:=Table[FromDigits[Join[{i},IntegerDigits[n],{i}]],{i, {1,3,7,9}}]; Select[Range[30000],AllTrue[enclose[#],PrimeQ]&] (* The program uses the AllTrue function from Mathematica version 10 *) (* _Harvey P. Dale_, Mar 21 2015 *) %Y A059677 Cf. A032682, A059693, A059694. %K A059677 base,nonn %O A059677 1,1 %A A059677 _Harvey P. Dale_, Feb 05 2001 %E A059677 More terms from _Patrick De Geest_, Feb 07 2001