login
a(n) = value of Chebyshev T-polynomial T_n(20).
2

%I #26 Jan 20 2019 06:47:07

%S 1,20,799,31940,1276801,51040100,2040327199,81562047860,3260441587201,

%T 130336101440180,5210183616019999,208277008539359780,

%U 8325870157958371201,332826529309795488260,13304735302233861159199,531856585560044650879700

%N a(n) = value of Chebyshev T-polynomial T_n(20).

%H Seiichi Manyama, <a href="/A322890/b322890.txt">Table of n, a(n) for n = 0..624</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 (40, -1).

%F a(0) = 1, a(1) = 20 and a(n) = 40*a(n-1) - a(n-2) for n > 1.

%F From _Colin Barker_, Dec 30 2018: (Start)

%F G.f.: (1 - 20*x) / (1 - 40*x + x^2).

%F a(n) = ((20+sqrt(399))^(-n) * (1+(20+sqrt(399))^(2*n))) / 2.

%F (End)

%p seq(coeff(series((1-20*x)/(1-40*x+x^2),x,n+1), x, n), n = 0 .. 20); # _Muniru A Asiru_, Dec 31 2018

%t CoefficientList[Series[(1 - 20 x)/(1 - 40 x + x^2), {x, 0, 15}], x] (* or *)

%t Array[ChebyshevT[#, 20] &, 16, 0] (* _Michael De Vlieger_, Jan 01 2019 *)

%o (PARI) {a(n) = polchebyshev(n, 1, 20)}

%o (PARI) Vec((1 - 20*x) / (1 - 40*x + x^2) + O(x^20)) \\ _Colin Barker_, Dec 30 2018

%o (GAP) a:=[1,20];; for n in [3..20] do a[n]:=40*a[n-1]-a[n-2]; od; Print(a); # _Muniru A Asiru_, Dec 31 2018

%Y Column 20 of A322836.

%Y Cf. A041758.

%K nonn,easy

%O 0,2

%A _Seiichi Manyama_, Dec 29 2018