%I #9 Sep 08 2022 08:45:33
%S 13,17,29,113,157,269,337,797,1613,1777,1949,2129,3613,4637,5197,6737,
%T 7069,7757,8849,11677,12113,13009,13469,14897,15889,20177,20749,23117,
%U 24977
%N Primes of the form k^2 + 13.
%H Vincenzo Librandi, <a href="/A138375/b138375.txt">Table of n, a(n) for n = 1..4500</a>
%t Intersection[Table[n^2+13,{n,0,10^2}],Prime[Range[9*10^3]]] ...or... For[i=13,i<=13,a={};Do[If[PrimeQ[n^2+i],AppendTo[a,n^2+i]],{n,0,100}];Print["n^2+",i,",",a];i++ ]
%t Select[Table[n^2+13,{n,0,7000}],PrimeQ] (* _Vincenzo Librandi_, Nov 30 2011 *)
%o (Magma) [a: n in [0..700] | IsPrime(a) where a is n^2+13]; // _Vincenzo Librandi_, Nov 30 2011
%K nonn
%O 1,1
%A _Vladimir Joseph Stephan Orlovsky_, May 07 2008