OFFSET
1,1
COMMENTS
p=1, q=2(prime), a=3, b=4, c=5, s=12-+1 primes, ...
MATHEMATICA
lst={}; Do[p=n; q=p+1; a=q^2-p^2; c=q^2+p^2; b=2*p*q; ar=a*b/2; s=a+b+c; If[PrimeQ[s-1]&&PrimeQ[s+1], If[PrimeQ[q], AppendTo[lst, q]]], {n, 8!}]; lst
CROSSREFS
KEYWORD
nonn
AUTHOR
Vladimir Joseph Stephan Orlovsky, Jan 21 2009
STATUS
approved