This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.
a={1,8,8}; For[n=3, n<=1000, n++, sum=Plus@@a; If[PrimeQ[sum], Print[sum]]; a=RotateLeft[a]; a[[3]]=sum]
a={1,9,9}; For[n=3, n<=1000, n++, sum=Plus@@a; If[PrimeQ[sum], Print[n]]; a=RotateLeft[a]; a[[3]]=sum]
Flatten[Position[LinearRecurrence[{1,1,1},{1,9,9},200000],?PrimeQ]]-1 (* The prime number at index 195346 (term a(39) of this sequence) has 51699 digits, so this program takes a long time to run. *) (* _Harvey P. Dale, Dec 29 2019 *)
a={1,9,9}; For[n=3, n<=1000, n++, sum=Plus@@a; If[PrimeQ[sum], Print[sum]]; a=RotateLeft[a]; a[[3]]=sum]
a:=[5,7,9];; for n in [4..40] do a[n]:=a[n-1]+a[n-2]+a[n-3]; od; a; # G. C. Greubel, Apr 23 2019
I:=[5,7,9]; [n le 3 select I[n] else Self(n-1)+Self(n-2)+Self(n-3): n in [1..40]]; // Vincenzo Librandi, Feb 04 2016
LinearRecurrence[{1, 1, 1}, {5, 7, 9}, 40]
RecurrenceTable[{a[0]==5, a[1]==7, a[2]==9, a[n]==a[n-1]+a[n-2]+a[n-3]}, a, {n, 40}]
my(x='x+O('x^40)); Vec((5+2*x-3*x^2)/(1-x-x^2-x^3)) \\ G. C. Greubel, Apr 23 2019
((5+2*x-3*x^2)/(1-x-x^2-x^3)).series(x, 40).coefficients(x, sparse=False) # G. C. Greubel, Apr 23 2019
Comments