|
%I #14 Sep 08 2022 08:45:48
%S 1,-837,422809,-205297469,116802170481,-69673476119413,
%T 39794491851872649,-22150911964734611693,12419834337117692910305,
%U -7037064660459418136012197,3987785838055462331085793401,-2252091398491521818356890138525,1270709613993089447039294803101777
%N Expansion of 1/(1 + 837*x + 277760*x^2 + 83891456*x^3 + 7809531904*x^4).
%C Ratio limit is 496*-1.1388396294897187...;
%C the beta integer like rational pseudo-Pisot root.
%C This beta integer root is smaller than the lowest Salem number.
%H G. C. Greubel, <a href="/A167603/b167603.txt">Table of n, a(n) for n = 0..100</a>
%H <a href="/index/Rec#order_04">Index entries for linear recurrences with constant coefficients</a>, signature (-837, -277760, -83891456, -7809531904).
%F a(n+4) + 837*a(n+3) + 277760*a(n+2) + 83891456*a(n+1) + 7809531904*a(n) = 0. - _G. C. Greubel_, Jun 17 2016
%t LinearRecurrence[{-837, -277760, -83891456, -7809531904}, {1, -837, 422809, -205297469}, 50] (* _G. C. Greubel_, Jun 17 2016 *)
%o (PARI) x='x+O('x^50); Vec(1/(1 + 837*x + 277760*x^2 + 83891456*x^3 + 7809531904*x^4)) \\ _G. C. Greubel_, Nov 03 2018
%o (Magma) m:=50; R<x>:=PowerSeriesRing(Integers(), m); Coefficients(R!(1/(1 + 837*x + 277760*x^2 + 83891456*x^3 + 7809531904*x^4))); // _G. C. Greubel_, Nov 03 2018
%Y Cf. A143471, A143478.
%K sign,easy
%O 0,2
%A _Roger L. Bagula_, Nov 06 2009
%E New name by _Franck Maminirina Ramaharo_, Nov 02 2018
| 