OFFSET
0,4
MAPLE
A094896 := proc(n)
if isprime(4*n+1) and not isprime(4*n+3) then
4*n+1;
else
0;
end if;
end proc:
seq(A094896(n), n=0..86) ; # R. J. Mathar, Nov 15 2019
MATHEMATICA
a=Table[If[PrimeQ[4*n+1]==True&&PrimeQ[4*n+3]==False, 4*n+1, 0], {n, 0, 200}]
PROG
(Magma) [IsPrime(4*n+1) and not IsPrime(4*n+3) select 4*n+1 else 0:n in [0..86]]; // Marius A. Burtea, Nov 15 2019
CROSSREFS
KEYWORD
nonn
AUTHOR
Roger L. Bagula, Jun 14 2004
STATUS
approved