|
%I #22 Sep 08 2022 08:46:23
%S 1,18,647,23274,837217,30116538,1083358151,38970776898,1401864610177,
%T 50428155189474,1814011722210887,65253993844402458,
%U 2347329766676277601,84438617606501591178,3037442904067381004807,109263505928819214581874
%N Chebyshev T-polynomials T_n(18).
%H Colin Barker, <a href="/A322889/b322889.txt">Table of n, a(n) for n = 0..600</a>
%H Wikipedia, <a href="https://en.wikipedia.org/wiki/Chebyshev_polynomials">Chebyshev polynomials</a>.
%H <a href="/index/Ch#Cheby">Index entries for sequences related to Chebyshev polynomials.</a>
%H <a href="/index/Rec#order_02">Index entries for linear recurrences with constant coefficients</a>, signature (36, -1).
%F a(0) = 1, a(1) = 18 and a(n) = 36*a(n-1) - a(n-2) for n > 1.
%F From _Colin Barker_, Dec 30 2018: (Start)
%F G.f.: (1 - 18*x) / (1 - 36*x + x^2).
%F a(n) = ((18+sqrt(323))^(-n) * (1+(18+sqrt(323))^(2*n))) / 2.
%F (End)
%p seq(coeff(series((1-18*x)/(1-36*x+x^2),x,n+1), x, n), n = 0 .. 20); # _Muniru A Asiru_, Dec 31 2018
%t Array[ChebyshevT[#, 18] &, 16, 0] (* or *)
%t With[{k = 18}, CoefficientList[Series[(1 - k x)/(1 - 2 k x + x^2), {x, 0, 15}], x]] (* _Michael De Vlieger_, Jan 01 2019 *)
%o (PARI) {a(n) = polchebyshev(n, 1, 18)}
%o (PARI) Vec((1 - 18*x) / (1 - 36*x + x^2) + O(x^20)) \\ _Colin Barker_, Dec 30 2018
%o (GAP) a:=[1,18];; for n in [3..20] do a[n]:=36*a[n-1]-a[n-2]; od; Print(a); # _Muniru A Asiru_, Dec 31 2018
%o (Magma) I:=[1, 18]; [n le 2 select I[n] else 36*Self(n-1)-Self(n-2): n in [1..20]]; // _Vincenzo Librandi_, Jan 02 2019
%Y Column 18 of A322836.
%K nonn,easy
%O 0,2
%A _Seiichi Manyama_, Dec 29 2018
| 