cp's OEIS Frontend

A122055 A Somos 9-type recurrence: a(n) = (3*a(n-1)*a(n-8) - a(n-4)*a(n-5))/a(n-9), with a(0)=...=a(8)=1.

%I A122055 #20 Jun 05 2024 19:33:58
%S A122055 1,1,1,1,1,1,1,1,1,2,5,14,41,121,353,989,2393,9397,49121,342793,
%T A122055 2842633,24619238,211654405,1731275594,11581792513,107195509553,
%U A122055 1126517154817,16124559341513,342648008481505,8465982933121657,213444061953471233
%N A122055 A Somos 9-type recurrence: a(n) = (3*a(n-1)*a(n-8) - a(n-4)*a(n-5))/a(n-9), with a(0)=...=a(8)=1.
%H A122055 G. C. Greubel, <a href="/A122055/b122055.txt">Table of n, a(n) for n = 0..200</a>
%F A122055 a(n) = (3*a(n-1)*a(n-8) - a(n-4)*a(n-5))/a(n-9).
%p A122055 a:= proc (n) option remember;
%p A122055    if n < 9 then 1
%p A122055    else (3*a(n-1)*a(n-8) - a(n-4)*a(n-5))/a(n-9)
%p A122055    fi;
%p A122055 end proc;
%p A122055 seq(a(n), n=0..35); # _G. C. Greubel_, Oct 03 2019
%t A122055 a[n_]:= a[n]= If[n<9, 1, (3*a[n-1]*a[n-8] -a[n-4]*a[n-5])/a[n-9]];
%t A122055 Table[a[n], {n, 0, 30}] (* modified by _G. C. Greubel_, Oct 03 2019 *)
%t A122055 nxt[{a_,b_,c_,d_,e_,f_,g_,h_,i_}]:={b,c,d,e,f,g,h,i,(3i*b-f*e)/a}; NestList[nxt,PadRight[{},9,1],30][[;;,1]] (* _Harvey P. Dale_, Jun 05 2024 *)
%o A122055 (PARI) m=35; v=concat([1,1,1,1,1,1,1,1,1], vector(m-9)); for(n=10, m, v[n] = (3*v[n-1]*v[n-8] - v[n-4]*v[n-5])/v[n-9] ); v \\ _G. C. Greubel_, Oct 03 2019
%o A122055 (Magma) [n le 10 select 1 else (3*Self(n-1)*Self(n-8) - Self(n-4)*Self(n-5))/Self(n-9): n in [1..35]]; // _G. C. Greubel_, Oct 03 2019
%o A122055 (Sage)
%o A122055 def a(n):
%o A122055     if n<10: return 1
%o A122055     else: return (3*a(n-1)*a(n-8) - a(n-4)*a(n-5))/a(n-9)
%o A122055 [a(n) for n in (0..35)] # _G. C. Greubel_, Oct 03 2019
%o A122055 (GAP)
%o A122055 a:= function(n)
%o A122055     if n<10 then return 1;
%o A122055     else return (3*a(n-1)*a(n-8) - a(n-4)*a(n-5))/a(n-9);
%o A122055     fi;
%o A122055   end;
%o A122055 List([0..35], n-> a(n) ); # _G. C. Greubel_, Oct 03 2019
%Y A122055 Cf. A122025.
%K A122055 nonn
%O A122055 0,10
%A A122055 _Roger L. Bagula_, Sep 13 2006
%E A122055 Name edited and offset changed by _G. C. Greubel_, Oct 03 2019