This site is supported by donations to The OEIS Foundation.

Thanks to everyone who made a donation during our annual appeal!
To see the list of donors, or make a donation, see the OEIS Foundation home page.

 Hints (Greetings from The On-Line Encyclopedia of Integer Sequences!)
 A271837 Somos's sequence {a(7,n)} defined in comment in A018896: a(0)=a(1)= ... = a(15) = 1; for n>=16, a(n) = (a(n-1)*a(n-15)+ a(n-8)^2)/a(n-16). 7
 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 3, 4, 5, 6, 7, 8, 9, 13, 22, 38, 63, 99, 148, 212, 505, 842, 1284, 1966, 3153, 5312, 9200, 15968, 51401, 141522, 319386, 631223, 1149722, 2003800, 3442200, 9402302, 20908517, 55671036, 164685883, 466783858 (list; graph; refs; listen; history; text; internal format)
 OFFSET 0,17 LINKS Seiichi Manyama, Table of n, a(n) for n = 0..463 Eric Weisstein's World of Mathematics, Somos Sequence MATHEMATICA a[k_, n_]:=a[k, n] = If[n>2*k+1, (a[k, (n-1)]*a[k, (n-2*k-1)]+(a[k, (n-k-1)])^2 )/a[k, (n-2*k-2)], 1]; Map[a[7, #]&, Range[0, 50]] (* Peter J. C. Moses, Apr 15 2016 *) PROG (PARI) {a(n) = if(n< 16, 1, (a(n-1)*a(n-15) + a(n-8)^2)/a(n-16))}; for(n=0, 50, print1(a(n), ", ")) \\ G. C. Greubel, Feb 21 2018 (MAGMA) [n le 16 select 1 else (Self(n-1)*Self(n-15) + Self(n-8)^2 )/Self(n-16): n in [1..50]]; // G. C. Greubel, Feb 21 2018 CROSSREFS Cf. A018896, A006125, A006720, A102276, A271341, A271831, A271835, A271838, A271839. Sequence in context: A261556 A033088 A304246 * A290950 A032577 A244155 Adjacent sequences:  A271834 A271835 A271836 * A271838 A271839 A271840 KEYWORD nonn AUTHOR Vladimir Shevelev and Peter J. C. Moses, Apr 15 2016 STATUS approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

Last modified January 16 03:55 EST 2019. Contains 319184 sequences. (Running on oeis4.)