|
%I #28 Sep 08 2022 08:46:16
%S 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,
%T 212,505,842,1284,1966,3153,5312,9200,15968,51401,141522,319386,
%U 631223,1149722,2003800,3442200,9402302,20908517,55671036,164685883,466783858
%N 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).
%H Seiichi Manyama, <a href="/A271837/b271837.txt">Table of n, a(n) for n = 0..463</a>
%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/SomosSequence.html">Somos Sequence</a>
%t 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 *)
%o (PARI) {a(n) = if(n< 16, 1, (a(n-1)*a(n-15) + a(n-8)^2)/a(n-16))};
%o for(n=0,50, print1(a(n), ", ")) \\ _G. C. Greubel_, Feb 21 2018
%o (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
%Y Cf. A018896, A006125, A006720, A102276, A271341, A271831, A271835, A271838, A271839.
%K nonn
%O 0,17
%A _Vladimir Shevelev_ and _Peter J. C. Moses_, Apr 15 2016
|
