login
This site is supported by donations to The OEIS Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A271831 Somos's sequence {a(6,n)} defined in comment in A018896: a(0)=a(1)= ... = a(13) = 1;  for n>=14, a(n) = (a(n-1)*a(n-13) + a(n-7)^2)/a(n-14). 8
1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 3, 4, 5, 6, 7, 8, 12, 21, 37, 62, 98, 147, 358, 609, 959, 1541, 2618, 4655, 8407, 28631, 81011, 186528, 376741, 706041, 1280174, 3598503, 8411236, 24021605, 74880071, 219318499, 580374907, 1400227135, 6308924342 (list; graph; refs; listen; history; text; internal format)
OFFSET

0,15

LINKS

Seiichi Manyama, Table of n, a(n) for n = 0..416

Eric Weisstein's World of Mathematics, Somos Sequence

MATHEMATICA

a[n_ /; 0 <= n <= 14] = 1; a[n_]:= a[n] = (a[n-1]*a[n-13] + a[n-7]^2)/a[n -14]; Table[a[n], {n, 0, 50}] (* G. C. Greubel, Feb 21 2018 *)

PROG

(PARI) {a(n) = if(n<= 14, 1, (a(n-1)*a(n-13) + a(n-7)^2)/a(n-14))};

for(n=1, 50, print1(a(n), ", ")) \\ G. C. Greubel, Feb 21 2018

(MAGMA) [n le 14 select 1 else (Self(n-1)*Self(n-13) + Self(n-7)^2)/Self(n-14): n in [1..50]]; // G. C. Greubel, Feb 21 2018

CROSSREFS

Cf. A018896, A006125, A006720, A102276, A271341, A271835, A271837, A271838, A271839.

Sequence in context: A305926 A033087 A084589 * A180628 A165306 A254116

Adjacent sequences:  A271828 A271829 A271830 * A271832 A271833 A271834

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.

License Agreements, Terms of Use, Privacy Policy. .

Last modified January 21 11:00 EST 2019. Contains 319351 sequences. (Running on oeis4.)