cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-5 of 5 results.

A340322 Decimal expansion of Integral_{x=0..Pi/2, y=0..Pi/2, z=0..Pi/2} log(4*cos(x)^2 + 4*cos(y)^2 + 4*cos(z)^2) dz dy dx.

Original entry on oeis.org

6, 4, 8, 5, 6, 9, 6, 4, 6, 5, 2, 1, 8, 4, 9, 7, 6, 9, 3, 7, 0, 8, 5, 8, 1, 3, 7, 2, 1, 0, 3, 3, 1, 5, 7, 6, 4, 1, 5, 2, 2, 6, 6, 3, 2, 5, 6, 1, 7, 9, 7, 6, 3, 1, 6, 8, 3, 1, 7, 3, 8, 8, 4, 2, 4, 5, 2, 5, 5, 5, 2, 3, 8, 7, 8, 4, 0, 0, 5, 6, 7, 8, 5, 4, 1, 8, 8, 8, 8, 7, 6, 4, 1, 9, 6, 8, 5, 7, 5, 5, 3, 9, 1, 7, 4
Offset: 1

Views

Author

Vaclav Kotesovec, Jan 04 2021

Keywords

Comments

Integral_{x=0..Pi/2, y=0..Pi/2} log(4*cos(x)^2 + 4*cos(y)^2) dy dx = G*Pi, where G is Catalan's constant A006752.

Examples

			6.485696465218497693708581372103315764152266325617976316831738842452555238784...

Links

Crossrefs

Cf. A071763, A097469, A340182, A340350, A340421, A340422.

Programs

  • Maple
    evalf(Integrate(log(4*cos(x)^2 + 4*cos(y)^2 + 4*cos(z)^2), x = 0..Pi/2, y = 0..Pi/2, z = 0..Pi/2));
  • PARI
    intnum(x = 0, Pi/2, intnum(y = 0, Pi/2, intnum(z = 0, Pi/2, log(4*cos(x)^2 + 4*cos(y)^2 + 4*cos(z)^2)))) \\ 20 valid digits

Formula

Equals limit_{n->infinity} Pi^3 * log(A340182(n)) / (8*n^3).
Equals Pi^3 * log(2)/8 + Integral_{x=0..Pi/2, y=0..Pi/2, z=0..Pi/2} log(3 + cos(2*x) + cos(2*y) + cos(2*z)) dz dy dx.

A071763 Number of spanning trees in n X n X n grid.

Original entry on oeis.org

1, 384, 8193540096000, 172685928902844729688524604506636288, 77746347057132811936046563068332100246216273086593103906734080000000000000
Offset: 1

Views

Author

Sharon Sela (sharonsela(AT)hotmail.com), Jun 04 2002

Keywords

Links

Crossrefs

Cf. A007341, A340182.

Programs

  • Mathematica
    Table[2^(n^3 - 1)/n^3 Product[Piecewise[{{1, i == j == k == 0}}, 3 - Cos[Pi i/n] - Cos[Pi j/n] - Cos[Pi k/n]], {i, 0, n - 1}, {j, 0, n - 1}, {k, 0, n - 1}], {n, 12}] // Round

Formula

a(n) = 2^(n^3-1) / n^3 * Product_{n1=0..n-1 n2=0..n-1 n3=0..n-1} (3- cos(Pi*n1/n) - cos(Pi*n2/n) - cos(Pi*n3/n) ) where n1, n2, n3 are not all 0.
Limit_{n->infinity} a(n)^(1/n^3) = exp(8 * A340322 / Pi^3) = 5.330202889205167421134597996649659520108446730592285502966091902480522584119... - Vaclav Kotesovec, Jan 05 2021

A340183 a(n) = Product_{1<=j,k,m<=n-1} (4*sin(j*Pi/(2*n))^2 + 4*sin(k*Pi/(2*n))^2 + 4*sin(m*Pi/(2*n))^2).

Original entry on oeis.org

1, 6, 1157625, 170875128460147163136, 448524809573174705684873233798538664686384705625
Offset: 1

Views

Author

Seiichi Manyama, Dec 31 2020

Keywords

Comments

(a(n)/(n*3^(n-1)))^(1/3) is an integer.

Crossrefs

Cf. A007341, A124647, A340181, A340182.

Programs

  • Mathematica
    Round[Table[2^((n-1)^3)* Product[3 - Cos[j*Pi/n] - Cos[k*Pi/n] - Cos[m*Pi/n], {j, 1, n-1}, {k, 1, n-1}, {m, 1, n-1}], {n, 1, 5}]] (* Vaclav Kotesovec, Jan 04 2021 *)
  • PARI
    default(realprecision, 500);
{a(n) = round(prod(j=1, n-1, prod(k=1, n-1, prod(m=1, n-1, 4*sin(j*Pi/(2*n))^2+4*sin(k*Pi/(2*n))^2+4*sin(m*Pi/(2*n))^2))))}

Formula

a(n) = Product_{1<=i,j,k<=n-1} (4*f(i*Pi/(2*n))^2 + 4*g(j*Pi/(2*n))^2 + 4*h(k*Pi/(2*n))^2), where f(x), g(x) and h(x) are sin(x) or cos(x).
Limit_{n->infinity} a(n)^(1/n^3) = exp(8*A340322/Pi^3). - Vaclav Kotesovec, Jan 05 2021

A340181 a(n) = Product_{1<=j,k,m<=n} (4*sin(j*Pi/(2*n+1))^2 + 4*sin(k*Pi/(2*n+1))^2 + 4*sin(m*Pi/(2*n+1))^2).

Original entry on oeis.org

1, 9, 7486875, 14334918272193811385583, 1483160703050490588200236172057973908184332257091136
Offset: 0

Views

Author

Seiichi Manyama, Dec 31 2020

Keywords

Comments

(a(n)/((2n + 1)*3^n))^(1/3) is an integer.

Crossrefs

Cf. A124647, A127605, A340182, A340183.

Programs

  • Mathematica
    Round[Table[2^(n^3)* Product[3 - Cos[2*j*Pi/(2*n + 1)] - Cos[2*k*Pi/(2*n + 1)] - Cos[2*m*Pi/(2*n + 1)], {j, 1, n}, {k, 1, n}, {m, 1, n}], {n, 0, 5}]] (* Vaclav Kotesovec, Jan 04 2021 *)
  • PARI
    default(realprecision, 500);
{a(n) = round(prod(j=1, n, prod(k=1, n, prod(m=1, n, 4*sin(j*Pi/(2*n+1))^2+4*sin(k*Pi/(2*n+1))^2+4*sin(m*Pi/(2*n+1))^2))))}

Formula

Limit_{n->infinity} a(n)^(1/n^3) = exp(8*A340322/Pi^3). - Vaclav Kotesovec, Jan 05 2021

A340396 a(n) = 2^(n^2 - 1) * Product_{j=1..n, k=1..n} (1 + sin(Pi*j/n)^2 + sin(Pi*k/n)^2).

Original entry on oeis.org

0, 1, 96, 93789, 1244160000, 241885578271872, 700566272328037500000, 30323548995402141685610526683, 19627362048402730985830806120284160000, 189995156103157091521654945902925881881155376920, 27506190205802587152768139358989866456457087869970721213256
Offset: 0

Views

Author

Vaclav Kotesovec, Jan 06 2021

Keywords

Crossrefs

Cf. A004003, A007726, A065072, A067518, A127605, A340052, A340165, A340182, A340183, A340293.

Programs

  • Mathematica
    Table[2^(n^2 - 1) * Product[1 + Sin[Pi*j/n]^2 + Sin[Pi*k/n]^2, {j, 1, n}, {k, 1, n}], {n, 0, 10}] // Round

Formula

a(n) = 2^(n^2-1) * Product_{j=1..n, k=1..n} (3 - cos(Pi*j/n)^2 - cos(Pi*k/n)^2).
a(n) = 2^(n^2-1) * Product_{j=1..n, k=1..n} (2-cos(2*Pi*j/n)/2-cos(2*Pi*k/n)/2).
a(n) ~ 2^(n^2-1) * exp(4*c*n^2/Pi^2), where c = Integral_{x=0..Pi/2, y=0..Pi/2} log(1 + sin(x)^2 + sin(y)^2) dy dx = -Pi^2*(log(2) + log(sqrt(2)-1)/2) + Pi * Integral_{x=0..Pi/2} log(1 + sqrt(1 + 1/(1 + sin(x)^2))) dx = A340421 = 1.627008991085721315763766677017604437985734719035793082916212355323520649...
Showing 1-5 of 5 results.

Started in 1964 by Neil J. A. Sloane | Maintained by The OEIS Foundation Inc.

Content is available under The OEIS End-User License Agreement.