|
|
A141528
|
|
Expansion of x/(1 + x + 41*x^2).
|
|
2
|
|
|
0, -1, 1, 40, -81, -1559, 4880, 59039, -259119, -2161480, 12785359, 75835321, -600035040, -2509213121, 27110649761, 75767088200, -1187303728401, -1919146887799, 50598599752240, 28086422647519, -2102629012489359, 951085683941080, 85256703828122639
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,4
|
|
LINKS
|
|
|
FORMULA
|
a(n) = (-1)^(n-1)*(p^n - q^n)/(p-q), where p = (1 + sqrt(163)*i)/2, q = (1 - sqrt(163)*i)/2.
a(n) = -a(n-1) -41*a(n-2), with a(0) = 0, a(1) = -1. - G. C. Greubel, Mar 29 2021
|
|
MATHEMATICA
|
p:= (1 +Sqrt[163]*I)/2; q:= (1 -Sqrt[163]*I)/2; f[n_]:= (-1)^(n-1)*(p^n -q^n)/(p-q); Table[Simplify[f[n]], {n, 0, 30}] (* modified by G. C. Greubel, Mar 29 2021 *)
LinearRecurrence[{-1, -41}, {0, -1}, 30] (* G. C. Greubel, Mar 29 2021 *)
|
|
PROG
|
(Magma)
R<x>:=PowerSeriesRing(Rationals(), 30);
Coefficients(R!( x/(1+x+41*x^2) )); // G. C. Greubel, Mar 29 2021
(Sage)
P.<x> = PowerSeriesRing(QQ, prec)
return P( x/(1+x+41*x^2) ).list()
|
|
CROSSREFS
|
|
|
KEYWORD
|
sign
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|