The OEIS is supported by the many generous donors to the OEIS Foundation.

 Hints (Greetings from The On-Line Encyclopedia of Integer Sequences!)
 A357802 Coefficients T(n,k) of x^(4*n)*r^(4*k)/(4*n)! in power series D(x,r) = 1 + r^4 * Integral S(x,r)^3 * C(x,r)^3 dx such that C(x,r)^4 - S(x,r)^4 = 1 and D(x,r)^4 - r^4*S(x,r)^4 = 1, as a triangle read by rows. 5

%I #10 Apr 12 2023 22:49:15

%S 1,0,6,0,6048,2268,0,35126784,56282688,7434504,0,679185948672,

%T 2514356038656,1409371197696,95227613712,0,33022883487154176,

%U 220415417637617664,284770675495950336,81696140755536384,3354162536029536,0,3393656235362623684608,35701050229143616880640,83079959422282198548480,54913653475645427527680,9583398717725834749440,264444869673131894208

%N Coefficients T(n,k) of x^(4*n)*r^(4*k)/(4*n)! in power series D(x,r) = 1 + r^4 * Integral S(x,r)^3 * C(x,r)^3 dx such that C(x,r)^4 - S(x,r)^4 = 1 and D(x,r)^4 - r^4*S(x,r)^4 = 1, as a triangle read by rows.

%C Equals a row reversal of triangle A357801.

%H Paul D. Hanna, <a href="/A357802/b357802.txt">Table of n, a(n) for n = 0..2925</a>

%F Generating function D(x,r) = Sum_{n>=0} Sum_{k=0..n} T(n,k) * x^(4*n) * r^(4*k) / (4*n)! and related functions S(x,r) and C(x,r) satisfy the following formulas.

%F For brevity, some formulas here will use S = S(x,r), C = C(x,r), and D = D(x,r).

%F (1.a) C(x,r)^4 - S(x,r)^4 = 1.

%F (1.b) D(x,r)^4 - r^4 * S(x,r)^4 = 1.

%F (1.c) D(x,r)^4 - r^4 * C(x,r)^4 = 1 - r^4.

%F Integral formulas.

%F (2.a) S(x,r) = Integral C(x,r)^3 * D(x,r)^3 dx.

%F (2.b) C(x,r) = 1 + Integral S(x,r)^3 * D(x,r)^3 dx.

%F (2.c) D(x,r) = 1 + r^4 * Integral S(x,r)^3 * C(x,r)^3 dx.

%F (2.d) D(x,r)^4 = 1 + r^4 * Integral 4 * S(x,r)^3 * C(x,r)^3 * D(x,r)^3 dx.

%F Derivatives.

%F (3.a) d/dx S(x,r) = C(x,r)^3 * D(x,r)^3.

%F (3.b) d/dx C(x,r) = S(x,r)^3 * D(x,r)^3.

%F (3.c) d/dx D(x,r) = r^4 * S(x,r)^3 * C(x,r)^3.

%F Exponential formulas.

%F (4.a) C + S = exp( Integral (C^2 - C*S + S^2) * D^3 dx ).

%F (4.b) D + r*S = exp( r * Integral (D^2 - r*D*S + r^2*S^2) * C^3 dx ).

%F (4.c) C - S = exp( -Integral (C^2 + C*S + S^2) * D^3 dx ).

%F (4.d) D - r*S = exp( -r * Integral (D^2 + r*D*S + r^2*S^2) * C^3 dx ).

%F (5.a) C^2 + S^2 = exp( 2 * Integral S*C * D^3 dx ).

%F (5.b) D^2 + r^2*S^2 = exp( 2*r^2 * Integral S*D * C^3 dx ).

%F (5.c) C^2 - S^2 = exp( -2 * Integral S*C * D^3 dx ).

%F (5.d) D^2 - r^2*S^2 = exp( -2*r^2 * Integral S*D * C^3 dx ).

%F Hyperbolic functions.

%F (6.a) C = sqrt(C^2 - S^2) * cosh( Integral (C^2 + S^2) * D^3 dx ).

%F (6.b) S = sqrt(C^2 - S^2) * sinh( Integral (C^2 + S^2) * D^3 dx ).

%F (6.c) D = sqrt(D^2 - r^2*S^2) * cosh( r * Integral (D^2 + r^2*S^2) * C^3 dx ).

%F (6.d) r*S = sqrt(D^2 - r^2*S^2) * sinh( r * Integral (D^2 + r^2*S^2) * C^3 dx ).

%F (7.a) C^2 = cosh( 2 * Integral S*C * D^3 dx ).

%F (7.b) S^2 = sinh( 2 * Integral S*C * D^3 dx ).

%F (7.c) D^2 = cosh( 2*r^2 * Integral S*D * C^3 dx ).

%F (7.d) r^2*S^2 = sinh( 2*r^2 * Integral S*D * C^3 dx ).

%F Other formulas.

%F (8) S(x,r) = Series_Reversion( Integral 1/((1 + x^4)*(1 + r^4*x^4))^(3/4) dx ).

%F (9.a) T(n,0) = T(n,n) = A153301(n).

%F (9.b) Sum_{k=0..n} T(n,k) = A357805(n), for n >= 0.

%F From _Paul D. Hanna_, Apr 12 2023 (Start):

%F Let F(x,r) = Integral 1/((1 + x^4)*(1 + r^4*x^4))^(3/4) dx, then

%F (10.a) S( F(x,r), r) = x,

%F (10.b) C( F(x,r), r) = (1 + x^4)^(1/4),

%F (10.c) D( F(x,r), r) = (1 + r^4*x^4)^(1/4). (End)

%e E.g.f.: D(x,r) = Sum_{n>=0} T(n,k) * x^(4*n) * r^(4*k) / (4*n)! begins:

%e D(x,r) = 1 + 6*r^4*x^4/4! + (6048*r^4 + 2268*r^8)*x^8/8! + (35126784*r^4 + 56282688*r^8 + 7434504*r^12)*x^12/12! + (679185948672*r^4 + 2514356038656*r^8 + 1409371197696*r^12 + 95227613712*r^16)*x^16/16! + (33022883487154176*r^4 + 220415417637617664*r^8 + 284770675495950336*r^12 + 81696140755536384*r^16 + 3354162536029536*r^20)*x^20/20! + (3393656235362623684608*r^4 + 35701050229143616880640*r^8 + 83079959422282198548480*r^12 + 54913653475645427527680*r^16 + 9583398717725834749440*r^20 + 264444869673131894208*r^24)*x^24/24! +

%e where D(x,r) = 1 + r^4 * Integral S(x,r)^3 * D(x,r)^3 dx.

%e TRIANGLE.

%e This triangle of coefficients T(n,k) of x^(4*n) * r^(4*k) / (4*n)! in D(x,r) for n >= 0, k = 0..n, begins:

%e n = 0: [1];

%e n = 1: [0, 6];

%e n = 2: [0, 6048, 2268];

%e n = 3: [0, 35126784, 56282688, 7434504];

%e n = 4: [0, 679185948672, 2514356038656, 1409371197696, 95227613712];

%e n = 5: [0, 33022883487154176, 220415417637617664, 284770675495950336, 81696140755536384, 3354162536029536];

%e n = 6: [0, 3393656235362623684608, 35701050229143616880640, 83079959422282198548480, 54913653475645427527680, 9583398717725834749440, 264444869673131894208]; ...

%e in which the main diagonal equals A153300.

%e RELATED SERIES.

%e S(x,r) = x + (18 + 18*r^4)*x^5/5! + (14364 + 58968*r^4 + 14364*r^8)*x^9/9! + (70203672 + 671650056*r^4 + 671650056*r^8 + 70203672*r^12)*x^13/13! + (1192064637456 + 20707300240704*r^4 + 47530354598496*r^8 + 20707300240704*r^12 + 1192064637456*r^16)*x^17/17! + (52269828456672288 + 1437626817559769760*r^4 + 5941554215913771840*r^8 + 5941554215913771840*r^12 + 1437626817559769760*r^16 + 52269828456672288*r^20)*x^21/21! + (4930307288899134335424 + 197041019249105562351744*r^4 + 1283341580573615116868160*r^8 + 2308585363008068715943680*r^12 + 1283341580573615116868160*r^16 + 197041019249105562351744*r^20 + 4930307288899134335424*r^24)*x^25/25! + ...

%e where D(x,r)^4 - r^4*S(x,r)^4 = 1.

%e C(x,r) = 1 + 6*x^4/4! + (2268 + 6048*r^4)*x^8/8! + (7434504 + 56282688*r^4 + 35126784*r^8)*x^12/12! + (95227613712 + 1409371197696*r^4 + 2514356038656*r^8 + 679185948672*r^12)*x^16/16! + (3354162536029536 + 81696140755536384*r^4 + 284770675495950336*r^8 + 220415417637617664*r^12 + 33022883487154176*r^16)*x^20/20! + (264444869673131894208 + 9583398717725834749440*r^4 + 54913653475645427527680*r^8 + 83079959422282198548480*r^12 + 35701050229143616880640*r^16 + 3393656235362623684608*r^20)*x^24/24! + ...

%e where D(x,r)^4 - r^4 * C(x,r)^4 = 1 - r^4.

%o (PARI) {T(n, k) = my(S=x, C=1, D=1); for(i=0, n,

%o S = intformal( C^3*D^3 +O(x^(4*n+4)));

%o C = 1 + intformal( S^3*D^3);

%o D = 1 + r^4*intformal( S^3*C^3); );

%o (4*n)!*polcoeff( polcoeff(D, 4*n, x), 4*k, r)}

%o for(n=0, 10, for(k=0, n, print1( T(n, k), ", ")); print(""))

%o (PARI) /* Using Series Reversion (faster) */

%o {T(n, k) = my(S = serreverse( intformal( 1/((1 + x^4)^3*(1 + r^4*x^4)^3 +O(x^(4*n+4)) )^(1/4) )) );

%o (4*n)!*polcoeff( polcoeff( (1 + r^4*S^4)^(1/4), 4*n, x), 4*k, r)}

%o for(n=0, 10, for(k=0, n, print1( T(n, k), ", ")); print(""))

%Y Cf. A153300 (diagonal), A357805 (row sums), A357800 (S(x,r)), A357801 (C(x,r)).

%Y Cf. A357542.

%K nonn,tabl

%O 0,3

%A _Paul D. Hanna_, Oct 14 2022

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

Last modified August 5 18:20 EDT 2024. Contains 374954 sequences. (Running on oeis4.)