|
Search: id:A143471
|
|
|
| A143471 |
|
A McMullen transform involving x->x+1/x of a Salem polynomial gives the polynomial used to get this expansion sequence: p(x)=1 + 10 x^2 - x^3 + 45 x^4 - 8 x^5 + 120 x^6 - 27 x^7 + 210x^8 - 48 x^9 + 253 x^10 - 48 x^11 + 210 x^12 - 27 x^13 + 120 x^14 - 8 x^15 + 45 x^16 - x^17 +10 x^18 + x^20. |
|
+0
1
|
|
| 1, 0, -10, 1, 55, -12, -219, 77, 701, -351, -1900, 1277, 4494, -3966, -9485, 11058, 18342, -29012, -34057, 75053, 65836, -198845, -144194, 547462, 359314, -1548522, -937883, 4396415, 2346732, -12282817, -5272447
(list; graph; listen)
|
|
|
OFFSET
|
1,3
|
|
|
FORMULA
|
q(x)=x^10 - x^7 - x^5 - x^3 + 1; p(x)=x^10*q(x+1/x); p(x)=p(x)=1 + 10 x^2 - x^3 + 45 x^4 - 8 x^5 + 120 x^6 - 27 x^7 + 210x^8 - 48 x^9 + 253 x^10 - 48 x^11 + 210 x^12 - 27 x^13 + 120 x^14 - 8 x^15 + 45 x^16 - x^17 +10 x^18 + x^20; a(n)=Coefficient_Expansion(x^20*p(1/x)).
|
|
MATHEMATICA
|
f[x_] = x^10 - x^7 - x^5 - x^3 + 1; h[x_] = ExpandAll[x^10*f[x + 1/x]]; g[x] = ExpandAll[x^20*h[1/x]]; a = Table[SeriesCoefficient[Series[1/g[x], {x, 0, 30}], n], {n, 0, 30}]
|
|
CROSSREFS
|
Sequence in context: A049326 A146537 A050304 this_sequence A009209 A009227 A030526
Adjacent sequences: A143468 A143469 A143470 this_sequence A143472 A143473 A143474
|
|
KEYWORD
|
uned,probation,sign
|
|
AUTHOR
|
Roger L. Bagula and Gary W. Adamson (rlbagulatftn(AT)yahoo.com), Oct 24 2008
|
|
|
Search completed in 0.005 seconds
|
