OFFSET
1,1
COMMENTS
Primes > 3 in A014209. - Klaus Purath, Dec 10 2020
MAPLE
A248697:=n->`if`(isprime(n+(n+3)^2), n+(n+3)^2, NULL): seq(A248697(n), n=1..5*10^2); # Wesley Ivan Hurt, Oct 11 2014
MATHEMATICA
f[x_] := x + (x + 3)^2;
n = 50; result = {}; counter = 0; number = 0;
While[counter < n,
value = f[number];
If[PrimeQ[value] == True, AppendTo[result, value]; counter = counter + 1];
number = number + 1]; result
Select[Table[n + (n + 3)^2, {n, 0, 300}], PrimeQ] (* Vincenzo Librandi, Oct 12 2014 *)
PROG
(Magma) [a: n in [0..250] | IsPrime(a) where a is n^2+7*n+9]; // Vincenzo Librandi, Oct 12 2014
(PARI) for(n=1, 10^3, if(isprime(n^2+7*n+9), print1(n^2+7*n+9, ", "))) \\ Derek Orr, Oct 12 2014
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Michael Savoric, Oct 11 2014
STATUS
approved