|
%I #22 Sep 08 2022 08:46:25
%S 1,1,0,3,0,-1,15,0,-6,0,105,0,-45,0,1,945,0,-420,0,15,0,10395,0,-4725,
%T 0,210,0,-1,135135,0,-62370,0,3150,0,-28,0,2027025,0,-945945,0,51975,
%U 0,-630,0,1,34459425,0,-16216200,0,945945,0,-13860,0,45,0,654729075,0,-310134825,0,18918900,0,-315315,0,1485,0,-1
%N Triangle, read by rows, of Lambert's denominator polynomials related to convergents of tan(x).
%C Lambert's numerator polynomials related to convergents of tan(x), g(n, x), are given in A334824.
%H G. C. Greubel, <a href="/A334823/b334823.txt">Rows n = 0..100 of the triangle, flattened</a>
%H J.-H. Lambert, <a href="http://www.kuttaka.org/~JHL/L1768b.pdf">Mémoire sur quelques propriétés remarquables des quantités transcendantes et logarithmiques</a> (Memoir on some properties that can be traced from circular transcendent and logarithmic quantities), Histoire de l’Académie royale des sciences et belles-lettres (1761), Berlin. See <a href="http://www.bibnum.education.fr/mathematiques/theorie-des-nombres/lambert-et-l-irrationalite-de-p-1761">also</a>.
%F Equals the coefficients of the polynomials, f(n, x), defined by: (Start)
%F f(n, x) = Sum_{k=0..floor(n/2)} ((-1)^k*(2*n-2*k)!/((2*k)!*(n-2*k)!))*(x/2)^(n-2*k).
%F f(n, x) = ((2*n)!/n!)*(x/2)^n*Hypergeometric2F3(-n/2, (1-n)/2; 1/2, -n, -n+1/2; -1/x^2).
%F f(n, x) = ((-i)^n/2)*(y(n, i*x) + (-1)^n*y(n, -i*x)), where y(n, x) are the Bessel Polynomials.
%F f(n, x) = (2*n-1)*x*f(n-1, x) - f(n-2, x).
%F E.g.f. of f(n, x): cos((1 - sqrt(1-2*x*t))/2)/sqrt(1-2*x*t).
%F f(n, 1) = (-1)^n*f(n, -1) = A053983(n) = (-1)^(n+1)*A053984(-n-1) = (-1)^(n+1) * g(-n-1, 1).
%F f(n, 2) = (-1)^n*f(n, -2) = A053988(n+1). (End)
%F As a number triangle:
%F T(n, k) = i^k*(2*n-k)!*(1+(-1)^k)/(2^(n-k+1)*(k)!*(n-k)!), where i = sqrt(-1).
%F T(n, 0) = A001147(n).
%e Polynomials:
%e f(0, x) = 1;
%e f(1, x) = x;
%e f(2, x) = 3*x^2 - 1;
%e f(3, x) = 15*x^3 - 6*x;
%e f(4, x) = 105*x^4 - 45*x^2 + 1;
%e f(5, x) = 945*x^5 - 420*x^3 + 15*x;
%e f(6, x) = 10395*x^6 - 4725*x^4 + 210*x^2 - 1;
%e f(7, x) = 135135*x^7 - 62370*x^5 + 3150*x^3 - 28*x;
%e f(8, x) = 2027025*x^8 - 945945*x^6 + 51975*x^4 - 630*x^2 + 1.
%e Triangle of coefficients begins as:
%e 1;
%e 1, 0;
%e 3, 0, -1;
%e 15, 0, -6, 0;
%e 105, 0, -45, 0, 1;
%e 945, 0, -420, 0, 15, 0;
%e 10395, 0, -4725, 0, 210, 0, -1;
%e 135135, 0, -62370, 0, 3150, 0, -28, 0;
%e 2027025, 0, -945945, 0, 51975, 0, -630, 0, 1.
%p T:= (n, k) -> I^k*(2*n-k)!*(1+(-1)^k)/(2^(n-k+1)*(k)!*(n-k)!);
%p seq(seq(T(n, k), k = 0 .. n), n = 0 .. 10);
%t (* First program *)
%t y[n_, x_]:= Sqrt[2/(Pi*x)]*E^(1/x)*BesselK[-n -1/2, 1/x];
%t f[n_, k_]:= Coefficient[((-I)^n/2)*(y[n, I*x] + (-1)^n*y[n, -I*x]), x, k];
%t Table[f[n, k], {n,0,10}, {k,n,0,-1}]//Flatten
%t (* Second program *)
%t Table[ I^k*(2*n-k)!*(1+(-1)^k)/(2^(n-k+1)*(k)!*(n-k)!), {n,0,10}, {k,0,n}]//Flatten
%o (Magma)
%o C<i> := ComplexField();
%o T:= func< n, k| Round( i^k*Factorial(2*n-k)*(1+(-1)^k)/(2^(n-k+1)*Factorial(k)*Factorial(n-k)) ) >;
%o [T(n,k): k in [0..n], n in [0..10]];
%o (Sage) [[ i^k*factorial(2*n-k)*(1+(-1)^k)/(2^(n-k+1)*factorial(k)*factorial(n-k)) for k in (0..n)] for n in (0..10)]
%Y Cf. A001497, A001498, A053983, A053984, A053987, A053988, A094674, A334824.
%Y Columns k: A001147 (k=0), A001879 (k=2), A001880 (k=4), A038121 (k=6).
%K tabl,sign
%O 0,4
%A _G. C. Greubel_, May 12 2020, following a suggestion from _Michel Marcus_
|
