|
%I #27 Feb 28 2023 23:48:23
%S 1,2,36,224,38416,2540032,4115479104,3044533460992,48656376372265216,
%T 387018647188487143424,62441634466575620320306176,
%U 5221063878050546380074377019392,8590392749565593082105293619707908096,7476351474500749779460880888573410601336832
%N a(n) = sqrt(Product_{1<=j,k<=n} (4*sin((2*j-1)*Pi/(2*n))^2 + 4*sin((2*k-1)*Pi/n)^2)).
%H Seiichi Manyama, <a href="/A341535/b341535.txt">Table of n, a(n) for n = 0..62</a>
%F a(n) ~ 2^(1/4)*(1 + sqrt(2)*(1 + (-1)^n)/2) * exp(2*G*n^2/Pi), where G is Catalan's constant A006752. - _Vaclav Kotesovec_, Feb 14 2021
%F If n is odd, a(n) = 2*A341478(n). - _Seiichi Manyama_, Feb 19 2021
%t Table[Sqrt[Product[4*Sin[(2*j - 1)*Pi/(2*n)]^2 + 4*Sin[(2*k - 1)*Pi/n]^2, {k, 1, n}, {j, 1, n}]], {n, 0, 20}] // Round (* _Vaclav Kotesovec_, Feb 14 2021 *)
%o (PARI) default(realprecision, 120);
%o a(n) = round(sqrt(prod(j=1, n, prod(k=1, n, 4*sin((2*j-1)*Pi/(2*n))^2+4*sin((2*k-1)*Pi/n)^2))));
%Y Main diagonal of A341533.
%Y Cf. A341478, A341479, A341782.
%K nonn
%O 0,2
%A _Seiichi Manyama_, Feb 13 2021
|
