|
%I #10 Dec 04 2022 18:37:14
%S 62527434837271029229
%N Primes of the form GegenbauerC[n,2,3].
%C No other terms < 10^10000. - _Robert Israel_, Apr 27 2020
%H <a href="/index/O#oneterm">Index entries for one-term sequences</a>
%p f:= gfun:-rectoproc({(-n - 3)*a(n + 3) + (24 + 6*n)*a(n + 2) + (-5 - n)*a(n + 1), a(0) = 1, a(1) = 12, a(2) = 106},a(n),remember):
%p select(isprime, map(f, [$0..4000])); # _Robert Israel_, Apr 27 2020
%t lst={};Do[p=GegenbauerC[n,2,3];If[PrimeQ[p],AppendTo[lst,p]],{n,15^3}];lst
%t Select[GegenbauerC[Range[25],2,3],PrimeQ] (* _Harvey P. Dale_, Dec 04 2022 *)
%Y Primes in A144133.
%K nonn,bref
%O 1,1
%A _Vladimir Joseph Stephan Orlovsky_, Sep 11 2008
|
