login
a(n) = [x^n*y^n/(2*n)!] = (cosh(x)*cosh(y) + sinh(x) + sinh(y)) / (1 - sinh(x)*sinh(y)), for n >= 0.
4

%I #11 Jan 01 2019 10:00:47

%S 1,2,30,1040,68390,7322112,1161583500,256124504064,74951279707590,

%T 28111587450552320,13150168194612783620,7506080397628737060864,

%U 5134940145279960796279772,4147521174738722818202009600,3905319571990498174351316679000,4240383868623169668722969322455040,5259293239338602314365337386168258630,7389696615368349382565487880348756869120

%N a(n) = [x^n*y^n/(2*n)!] = (cosh(x)*cosh(y) + sinh(x) + sinh(y)) / (1 - sinh(x)*sinh(y)), for n >= 0.

%C a(n) = A322620(n,n) for n >= 0.

%C a(n) = A322621(n,n) for n >= 0.

%C a(n) = binomial(2*n,n) * A322195(n) for n >= 0.

%H Paul D. Hanna, <a href="/A322624/b322624.txt">Table of n, a(n) for n = 0..100</a>

%o (PARI) {A322620(n, k) = my(X=x+x*O(x^n), Y=y+y*O(y^k));

%o C = cosh(X)*cosh(Y)/(1 - sinh(X)*sinh(Y));

%o S = (sinh(X) + sinh(Y))/(1 - sinh(X)*sinh(Y));

%o (n+k)!*polcoeff(polcoeff( C + S, n, x), k, y)}

%o for(n=0,20, print1( A322620(n,n),", "))

%Y Cf. A322620, A322621, A322195.

%K nonn

%O 0,2

%A _Paul D. Hanna_, Dec 29 2018