|
%I #22 Mar 13 2019 21:51:39
%S 1,0,1,0,2,5,0,16,28,61,0,272,440,662,1385,0,7936,12448,17176,24568,
%T 50521,0,353792,546560,727232,949520,1326122,2702765,0,22368256,
%U 34259968,44720896,56140288,71350336,98329108,199360981,0,1903757312,2900372480,3742967552,4600173440,5610570992,7020926600,9596075582,19391512145,0,209865342976,318605529088,408133590016,495154244608,590470281856,708137588128,877465887496,1192744081648,2404879675441
%N E.g.f. C(y,x) = cosh( Integral C(x,y)*C(y,x) dy ), where C(x,y) = cosh( Integral C(x,y)*C(y,x) dx ).
%C Row reversal of triangle A324609.
%C Related identity: (1 + sin(z))/cos(z) = exp( Integral 1/cos(z) dz ).
%H Paul D. Hanna, <a href="/A324611/b324611.txt">Table of n, a(n) for n = 0..495 terms in rows 0..30 of this triangle in flattened form.</a>
%F E.g.f. Cy = C(y,x) and related functions Sy = S(y,x), Cx = C(x,y), and Sx = S(x,y) satisfy the following relations.
%F (1a) Cx = 1 + Integral Sx * Cx*Cy dx.
%F (1b) Sx = Integral Cx * Cx*Cy dx.
%F (1c) Cy = 1 + Integral Sy * Cx*Cy dy.
%F (1d) Sy = Integral Cy * Cx*Cy dy.
%F (2a) Cx^2 - Sx^2 = 1.
%F (2b) Cy^2 - Sy^2 = 1.
%F (3a) Cx = cosh( Integral Cx*Cy dx ).
%F (3b) Sx = sinh( Integral Cx*Cy dx ).
%F (3c) Cy = cosh( Integral Cx*Cy dy ).
%F (3d) Sy = sinh( Integral Cx*Cy dy ).
%F (4a) Cx + Sx = exp( Integral Cx*Cy dx ).
%F (4b) Cy + Sy = exp( Integral Cx*Cy dy ).
%F (5a) (Cx + Sx)*(Cy + Sy) = (1 + sin(x+y))/cos(x+y).
%F (5b) (Cx + Sx)*(Cy - Sy) = (1 + sin(x-y))/cos(x-y).
%F (6a) Cx*Cy + Sx*Sy = 1/cos(x+y).
%F (6b) Cx*Sy + Sx*Cy = tan(x+y).
%F (7a) Cx*Cy = ( 1/cos(x+y) + 1/cos(x-y) )/2.
%F (7b) Sx*Sy = ( 1/cos(x+y) - 1/cos(x-y) )/2.
%F (7c) Cx*Sy = ( tan(x+y) - tan(x-y) )/2.
%F (7d) Sx*Cy = ( tan(x+y) + tan(x-y) )/2.
%F (8a) Cx*Cy = cos(x)*cos(y) / (cos(x+y)*cos(x-y)).
%F (8b) Sx*Sy = sin(x)*sin(y) / (cos(x+y)*cos(x-y)).
%F (8c) Cx*Sy = cos(y)*sin(y) / (cos(x+y)*cos(x-y)).
%F (8d) Sx*Cy = sin(x)*cos(x) / (cos(x+y)*cos(x-y)).
%F (9a) Cx + Sx = sqrt( (1 + sin(x+y))/cos(x+y) * (1 + sin(x-y))/cos(x-y) ).
%F (9b) Cy + Sy = sqrt( (1 + sin(x+y))/cos(x+y) * (1 - sin(x-y))/cos(x-y) ).
%F (9c) Cx - Sx = sqrt( (1 - sin(x+y))/cos(x+y) * (1 - sin(x-y))/cos(x-y) ).
%F (9d) Cy - Sy = sqrt( (1 - sin(x+y))/cos(x+y) * (1 + sin(x-y))/cos(x-y) ).
%F Let E(y,x) = sqrt( (1 + sin(x+y))/cos(x+y) * (1 - sin(x-y))/cos(x-y) ) then
%F (10a) E(y,x) = C(y,x) + S(y,x) where E(-y,x) = 1/E(y,x),
%F (10b) C(y,x) = (E(y,x) + E(-y,x))/2,
%F (10c) S(y,x) = (E(y,x) - E(-y,x))/2.
%F PARTICULAR ARGUMENTS.
%F E.g.f. at x = 0: C(y,x=0) = 1/cos(y).
%F E.g.f. at x = y: C(y,x=y) = cos(y)/sqrt(cos(2*y)).
%F FORMULAS INVOLVING TERMS.
%F T(n,n) = A000364(n) for n >= 0, where A000364 is the secant numbers.
%F T(n,1) = A000182(n) for n >= 1, where A000182 is the tangent numbers.
%e E.g.f.: C(y,x) = 1 + (1*y^2/2!) + (2*x^2*y^2/(2!*2!) + 5*y^4/4!) + (16*x^4*y^2/(4!*2!) + 28*x^2*y^4/(2!*4!) + 61*y^6/6!) + (272*x^6*y^2/(6!*2!) + 440*x^4*y^4/(4!*4!) + 662*x^2*y^6/(2!*6!) + 1385*y^8/8!) + (7936*x^8*y^2/(8!*2!) + 12448*x^6*y^4/(6!*4!) + 17176*x^4*y^6/(4!*6!) + 24568*x^2*y^8/(2!*8!) + 50521*y^10/10!) + (353792*x^10*y^2/(10!*2!) + 546560*x^8*y^4/(8!*4!) + 727232*x^6*y^6/(6!*6!) + 949520*x^4*y^8/(4!*8!) + 1326122*x^2*y^10/(2!*10!) + 2702765*y^12/12!) + ...
%e such that C(y,x) = cosh( Integral C(x,y)*C(y,x) dy ).
%e Explicitly,
%e C(y,x) = ( sqrt( (1 + sin(x+y))/cos(x+y) * (1 - sin(x-y))/cos(x-y) ) + sqrt( (1 - sin(x+y))/cos(x+y) * (1 + sin(x-y))/cos(x-y) ) )/2.
%e This triangle of coefficients T(n,k) of x^(2*n-2*k)*y^(2*k)/((2*n-2*k)!*(2*k)!) in e.g.f. C(y,x) begins
%e 1;
%e 0, 1;
%e 0, 2, 5;
%e 0, 16, 28, 61;
%e 0, 272, 440, 662, 1385;
%e 0, 7936, 12448, 17176, 24568, 50521;
%e 0, 353792, 546560, 727232, 949520, 1326122, 2702765;
%e 0, 22368256, 34259968, 44720896, 56140288, 71350336, 98329108, 199360981;
%e 0, 1903757312, 2900372480, 3742967552, 4600173440, 5610570992, 7020926600, 9596075582, 19391512145;
%e 0, 209865342976, 318605529088, 408133590016, 495154244608, 590470281856, 708137588128, 877465887496, 1192744081648, 2404879675441; ...
%e RELATED SERIES.
%e C(x,y) = 1 + (1*x^2/2!) + (5*x^4/4! + 2*x^2*y^2/(2!*2!)) + (61*x^6/6! + 28*x^4*y^2/(4!*2!) + 16*x^2*y^4/(2!*4!)) + (1385*x^8/8! + 662*x^6*y^2/(6!*2!) + 440*x^4*y^4/(4!*4!) + 272*x^2*y^6/(2!*6!)) + (50521*x^10/10! + 24568*x^8*y^2/(8!*2!) + 17176*x^6*y^4/(6!*4!) + 12448*x^4*y^6/(4!*6!) + 7936*x^2*y^8/(2!*8!)) + (2702765*x^12/12! + 1326122*x^10*y^2/(10!*2!) + 949520*x^8*y^4/(8!*4!) + 727232*x^6*y^6/(6!*6!) + 546560*x^4*y^8/(4!*8!) + 353792*x^2*y^10/(2!*10!)) + ...
%e such that C(x,y) = cosh( Integral C(x,y)*C(y,x) dx ).
%e S(x,y) = x + (2*x^3/3! + 1*x*y^2/2!) + (16*x^5/5! + 8*x^3*y^2/(3!*2!) + 5*x*y^4/4!) + (272*x^7/7! + 136*x^5*y^2/(5!*2!) + 94*x^3*y^4/(3!*4!) + 61*x*y^6/6!) + (7936*x^9/9! + 3968*x^7*y^2/(7!*2!) + 2840*x^5*y^4/(5!*4!) + 2108*x^3*y^6/(3!*6!) + 1385*x*y^8/8!) + (353792*x^11/11! + 176896*x^9*y^2/(9!*2!) + 128704*x^7*y^4/(7!*4!) + 100096*x^5*y^6/(5!*6!) + 76474*x^3*y^8/(3!*8!) + 50521*x*y^10/10!) + (22368256*x^13/13! + 11184128*x^11*y^2/(11!*2!) + 8211200*x^9*y^4/(9!*4!) + 6531968*x^7*y^6/(7!*6!) + 5261120*x^5*y^8/(5!*8!) + 4079408*x^3*y^10/(3!*10!) + 2702765*x*y^12/12!) + ...
%e such that C(x,y)^2 - S(x,y)^2 = 1.
%e S(y,x) = y + (1*x^2*y/2! + 2*y^3/3!) + (5*x^4*y/4! + 8*x^2*y^3/(2!*3!) + 16*y^5/5!) + (61*x^6*y/6! + 94*x^4*y^3/(4!*3!) + 136*x^2*y^5/(2!*5!) + 272*y^7/7!) + (1385*x^8*y/8! + 2108*x^6*y^3/(6!*3!) + 2840*x^4*y^5/(4!*5!) + 3968*x^2*y^7/(2!*7!) + 7936*y^9/9!) + (50521*x^10*y/10! + 76474*x^8*y^3/(8!*3!) + 100096*x^6*y^5/(6!*5!) + 128704*x^4*y^7/(4!*7!) + 176896*x^2*y^9/(2!*9!) + 353792*y^11/11!) + (2702765*x^12*y/12! + 4079408*x^10*y^3/(10!*3!) + 5261120*x^8*y^5/(8!*5!) + 6531968*x^6*y^7/(6!*7!) + 8211200*x^4*y^9/(4!*9!) + 11184128*x^2*y^11/(2!*11!) + 22368256*y^13/13!) + ...
%e such that C(y,x)^2 - S(y,x)^2 = 1.
%o (PARI) {T(n,k) = my(Cx = 1 + x*O(x^(2*n)), Cy = 1 + y*O(y^(2*n)));
%o for(i=1,2*n,
%o Cx = cosh(intformal(Cx*Cy,x));
%o Cy = cosh(intformal(Cx*Cy,y)););
%o Sx = sinh(intformal(Cx*Cy,x));
%o Sy = sinh(intformal(Cx*Cy,y));
%o (2*n-2*k)!*(2*k)! * polcoeff(polcoeff(Cy,2*n-2*k,x),2*k,y)}
%o for(n=0,10,for(k=0,n, print1( T(n,k),", "));print(""))
%Y Cf. A324612 (S(y,x)), A324609 (C(x,y)), A324610 (S(x,y)).
%Y Cf. A000364 (T(n,n)), A000182 (T(n,1)).
%K nonn,tabl
%O 0,5
%A _Paul D. Hanna_, Mar 09 2019
|
