OFFSET
0,19
LINKS
Seiichi Manyama, Table of n, a(n) for n = 0..510
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[8, #]&, Range[0, 50]] (* Peter J. C. Moses, Apr 15 2016 *)
RecurrenceTable[{a[0]==a[1]==a[2]==a[3]==a[4]==a[5]==a[6]==a[7]==a[8]==a[9]== a[10]==a[11]==a[12]==a[13]==a[14]==a[15]==a[16]==a[17]==1, a[n]==(a[n-1]a[n-17]+ a[n-9]^2)/a[n-18]}, a, {n, 60}] (* Harvey P. Dale, Jun 30 2023 *)
PROG
(PARI) {a(n) = if(n< 18, 1, (a(n-1)*a(n-17) + a(n-9)^2)/a(n-18))};
for(n=0, 40, print1(a(n), ", ")) \\ G. C. Greubel, Feb 21 2018
(Magma) [n le 18 select 1 else (Self(n-1)*Self(n-17) + Self(n-9)^2 )/Self(n-18): n in [1..40]]; // G. C. Greubel, Feb 21 2018
CROSSREFS
KEYWORD
nonn
AUTHOR
Vladimir Shevelev and Peter J. C. Moses, Apr 15 2016
STATUS
approved