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-10 of 85 results. Next

A000796 Decimal expansion of Pi (or digits of Pi).

Original entry on oeis.org

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

Views

Author

Keywords

Comments

Sometimes called Archimedes's constant.
Ratio of a circle's circumference to its diameter.
Also area of a circle with radius 1.
Also surface area of a sphere with diameter 1.
A useful mnemonic for remembering the first few terms: How I want a drink, alcoholic of course, after the heavy lectures involving quantum mechanics ...
Also ratio of surface area of sphere to one of the faces of the circumscribed cube. Also ratio of volume of a sphere to one of the six inscribed pyramids in the circumscribed cube. - Omar E. Pol, Aug 09 2012
Also surface area of a quarter of a sphere of radius 1. - Omar E. Pol, Oct 03 2013
Also the area under the peak-shaped even function f(x)=1/cosh(x). Proof: for the upper half of the integral, write f(x) = (2*exp(-x))/(1+exp(-2x)) = 2*Sum_{k>=0} (-1)^k*exp(-(2k+1)*x) and integrate term by term from zero to infinity. The result is twice the Gregory series for Pi/4. - Stanislav Sykora, Oct 31 2013
A curiosity: a 144 X 144 magic square of 7th powers was recently constructed by Toshihiro Shirakawa. The magic sum = 3141592653589793238462643383279502884197169399375105, which is the concatenation of the first 52 digits of Pi. See the MultiMagic Squares link for details. - Christian Boyer, Dec 13 2013 [Comment revised by N. J. A. Sloane, Aug 27 2014]
x*Pi is also the surface area of a sphere whose diameter equals the square root of x. - Omar E. Pol, Dec 25 2013
Also diameter of a sphere whose surface area equals the volume of the circumscribed cube. - Omar E. Pol, Jan 13 2014
From Daniel Forgues, Mar 20 2015: (Start)
An interesting anecdote about the base-10 representation of Pi, with 3 (integer part) as first (index 1) digit:
358 0
359 3
360 6
361 0
362 0
And the circle is customarily subdivided into 360 degrees (although Pi radians yields half the circle)...
(End)
Sometimes referred to as Archimedes's constant, because the Greek mathematician computed lower and upper bounds of Pi by drawing regular polygons inside and outside a circle. In Germany it was called the Ludolphian number until the early 20th century after the Dutch mathematician Ludolph van Ceulen (1540-1610), who calculated up to 35 digits of Pi in the late 16th century. - Martin Renner, Sep 07 2016
As of the beginning of 2019 more than 22 trillion decimal digits of Pi are known. See the Wikipedia article "Chronology of computation of Pi". - Harvey P. Dale, Jan 23 2019
On March 14, 2019, Emma Haruka Iwao announced the calculation of 31.4 trillion digits of Pi using Google Cloud's infrastructure. - David Radcliffe, Apr 10 2019
Also volume of three quarters of a sphere of radius 1. - Omar E. Pol, Aug 16 2019
On August 5, 2021, researchers from the University of Applied Sciences of the Grisons in Switzerland announced they had calculated 62.8 trillion digits. Guinness World Records has not verified this yet. - Alonso del Arte, Aug 23 2021
The Hermite-Lindemann (1882) theorem states, that if z is a nonzero algebraic number, then e^z is a transcendent number. The transcendence of Pi then results from Euler's relation: e^(i*Pi) = -1. - Peter Luschny, Jul 21 2023
The 10000 words of the book "Not A Wake" by Michael Keith, written in Pilish, match in length the first 10000 digits of Pi. - Paolo Xausa, Aug 07 2025

Examples

			3.1415926535897932384626433832795028841971693993751058209749445923078164062\
862089986280348253421170679821480865132823066470938446095505822317253594081\
284811174502841027019385211055596446229489549303819...
		

References

  • Mohammad K. Azarian, A Summary of Mathematical Works of Ghiyath ud-din Jamshid Kashani, Journal of Recreational Mathematics, Vol. 29(1), pp. 32-42, 1998.
  • J. Arndt & C. Haenel, Pi Unleashed, Springer NY 2001.
  • P. Beckmann, A History of Pi, Golem Press, Boulder, CO, 1977.
  • Florian Cajori, A History of Mathematical Notations, Dover edition (2012), par. 396.
  • John H. Conway and Richard K. Guy, The Book of Numbers, New York: Springer-Verlag, 1996. See pp. 24, 237-239.
  • J.-P. Delahaye, Le fascinant nombre pi, Pour la Science, Paris 1997.
  • P. Eyard and J.-P. Lafon, The Number Pi, Amer. Math. Soc., 2004.
  • S. R. Finch, Mathematical Constants, Encyclopedia of Mathematics and its Applications, vol. 94, Cambridge University Press, Section 1.4.
  • Le Petit Archimede, Special Issue On Pi, Supplement to No. 64-5, May 1980 ADCS Amiens.
  • Clifford A. Pickover, A Passion for Mathematics, Wiley, 2005; see p. 31.
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
  • Jerome Spanier and Keith B. Oldham, "Atlas of Functions", Hemisphere Publishing Corp., 1987, chapter 1, equations 1:7:1, 1:7:2 at pages 12-13.
  • David Wells, The Penguin Dictionary of Curious and Interesting Numbers. Penguin Books, NY, 1986, Revised edition 1987. See pp. 48-55.

Crossrefs

Cf. A001203 (continued fraction).
Pi in base b: A004601 (b=2), A004602 (b=3), A004603 (b=4), A004604 (b=5), A004605 (b=6), A004606 (b=7), A006941 (b=8), A004608 (b=9), this sequence (b=10), A068436 (b=11), A068437 (b=12), A068438 (b=13), A068439 (b=14), A068440 (b=15), A062964 (b=16), A224750 (b=26), A224751 (b=27), A060707 (b=60). - Jason Kimberley, Dec 06 2012
Decimal expansions of expressions involving Pi: A002388 (Pi^2), A003881 (Pi/4), A013661 (Pi^2/6), A019692 (2*Pi=tau), A019727 (sqrt(2*Pi)), A059956 (6/Pi^2), A060294 (2/Pi), A091925 (Pi^3), A092425 (Pi^4), A092731 (Pi^5), A092732 (Pi^6), A092735 (Pi^7), A092736 (Pi^8), A163973 (Pi/log(2)).
Cf. A001901 (Pi/2; Wallis), A002736 (Pi^2/18; Euler), A007514 (Pi), A048581 (Pi; BBP), A054387 (Pi; Newton), A092798 (Pi/2), A096954 (Pi/4; Machin), A097486 (Pi), A122214 (Pi/2), A133766 (Pi/4 - 1/2), A133767 (5/6 - Pi/4), A166107 (Pi; MGL).

Programs

  • Haskell
    -- see link: Literate Programs
    import Data.Char (digitToInt)
    a000796 n = a000796_list (n + 1) !! (n + 1)
    a000796_list len = map digitToInt $ show $ machin' `div` (10 ^ 10) where
       machin' = 4 * (4 * arccot 5 unity - arccot 239 unity)
       unity = 10 ^ (len + 10)
       arccot x unity = arccot' x unity 0 (unity `div` x) 1 1 where
         arccot' x unity summa xpow n sign
          | term == 0 = summa
          | otherwise = arccot'
            x unity (summa + sign * term) (xpow `div` x ^ 2) (n + 2) (- sign)
          where term = xpow `div` n
    -- Reinhard Zumkeller, Nov 24 2012
    
  • Haskell
    -- See Niemeijer link and also Gibbons link.
    a000796 n = a000796_list !! (n-1) :: Int
    a000796_list = map fromInteger $ piStream (1, 0, 1)
       [(n, a*d, d) | (n, d, a) <- map (\k -> (k, 2 * k + 1, 2)) [1..]] where
       piStream z xs'@(x:xs)
         | lb /= approx z 4 = piStream (mult z x) xs
         | otherwise = lb : piStream (mult (10, -10 * lb, 1) z) xs'
         where lb = approx z 3
               approx (a, b, c) n = div (a * n + b) c
               mult (a, b, c) (d, e, f) = (a * d, a * e + b * f, c * f)
    -- Reinhard Zumkeller, Jul 14 2013, Jun 12 2013
    
  • Macsyma
    py(x) := if equal(6,6+x^2) then 2*x else (py(x:x/3),3*%%-4*(%%-x)^3); py(3.); py(dfloat(%)); block([bfprecision:35], py(bfloat(%))) /* Bill Gosper, Sep 09 2002 */
    
  • Magma
    pi:=Pi(RealField(110)); Reverse(Intseq(Floor(10^105*pi))); // Bruno Berselli, Mar 12 2013
    
  • Maple
    Digits := 110: Pi*10^104:
    ListTools:-Reverse(convert(floor(%), base, 10)); # Peter Luschny, Oct 29 2019
  • Mathematica
    RealDigits[ N[ Pi, 105]] [[1]]
    Table[ResourceFunction["NthDigit"][Pi, n], {n, 1, 102}] (* Joan Ludevid, Jun 22 2022; easy to compute a(10000000)=7 with this function; requires Mathematica 12.0+ *)
  • PARI
    { default(realprecision, 20080); x=Pi; for (n=1, 20000, d=floor(x); x=(x-d)*10; write("b000796.txt", n, " ", d)); } \\ Harry J. Smith, Apr 15 2009
    
  • PARI
    A796=[]; A000796(n)={if(n>#A796, localprec(n*6\5+29); A796=digits(Pi\.1^(precision(Pi)-3))); A796[n]} \\ NOTE: as the other programs, this returns the n-th term of the sequence, with n = 1, 2, 3, ... and not n = 1, 0, -1, -2, .... - M. F. Hasler, Jun 21 2022
    
  • PARI
    first(n)= default(realprecision, n+10); digits(floor(Pi*10^(n-1))) \\ David A. Corneth, Jun 21 2022
    
  • PARI
    lista(nn, p=20)= {my(u=10^(nn+p+1), f(x, u)=my(n=1, q=u\x, r=q, s=1, t); while(t=(q\=(x*x))\(n+=2), r+=(s=-s)*t); r*4); digits((4*f(5, u)-f(239, u))\10^(p+2)); } \\ Machin-like, with p > the maximal number of consecutive 9-digits to be expected (A048940) - Ruud H.G. van Tol, Dec 26 2024
    
  • Python
    from sympy import pi, N; print(N(pi, 1000)) # David Radcliffe, Apr 10 2019
    
  • Python
    from mpmath import mp
    def A000796(n):
        if n >= len(A000796.str): mp.dps = n*6//5+50; A000796.str = str(mp.pi-5/mp.mpf(10)**mp.dps)
        return int(A000796.str[n if n>1 else 0])
    A000796.str = '' # M. F. Hasler, Jun 21 2022
    
  • SageMath
    m=125
    x=numerical_approx(pi, digits=m+5)
    a=[ZZ(i) for i in x.str(skip_zeroes=True) if i.isdigit()]
    a[:m] # G. C. Greubel, Jul 18 2023

Formula

Pi = 4*Sum_{k>=0} (-1)^k/(2k+1) [Madhava-Gregory-Leibniz, 1450-1671]. - N. J. A. Sloane, Feb 27 2013
From Johannes W. Meijer, Mar 10 2013: (Start)
2/Pi = (sqrt(2)/2) * (sqrt(2 + sqrt(2))/2) * (sqrt(2 + sqrt(2 + sqrt(2)))/2) * ... [Viete, 1593]
2/Pi = Product_{k>=1} (4*k^2-1)/(4*k^2). [Wallis, 1655]
Pi = 3*sqrt(3)/4 + 24*(1/12 - Sum_{n>=2} (2*n-2)!/((n-1)!^2*(2*n-3)*(2*n+1)*2^(4*n-2))). [Newton, 1666]
Pi/4 = 4*arctan(1/5) - arctan(1/239). [Machin, 1706]
Pi^2/6 = 3*Sum_{n>=1} 1/(n^2*binomial(2*n,n)). [Euler, 1748]
1/Pi = (2*sqrt(2)/9801) * Sum_{n>=0} (4*n)!*(1103+26390*n)/((n!)^4*396^(4*n)). [Ramanujan, 1914]
1/Pi = 12*Sum_{n>=0} (-1)^n*(6*n)!*(13591409 + 545140134*n)/((3*n)!*(n!)^3*(640320^3)^(n+1/2)). [David and Gregory Chudnovsky, 1989]
Pi = Sum_{n>=0} (1/16^n) * (4/(8*n+1) - 2/(8*n+4) - 1/(8*n+5) - 1/(8*n+6)). [Bailey-Borwein-Plouffe, 1989] (End)
Pi = 4 * Sum_{k>=0} 1/(4*k+1) - 1/(4*k+3). - Alexander R. Povolotsky, Dec 25 2008
Pi = 4*sqrt(-1*(Sum_{n>=0} (i^(2*n+1))/(2*n+1))^2). - Alexander R. Povolotsky, Jan 25 2009
Pi = Integral_{x=-oo..oo} dx/(1+x^2). - Mats Granvik and Gary W. Adamson, Sep 23 2012
Pi - 2 = 1/1 + 1/3 - 1/6 - 1/10 + 1/15 + 1/21 - 1/28 - 1/36 + 1/45 + ... [Jonas Castillo Toloza, 2007], that is, Pi - 2 = Sum_{n>=1} (1/((-1)^floor((n-1)/2)*(n^2+n)/2)). - José de Jesús Camacho Medina, Jan 20 2014
Pi = 3 * Product_{t=img(r),r=(1/2+i*t) root of zeta function} (9+4*t^2)/(1+4*t^2) <=> RH is true. - Dimitris Valianatos, May 05 2016
From Ilya Gutkovskiy, Aug 07 2016: (Start)
Pi = Sum_{k>=1} (3^k - 1)*zeta(k+1)/4^k.
Pi = 2*Product_{k>=2} sec(Pi/2^k).
Pi = 2*Integral_{x>=0} sin(x)/x dx. (End)
Pi = 2^{k + 1}*arctan(sqrt(2 - a_{k - 1})/a_k) at k >= 2, where a_k = sqrt(2 + a_{k - 1}) and a_1 = sqrt(2). - Sanjar Abrarov, Feb 07 2017
Pi = Integral_{x = 0..2} sqrt(x/(2 - x)) dx. - Arkadiusz Wesolowski, Nov 20 2017
Pi = lim_{n->oo} 2/n * Sum_{m=1,n} ( sqrt( (n+1)^2 - m^2 ) - sqrt( n^2 - m^2 ) ). - Dimitri Papadopoulos, May 31 2019
From Peter Bala, Oct 29 2019: (Start)
Pi = Sum_{n >= 0} 2^(n+1)/( binomial(2*n,n)*(2*n + 1) ) - Euler.
More generally, Pi = (4^x)*x!/(2*x)! * Sum_{n >= 0} 2^(n+1)*(n+x)!*(n+2*x)!/(2*n+2*x+1)! = 2*4^x*x!^2/(2*x+1)! * hypergeom([2*x+1,1], [x+3/2], 1/2), valid for complex x not in {-1,-3/2,-2,-5/2,...}. Here, x! is shorthand notation for the function Gamma(x+1). This identity may be proved using Gauss's second summation theorem.
Setting x = 3/4 and x = -1/4 (resp. x = 1/4 and x = -3/4) in the above identity leads to series representations for the constant A085565 (resp. A076390). (End)
Pi = Im(log(-i^i)) = log(i^i)*(-2). - Peter Luschny, Oct 29 2019
From Amiram Eldar, Aug 15 2020: (Start)
Equals 2 + Integral_{x=0..1} arccos(x)^2 dx.
Equals Integral_{x=0..oo} log(1 + 1/x^2) dx.
Equals Integral_{x=0..oo} log(1 + x^2)/x^2 dx.
Equals Integral_{x=-oo..oo} exp(x/2)/(exp(x) + 1) dx. (End)
Equals 4*(1/2)!^2 = 4*Gamma(3/2)^2. - Gary W. Adamson, Aug 23 2021
From Peter Bala, Dec 08 2021: (Start)
Pi = 32*Sum_{n >= 1} (-1)^n*n^2/((4*n^2 - 1)*(4*n^2 - 9))= 384*Sum_{n >= 1} (-1)^(n+1)*n^2/((4*n^2 - 1)*(4*n^2 - 9)*(4*n^2 - 25)).
More generally, it appears that for k = 1,2,3,..., Pi = 16*(2*k)!*Sum_{n >= 1} (-1)^(n+k+1)*n^2/((4*n^2 - 1)* ... *(4*n^2 - (2*k+1)^2)).
Pi = 32*Sum_{n >= 1} (-1)^(n+1)*n^2/(4*n^2 - 1)^2 = 768*Sum_{n >= 1} (-1)^(n+1)*n^2/((4*n^2 - 1)^2*(4*n^2 - 9)^2).
More generally, it appears that for k = 0,1,2,..., Pi = 16*Catalan(k)*(2*k)!*(2*k+2)!*Sum_{n >= 1} (-1)^(n+1)*n^2/((4*n^2 - 1)^2* ... *(4*n^2 - (2*k+1)^2)^2).
Pi = (2^8)*Sum_{n >= 1} (-1)^(n+1)*n^2/(4*n^2 - 1)^4 = (2^17)*(3^5)*Sum_{n >= 2} (-1)^n*n^2*(n^2 - 1)/((4*n^2 - 1)^4*(4*n^2 - 9)^4) = (2^27)*(3^5)*(5^5)* Sum_{n >= 3} (-1)^(n+1)*n^2*(n^2 - 1)*(n^2 - 4)/((4*n^2 - 1)^4*(4*n^2 - 9)^4*(4*n^2 - 25)^4). (End)
For odd n, Pi = (2^(n-1)/A001818((n-1)/2))*gamma(n/2)^2. - Alan Michael Gómez Calderón, Mar 11 2022
Pi = 4/phi + Sum_{n >= 0} (1/phi^(12*n)) * ( 8/((12*n+3)*phi^3) + 4/((12*n+5)*phi^5) - 4/((12*n+7)*phi^7) - 8/((12*n+9)*phi^9) - 4/((12*n+11)*phi^11) + 4/((12*n+13)*phi^13) ) where phi = (1+sqrt(5))/2. - Chittaranjan Pardeshi, May 16 2022
Pi = sqrt(3)*(27*S - 36)/24, where S = A248682. - Peter Luschny, Jul 22 2022
Equals Integral_{x=0..1} 1/sqrt(x-x^2) dx. - Michal Paulovic, Sep 24 2023
From Peter Bala, Oct 28 2023: (Start)
Pi = 48*Sum_{n >= 0} (-1)^n/((6*n + 1)*(6*n + 3)*(6*n + 5)).
More generally, it appears that for k >= 0 we have Pi = A(k) + B(k)*Sum_{n >= 0} (-1)^n/((6*n + 1)*(6*n + 3)*...*(6*n + 6*k + 5)), where A(k) is a rational approximation to Pi and B(k) = (3 * 2^(3*k+3) * (3*k + 2)!) / (2^(3*k+1) - (-1)^k). The first few values of A(k) for k >= 0 are [0, 256/85, 65536/20955, 821559296/261636375, 6308233216/2008080987, 908209489444864/289093830828075, ...].
Pi = 16/5 - (288/5)*Sum_{n >= 0} (-1)^n * (6*n + 1)/((6*n + 1)*(6*n + 3)*...*(6*n + 9)).
More generally, it appears that for k >= 0 we have Pi = C(k) + D(k)*Sum_{n >= 0} (-1)^n* (6*n + 1)/((6*n + 1)*(6*n + 3)*...*(6*n + 6*k + 3)), where C(k) and D(k) are rational numbers. The case k = 0 is the Madhava-Gregory-Leibniz series for Pi.
Pi = 168/53 + (288/53)*Sum_{n >= 0} (-1)^n * (42*n^2 + 25*n)/((6*n + 1)*(6*n + 3)*(6*n + 5)*(6*n + 7)).
More generally, it appears that for k >= 1 we have Pi = E(k) + F(k)*Sum_{n >= 0} (-1)^n * (6*(6*k + 1)*n^2 + (24*k + 1)*n)/((6*n + 1)*(6*n + 3)*...*(6*n + 6*k + 1)), where E(k) and F(k) are rational numbers. (End)
From Peter Bala, Nov 10 2023: (Start)
The series representation Pi = 4 * Sum_{k >= 0} 1/(4*k + 1) - 1/(4*k + 3) given above by Alexander R. Povolotsky, Dec 25 2008, is the case n = 0 of the more general result (obtained by the WZ method): for n >= 0, there holds
Pi = Sum_{j = 0.. n-1} 2^(j+1)/((2*j + 1)*binomial(2*j,j)) + 8*(n+1)!*Sum_{k >= 0} 1/((4*k + 1)*(4*k + 3)*...*(4*k + 2*n + 3)).
Letting n -> oo gives the rapidly converging series Pi = Sum_{j >= 0} 2^(j+1)/((2*j + 1)*binomial(2*j,j)) due to Euler.
More generally, it appears that for n >= 1, Pi = 1/(2*n-1)!!^2 * Sum_{j >= 0} (Product_{i = 0..2*n-1} j - i) * 2^(j+1)/((2*j + 1)*binomial(2*j,j)).
For any integer n, Pi = (-1)^n * 4 * Sum_{k >= 0} 1/(4*k + 1 + 2*n) - 1/(4*k + 3 - 2*n). (End)
Pi = Product_{k>=1} ((k^3*(k + 2)*(2*k + 1)^2)/((k + 1)^4*(2*k - 1)^2))^k. - Antonio Graciá Llorente, Jun 13 2024
Equals Integral_{x=0..2} sqrt(8 - x^2) dx - 2 (see Ambrisi and Rizzi). - Stefano Spezia, Jul 21 2024
Equals 3 + 4*Sum_{k>0} (-1)^(k+1)/(4*k*(1+k)*(1+2*k)) (see Wells at p. 53). - Stefano Spezia, Aug 31 2024
Equals 4*Integral_{x=0..1} sqrt(1 - x^2) dx = lim_{n->oo} (4/n^2)*Sum_{k=0..n} sqrt(n^2 - k^2) (see Finch). - Stefano Spezia, Oct 19 2024
Equals Beta(1/2,1/2) (see Shamos). - Stefano Spezia, Jun 03 2025
From Kritsada Moomuang, Jun 18 2025: (Start)
Equals 2 + Integral_{x=0..1} 1/(sqrt(x)*(1 + sqrt(1 - x))) dx.
Equals 2 + Integral_{x=0..1} log(1 + sqrt(1 - x))/sqrt(x) dx. (End)
Pi = 2*arccos(1/phi) + arccos(1/phi^3) = 4*arcsin(1/phi) + 2*arcsin(1/phi^3) where phi = (1+sqrt(5))/2. - Chittaranjan Pardeshi, Jul 02 2025
Pi = Sum_{n >= 0} zeta(2*n)*(2^(2*n - 1) - 1)/2^(4*n - 3). - Andrea Pinos, Jul 29 2025

Extensions

Additional comments from William Rex Marshall, Apr 20 2001

A002415 4-dimensional pyramidal numbers: a(n) = n^2*(n^2-1)/12.

Original entry on oeis.org

0, 0, 1, 6, 20, 50, 105, 196, 336, 540, 825, 1210, 1716, 2366, 3185, 4200, 5440, 6936, 8721, 10830, 13300, 16170, 19481, 23276, 27600, 32500, 38025, 44226, 51156, 58870, 67425, 76880, 87296, 98736, 111265, 124950, 139860, 156066, 173641, 192660, 213200, 235340
Offset: 0

Views

Author

Keywords

Comments

Also number of ways to legally insert two pairs of parentheses into a string of m := n-1 letters. (There are initially 2C(m+4,4) (A034827) ways to insert the parentheses, but we must subtract 2(m+1) for illegal clumps of 4 parentheses, 2m(m+1) for clumps of 3 parentheses, C(m+1,2) for 2 clumps of 2 parentheses and (m-1)C(m+1,2) for 1 clump of 2 parentheses, giving m(m+1)^2(m+2)/12 = n^2*(n^2-1)/12.) See also A000217.
E.g., for n=2 there are 6 ways: ((a))b, ((a)b), ((ab)), (a)(b), (a(b)), a((b)).
Let M_n denote the n X n matrix M_n(i,j)=(i+j); then the characteristic polynomial of M_n is x^(n-2) * (x^2-A002378(n)*x - a(n)). - Benoit Cloitre, Nov 09 2002
Let M_n denote the n X n matrix M_n(i,j)=(i-j); then the characteristic polynomial of M_n is x^n + a(n)x^(n-2). - Michael Somos, Nov 14 2002 [See A114327 for the infinite matrix M in triangular form. - Wolfdieter Lang, Feb 05 2018]
Number of permutations of [n] which avoid the pattern 132 and have exactly 2 descents. - Mike Zabrocki, Aug 26 2004
Number of tilings of a <2,n,2> hexagon.
a(n) is the number of squares of side length at least 1 having vertices at the points of an n X n unit grid of points (the vertices of an n-1 X n-1 chessboard). [For a proof, see Comments in A051602. - N. J. A. Sloane, Sep 29 2021] For example, on the 3 X 3 grid (the vertices of a 2 X 2 chessboard) there are four 1 X 1 squares, one (skew) sqrt(2) X sqrt(2) square, and one 3 X 3 square, so a(3)=6. On the 4 X 4 grid (the vertices of a 3 X 3 chessboard) there are 9 1 X 1 squares, 4 2 X 2 squares, 1 3 X 3 square, 4 sqrt(2) X sqrt(2) squares, and 2 sqrt(5) X sqrt(5) squares, so a(4) = 20. See also A024206, A108279. [Comment revised by N. J. A. Sloane, Feb 11 2015]
Kekulé numbers for certain benzenoids. - Emeric Deutsch, Jun 12 2005
Number of distinct components of the Riemann curvature tensor. - Gene Ward Smith, Apr 24 2006
a(n) is the number of 4 X 4 matrices (symmetrical about each diagonal) M = [a,b,c,d;b,e,f,c;c,f,e,b;d,c,b,a] with a+b+c+d=b+e+f+c=n+2; (a,b,c,d,e,f natural numbers). - Philippe Deléham, Apr 11 2007
If a 2-set Y and an (n-2)-set Z are disjoint subsets of an n-set X then a(n-3) is the number of 5-subsets of X intersecting both Y and Z. - Milan Janjic, Sep 19 2007
a(n) is the number of Dyck (n+1)-paths with exactly n-1 peaks. - David Callan, Sep 20 2007
Starting (1,6,20,50,...) = third partial sums of binomial transform of [1,2,0,0,0,...]. a(n) = Sum_{i=0..n} C(n+3,i+3)*b(i), where b(i)=[1,2,0,0,0,...]. - Borislav St. Borisov (b.st.borisov(AT)abv.bg), Mar 05 2009
4-dimensional square numbers. - Borislav St. Borisov (b.st.borisov(AT)abv.bg), Mar 05 2009
Equals row sums of triangle A177877; a(n), n > 1 = (n-1) terms in (1,2,3,...) dot (...,3,2,1) with additive carryovers. Example: a(4) = 20 = (1,2,3) dot (3,2,1) with carryovers = (1*3) + (2*2 + 3) + (3*1 + 7) = (3 + 7 + 10).
Convolution of the triangular numbers A000217 with the odd numbers A004273.
a(n+2) is the number of 4-tuples (w,x,y,z) with all terms in {0,...,n} and w-x=max{w,x,y,z}-min{w,x,y,z}. - Clark Kimberling, May 28 2012
The second level of finite differences is a(n+2) - 2*a(n+1) + a(n) = (n+1)^2, the squares. - J. M. Bergot, May 29 2012
Because the differences of this sequence give A000330, this is also the number of squares in an n+1 X n+1 grid whose sides are not parallel to the axes.
a(n+2) gives the number of 2*2 arrays that can be populated with 0..n such that rows and columns are nondecreasing. - Jon Perry, Mar 30 2013
For n consecutive numbers 1,2,3,...,n, the sum of all ways of adding the k-tuples of consecutive numbers for n=a(n+1). As an example, let n=4: (1)+(2)+(3)+(4)=10; (1+2)+(2+3)+(3+4)=15; (1+2+3)+(2+3+4)=15; (1+2+3+4)=10 and the sum of these is 50=a(4+1)=a(5). - J. M. Bergot, Apr 19 2013
If P(n,k) = n*(n+1)*(k*n-k+3)/6 is the n-th (k+2)-gonal pyramidal number, then a(n) = P(n,k)*P(n-1,k-1) - P(n-1,k)*P(n,k-1). - Bruno Berselli, Feb 18 2014
For n > 1, a(n) = 1/6 of the area of the trapezoid created by the points (n,n+1), (n+1,n), (1,n^2+n), (n^2+n,1). - J. M. Bergot, May 14 2014
For n > 3, a(n) is twice the area of a triangle with vertices at points (C(n,4),C(n+1,4)), (C(n+1,4),C(n+2,4)), and (C(n+2,4),C(n+3,4)). - J. M. Bergot, Jun 03 2014
a(n) is the dimension of the space of metric curvature tensors (those having the symmetries of the Riemann curvature tensor of a metric) on an n-dimensional real vector space. - Daniel J. F. Fox, Dec 15 2018
Coefficients in the terminating series identity 1 - 6*n/(n + 5) + 20*n*(n - 1)/((n + 5)*(n + 6)) - 50*n*(n - 1)*(n - 2)/((n + 5)*(n + 6)*(n + 7)) + ... = 0 for n = 1,2,3,.... Cf. A000330 and A005585. - Peter Bala, Feb 18 2019

Examples

			a(7) = 6*21 - (6*0 + 4*1 + 2*3 + 0*6 - 2*10 - 4*15) = 196. - _Bruno Berselli_, Jun 22 2013
G.f. = x^2 + 6*x^3 + 20*x^4 + 50*x^5 + 105*x^6 + 196*x^7 + 336*x^8 + ...
		

References

  • O. D. Anderson, Find the next sequence, J. Rec. Math., 8 (No. 4, 1975-1976), 241.
  • A. H. Beiler, Recreations in the Theory of Numbers, Dover, NY, 1964, p. 195.
  • S. J. Cyvin and I. Gutman, Kekulé structures in benzenoid hydrocarbons, Lecture Notes in Chemistry, No. 46, Springer, New York, 1988 (p.165).
  • R. Euler and J. Sadek, "The Number of Squares on a Geoboard", Journal of Recreational Mathematics, 251-5 30(4) 1999-2000 Baywood Pub. NY
  • S. Mukai, An Introduction to Invariants and Moduli, Cambridge, 2003; see p. 238.
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

a(n) = ((-1)^n)*A053120(2*n, 4)/8 (one-eighth of fifth unsigned column of Chebyshev T-triangle, zeros omitted). Cf. A001296.
Second row of array A103905.
Third column of Narayana numbers A001263.
Partial sums of A000330.
The expression binomial(m+n-1,n)^2-binomial(m+n,n+1)*binomial(m+n-2,n-1) for the values m = 2 through 14 produces sequences A000012, A000217, A002415, A006542, A006857, A108679, A134288, A134289, A134290, A134291, A140925, A140935, A169937.
Cf. A220212 for a list of sequences produced by the convolution of the natural numbers (A000027) with the k-gonal numbers.

Programs

  • GAP
    List([0..45],n->Binomial(n^2,2)/6); # Muniru A Asiru, Dec 15 2018
  • Magma
    [n^2*(n^2-1)/12: n in [0..50]]; // Wesley Ivan Hurt, May 14 2014
    
  • Maple
    A002415 := proc(n) binomial(n^2,2)/6 ; end proc: # Zerinvary Lajos, Jan 07 2008
  • Mathematica
    Table[(n^4 - n^2)/12, {n, 0, 40}] (* Zerinvary Lajos, Mar 21 2007 *)
    LinearRecurrence[{5,-10,10,-5,1},{0,0,1,6,20},40] (* Harvey P. Dale, Nov 29 2011 *)
  • PARI
    a(n) = n^2 * (n^2 - 1) / 12;
    
  • PARI
    x='x+O('x^200); concat([0, 0], Vec(x^2*(1+x)/(1-x)^5)) \\ Altug Alkan, Mar 23 2016
    

Formula

G.f.: x^2*(1+x)/(1-x)^5. - Simon Plouffe in his 1992 dissertation
a(n) = Sum_{i=0..n} (n-i)*i^2 = a(n-1) + A000330(n-1) = A000217(n)*A000292(n-2)/n = A000217(n)*A000217(n-1)/3 = A006011(n-1)/3, convolution of the natural numbers with the squares. - Henry Bottomley, Oct 19 2000
a(n)+1 = A079034(n). - Mario Catalani (mario.catalani(AT)unito.it), Feb 12 2003
a(n) = 2*C(n+2, 4) - C(n+1, 3). - Paul Barry, Mar 04 2003
a(n) = C(n+2, 4) + C(n+1, 4). - Paul Barry, Mar 13 2003
a(n) = Sum_{k=1..n} A000330(n-1). - Benoit Cloitre, Jun 15 2003
a(n) = n*C(n+1,3)/2 = C(n+1,3)*C(n+1,2)/(n+1). - Mitch Harris, Jul 06 2006
a(n) = A006011(n)/3 = A008911(n)/2 = A047928(n-1)/12 = A083374(n)/6. - Zerinvary Lajos, May 09 2007
a(n) = (1/2)*Sum_{1 <= x_1, x_2 <= n} (det V(x_1,x_2))^2 = (1/2)*Sum_{1 <= i,j <= n} (i-j)^2, where V(x_1,x_2) is the Vandermonde matrix of order 2. - Peter Bala, Sep 21 2007
a(n) = C(n+1,3) + 2*C(n+1,4). - Borislav St. Borisov (b.st.borisov(AT)abv.bg), Mar 05 2009
a(n) = (1/48)*sinh(2*arccosh(n))^2. - Artur Jasinski, Feb 10 2010
a(n) = n*A000292(n-1)/2. - Tom Copeland, Sep 13 2011
a(n) = 5*a(n-1) - 10*a(n-2) + 10*a(n-3) - 5*a(n-4) + a(n-5), n > 4. - Harvey P. Dale, Nov 29 2011
a(n) = (n-1)*A000217(n-1) - Sum_{i=0..n-2} (n-1-2*i)*A000217(i) for n > 1. - Bruno Berselli, Jun 22 2013
a(n) = C(n,2)*C(n+1,3) - C(n,3)*C(n+1,2). - J. M. Bergot, Sep 17 2013
a(n) = Sum_{k=1..n} ( (2k-n)* k(k+1)/2 ). - Wesley Ivan Hurt, Sep 26 2013
a(n) = floor(n^2/3) + 3*Sum_{k=1..n} k^2*floor((n-k+1)/3). - Mircea Merca, Feb 06 2014
Euler transform of length 2 sequence [6, -1]. - Michael Somos, May 28 2014
G.f. x^2*2F1(3,4;2;x). - R. J. Mathar, Aug 09 2015
Sum_{n>=2} 1/a(n) = 21 - 2*Pi^2 = 1.260791197821282762331... . - Vaclav Kotesovec, Apr 27 2016
a(n) = A080852(2,n-2). - R. J. Mathar, Jul 28 2016
a(n) = A046092(n) * A046092(n-1)/48 = A000217(n) * A000217(n-1)/3. - Bruce J. Nicholson, Jun 06 2017
E.g.f.: (1/12)*exp(x)*x^2*(6 + 6*x + x^2). - Stefano Spezia, Dec 07 2018
Sum_{n>=2} (-1)^n/a(n) = Pi^2 - 9 (See A002388). - Amiram Eldar, Jun 28 2020

Extensions

Typo in link fixed by Matthew Vandermast, Nov 22 2010
Redundant comment deleted and more detail on relationship with A000330 added by Joshua Zucker, Jan 01 2013

A002161 Decimal expansion of square root of Pi.

Original entry on oeis.org

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

Views

Author

Keywords

Comments

Also Gamma(1/2). - Franklin T. Adams-Watters, Apr 07 2006
The integral of the Gaussian function exp(-x^2) over the real line. - Richard Chapling (r.chappers(AT)gmail.com), Jun 05 2008
Also equals the average distance between two points in two dimensions where coordinates are independent normally distributed random variables with mean 0 and variance 1. - Jean-François Alcover, Oct 31 2014, after Steven Finch
Also diameter of a sphere whose surface area equals Pi^2. More generally, the square root of x is also the diameter of a sphere whose surface area equals x*Pi. - Omar E. Pol, Nov 11 2018
Convergents of continued fractions: 7/4, 16/9, 23/13, 39/22, 257/145, 296/167, 8545/4821, ... - R. J. Mathar, Jan 29 2025

Examples

			1.7724538509055160272981674833411451827975494561223871282138...
		

References

  • George Boros and Victor H. Moll, Irresistible integrals, Cambridge University Press (2006), p. 190.
  • Steven R. Finch, Mathematical Constants, Encyclopedia of Mathematics and its Applications, vol. 94, Cambridge University Press, 2003, Section 1.5.4, p. 33.
  • W. E. Mansell, Tables of Natural and Common Logarithms. Royal Society Mathematical Tables, Vol. 8, Cambridge Univ. Press, 1964, p. XVIII.
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
  • Jerome Spanier and Keith B. Oldham, "Atlas of Functions", Hemisphere Publishing Corp., 1987, chapter 43, page 413.
  • David Wells, The Penguin Dictionary of Curious and Interesting Numbers. Penguin Books, NY, 1986, Revised edition 1987. See p. 40.

Crossrefs

Cf. decimal expansions of Gamma(1/k): A073005 (k=3), A068466 (k=4), A175380 (k=5), A175379 (k=6), A220086 (k=7), A203142 (k=8).

Programs

  • Magma
    R:= RealField(100); Sqrt(Pi(R));  // G. C. Greubel, Mar 10 2018
  • Maple
    evalf(sqrt(Pi),120); # Muniru A Asiru, Nov 11 2018
  • Mathematica
    RealDigits[N[Sqrt[Pi], 120]][[1]] (* Richard Chapling (r.chappers(AT)gmail.com), Jun 05 2008 *)
  • PARI
    default(realprecision, 20080); x=sqrt(Pi); for (n=1, 20000, d=floor(x); x=(x-d)*10; write("b002161.txt", n, " ", d)); \\ Harry J. Smith, May 01 2009
    

Formula

Equals (1/2) * Sum_{n>=0} ((-1)^n * (4*n+1) * (1/8)^(n+1) * (2^(n+1))^3 * Gamma(n+1/2)^3 / Gamma(n+1)^3). - Alexander R. Povolotsky, Mar 25 2013
Equals Integral_{x=0..1} 1/sqrt(-log(x)) dx. - Jean-François Alcover, Apr 29 2013
Equals Sum_{k>=0} (k+1/2)!/(k+2)!. - Amiram Eldar, Jun 19 2023
Equals Integral_{x=0..oo} exp(-x)/sqrt(x) dx. - Michal Paulovic, Sep 24 2023
Equals Integral_{x=0..oo} 4/(exp(x^2)*(2*x^2 + 1)^2) dx. - Kritsada Moomuang, Jun 05 2025

Extensions

More terms from Franklin T. Adams-Watters, Apr 07 2006

A019670 Decimal expansion of Pi/3.

Original entry on oeis.org

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

Views

Author

N. J. A. Sloane, Dec 11 1996

Keywords

Comments

With an offset of zero, also the decimal expansion of Pi/30 ~ 0.104719... which is the average arithmetic area of the 0-winding sectors enclosed by a closed Brownian planar path, of a given length t, according to Desbois, p. 1. - Jonathan Vos Post, Jan 23 2011
Polar angle (or apex angle) of the cone that subtends exactly one quarter of the full solid angle. See comments in A238238. - Stanislav Sykora, Jun 07 2014
60 degrees in radians. - M. F. Hasler, Jul 08 2016
Volume of a quarter sphere of radius 1. - Omar E. Pol, Aug 17 2019
Also smallest positive zero of Sum_{k>=1} cos(k*x)/k = -log(2*|sin(x/2)|). Proof of this identity: Sum_{k>=1} cos(k*x)/k = Re(Sum_{k>=1} exp(k*x*i)/k) = Re(-log(1-exp(x*i))) = -log(2*|sin(x/2)|), x != 2*m*Pi, where i = sqrt(-1). - Jianing Song, Nov 09 2019
The area of a circle circumscribing a unit-area regular dodecagon. - Amiram Eldar, Nov 05 2020

Examples

			Pi/3 = 1.04719755119659774615421446109316762806572313312503527365831486...
From _Peter Bala_, Nov 16 2016: (Start)
Case n = 1. Pi/3 = 18 * Sum_{k >= 0} (-1)^(k+1)( 1/((6*k - 5)*(6*k + 1)*(6*k + 7)) + 1/((6*k - 1)*(6*k + 5)*(6*k + 11)) ).
Using the methods of Borwein et al. we can find the following asymptotic expansion for the tails of this series: for N divisible by 6 there holds Sum_{k >= N/6} (-1)^(k+1)( 1/((6*k - 5)*(6*k + 1)*(6*k + 7)) + 1/((6*k - 1)*(6*k + 5)*(6*k + 11)) ) ~ 1/N^3 + 6/N^5 + 1671/N ^7 - 241604/N^9 + ..., where the sequence [1, 0, 6, 0, 1671, 0, -241604, 0, ...] is the sequence of coefficients in the expansion of ((1/18)*cosh(2*x)/cosh(3*x)) * sinh(3*x)^2 = x^2/2! + 6*x^4/4! + 1671*x^6/6! - 241604*x^8/8! + .... Cf. A024235, A278080 and A278195. (End)
		

References

  • Steven R. Finch, Mathematical Constants, Encyclopedia of Mathematics and its Applications, vol. 94, Cambridge University Press, 2003, Section 8.3, p. 489.

Crossrefs

Integral_{x=0..oo} 1/(1+x^m) dx: A013661 (m=2), A248897 (m=3), A093954 (m=4), A352324 (m=5), this sequence (m=6), A352125 (m=8), A094888 (m=10).

Programs

Formula

A third of A000796, a sixth of A019692, the square root of A100044.
Sum_{k >= 0} (-1)^k/(6k+1) + (-1)^k/(6k+5). - Charles R Greathouse IV, Sep 08 2011
Product_{k >= 1}(1-(6k)^(-2))^(-1). - Fred Daniel Kline, May 30 2013
From Peter Bala, Feb 05 2015: (Start)
Pi/3 = Sum {k >= 0} binomial(2*k,k)*1/(2*k + 1)*(1/16)^k = 2F1(1/2,1/2;3/2;1/4). Similar series expansions hold for Pi^2 (A002388), Pi^3 (A091925) and Pi/(2*sqrt(2)) (A093954.)
The integer sequences A(n) := 4^n*(2*n + 1)! and B(n) := A(n)*( Sum {k = 0..n} binomial(2*k,k)*1/(2*k + 1)*(1/16)^k ) both satisfy the second-order recurrence equation u(n) = (20*n^2 + 4*n + 1)*u(n-1) - 8*(n - 1)*(2*n - 1)^3*u(n-2). From this observation we can obtain the continued fraction expansion Pi/3 = 1 + 1/(24 - 8*3^3/(89 - 8*2*5^3/(193 - 8*3*7^3/(337 - ... - 8*(n - 1)*(2*n - 1)^3/((20*n^2 + 4*n + 1) - ... ))))). Cf. A002388 and A093954. (End)
Equals Sum_{k >= 1} arctan(sqrt(3)*L(2k)/L(4k)) where L=A000032. See also A005248 and A056854. - Michel Marcus, Mar 29 2016
Equals Product_{n >= 1} A016910(n) / A136017(n). - Fred Daniel Kline, Jun 09 2016
Equals Integral_{x=-oo..oo} sech(x)/3 dx. - Ilya Gutkovskiy, Jun 09 2016
From Peter Bala, Nov 16 2016: (Start)
Euler's series transformation applied to the series representation Pi/3 = Sum_{k >= 0} (-1)^k/(6*k + 1) + (-1)^k/(6*k + 5) given above by Greathouse produces the faster converging series Pi/3 = (1/2) * Sum_{n >= 0} 3^n*n!*( 1/(Product_{k = 0..n} (6*k + 1)) + 1/(Product_{k = 0..n} (6*k + 5)) ).
The series given above by Greathouse is the case n = 0 of the more general result Pi/3 = 9^n*(2*n)! * Sum_{k >= 0} (-1)^(k+n)*( 1/(Product_{j = -n..n} (6*k + 1 + 6*j)) + 1/(Product_{j = -n..n} (6*k + 5 + 6*j)) ) for n = 0,1,2,.... Cf. A003881. See the example section for notes on the case n = 1.(End)
Equals Product_{p>=5, p prime} p/sqrt(p^2-1). - Dimitris Valianatos, May 13 2017
Equals A019699/4 or A019693/2. - Omar E. Pol, Aug 17 2019
Equals Integral_{x >= 0} (sin(x)/x)^4 = 1/2 + Sum_{n >= 0} (sin(n)/n)^4, by the Abel-Plana formula. - Peter Bala, Nov 05 2019
Equals Integral_{x=0..oo} 1/(1 + x^6) dx. - Bernard Schott, Mar 12 2022
Pi/3 = -Sum_{n >= 1} i/(n*P(n, 1/sqrt(-3))*P(n-1, 1/sqrt(-3))), where i = sqrt(-1) and P(n, x) denotes the n-th Legendre polynomial. The first twenty terms of the series gives the approximation Pi/3 = 1.04719755(06...) correct to 8 decimal places. - Peter Bala, Mar 16 2024
Equals Integral_{x >= 0} (2*x^2 + 1)/((x^2 + 1)*(4*x^2 + 1)) dx. - Peter Bala, Feb 12 2025

A111003 Decimal expansion of Pi^2/8.

Original entry on oeis.org

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

Views

Author

Sam Alexander, Oct 01 2005

Keywords

Comments

According to Beckmann, Euler discovered the formula for this number as a sum of squares of reciprocals of odd numbers, along with similar formulas for Pi^2/6 and Pi^2/12. - Alonso del Arte, Apr 01 2013
Equals the asymptotic mean of the abundancy index of the odd numbers. - Amiram Eldar, May 12 2023

Examples

			1.23370055013616982735431137498451889191421242590509882830166867202...
1 + 1/9 + 1/25 + 1/49 + 1/81 + 1/121 + 1/169 + 1/225 + ... - _Bruno Berselli_, Mar 06 2017
		

References

  • F. Aubonnet, D. Guinin and B. Joppin, Précis de Mathématiques, Analyse 2, Classes Préparatoires, Premier Cycle Universitaire, Bréal, 1990, Exercice 908, pages 82 and 91-92.
  • Petr Beckmann, A History of Pi, 5th Ed. Boulder, Colorado: The Golem Press (1982): p. 153.
  • George Boros and Victor H. Moll, Irresistible integrals, Cambridge University Press (2006), p. 122.
  • Calvin C. Clawson, The Beauty and Magic of Numbers. New York: Plenum Press (1996): 98.
  • L. B. W. Jolley, Summation of Series, Dover (1961).
  • David Wells, The Penguin Dictionary of Curious and Interesting Numbers. Penguin Books, NY, 1986, Revised edition 1987. See p. 54.

Crossrefs

Programs

Formula

Equals 1 + 1/(2*3) + (1/3)*(1*2)/(3*5) + (1/4)*(1*2*3)/(3*5*7) + ... [Jolley eq 276]
Equals Sum_{k >= 1} 1/(2*k - 1)^2 [Clawson and Wells]. - Alonso del Arte, Aug 15 2012
Equals 2*(Integral_{t=0..1} sqrt(1 - t^2) dt)^2. - Alonso del Arte, Mar 29 2013
Equals Sum_{k >= 1} 2^k/(k^2*binomial(2*k, k)). - Jean-François Alcover, Apr 29 2013
Equals Integral_{x=0..1} log((1+x^2)/(1-x^2))/x dx. - Bruno Berselli, May 13 2013
Equals limit_{p->0} Integral_{x=0..Pi/2} x*tan(x)^p dx. [Jean-François Alcover, May 17 2013, after Boros & Moll p. 230]
Equals A002388/8 = A102753/4 = A091476/2. - R. J. Mathar, Oct 13 2015
Equals Integral_{x>=0} x*K_0(x)*K_1(x)dx where K are modified Bessel functions [Gradsteyn-Ryzhik 6.576.4]. - R. J. Mathar, Oct 22 2015
Equals (3/4)*zeta(2) = (3/4)*A013661. - Wolfdieter Lang, Sep 02 2019
From Amiram Eldar, Jul 17 2020: (Start)
Equals -Integral_{x=0..1} log(x)/(1 - x^2) dx = Integral_{x>=1} log(x)/(x^2-1) dx.
Equals -Integral_{x=0..oo} log(x)/(1 - x^4) dx.
Equals Integral_{x=0..oo} arctan(x)/(1 + x^2) dx. (End)
Equals Integral_{x=0..1} log(1+x+x^2+x^3)/x dx (Aubonnet). - Bernard Schott, Feb 04 2022
Equals Sum_{n>=1} A309891(n)/n^2. - Friedjof Tellkamp, Jan 25 2025
Equals lambda(2), where lambda is the Dirichlet lambda function. - Michel Marcus, Aug 15 2025

Extensions

More terms from Robert G. Wilson v, Oct 04 2005

A091925 Decimal expansion of Pi^3.

Original entry on oeis.org

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

Views

Author

Mohammad K. Azarian, Mar 16 2004

Keywords

Comments

Surface area of the 6-dimensional unit sphere. - Stanislav Sykora, Nov 08 2013
Surface area of a sphere of diameter Pi equals the volume of the circumscribed cube. - Omar E. Pol, Dec 25 2013
Area of a circle of radius Pi. - Omar E. Pol, Jan 31 2016

Examples

			31.00627668029982017547631506710139520222528856588510769414453810380639...
		

Crossrefs

Cf. A000796, A002388, A058285 (continued fraction), A019670, A093954, A092731, A092735.

Programs

  • Magma
    R:= RealField(100); (Pi(R))^3; // G. C. Greubel, Mar 09 2018
  • Mathematica
    First@ RealDigits@ N[Pi^3, 120] (* Michael De Vlieger, Jan 31 2016 *)
  • PARI
    default(realprecision, 20080); x=Pi^3/10; for (n=2, 20000, d=floor(x); x=(x-d)*10; write("b091925.txt", n, " ", d)); \\ Harry J. Smith, Jun 22 2009
    

Formula

Sum_{k >= 0} binomial(2*k,k)/((2*k + 1)^3*16^k) = 7*Pi^3/216. (Kh. Hessami Pilehrood and T. Hessami Pilehrood).
From Peter Bala, Feb 05 2015: (Start)
The integer sequences A(n) := 2^n*(2*n + 1)!^3/n!^2 and B(n) := A(n)*( Sum {k = 0..n} binomial(2*k,k)*1/(2*k + 1)^3*(1/16)^k ) both satisfy the second order recurrence equation u(n) = (160*n^4 + 128*n^3 + 144*n^2 + 2)*u(n-1) - 32*(n - 1)*(2*n - 1)^7*u(n-2). From this observation we can obtain the continued fraction expansion 7/216*Pi^3 = 1 + 2/(432 - 32*3^7/(4162 - 32*2*5^7/(17714 - ... - 32*(n - 1)*(2*n - 1)^7/((160*n^4 + 128*n^3 + 144*n^2 + 2) - ... )))). Cf. A002388, A019670 and A093954. (End)
From Peter Bala, Oct 31 2019: (Start)
Pi^3 = (1/7) * Sum_{n >= 0} (-1)^n*( 1/(n + 1/6)^3 + 1/(n + 5/6)^3 ).
Pi^3 = (1/31) * Sum_{n >= 0} (-1)^n*( 1/(n + 1/10)^3 - 1/(n + 3/10)^3 - 1/(n + 7/10)^3 + 1/(n + 9/10)^3 ). Cf. A019692, A092731 and A092735. (End)
Equals Integral_{x=-oo..oo} x^2/(exp(x/2) + exp(-x/2)) dx. - Amiram Eldar, May 21 2021

A093954 Decimal expansion of Pi/(2*sqrt(2)).

Original entry on oeis.org

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

Views

Author

Eric W. Weisstein, Apr 19 2004

Keywords

Comments

The value is the length Pi*sqrt(2)/4 of the diagonal in the square with side length Pi/4 = Sum_{n>=0} (-1)^n/(2n+1) = A003881. The area of the circumcircle of this square is Pi*(Pi*sqrt(2)/8)^2 = Pi^3/32 = A153071. - Eric Desbiaux, Jan 18 2009
This is the value of the Dirichlet L-function of modulus m=8 at argument s=1 for the non-principal character (1,0,1,0,-1,0,-1,0). See arXiv:1008.2547. - R. J. Mathar, Mar 22 2011
Archimedes's-like scheme: set p(0) = sqrt(2), q(0) = 1; p(n+1) = 2*p(n)*q(n)/(p(n)+q(n)) (harmonic mean, i.e., 1/p(n+1) = (1/p(n) + 1/q(n))/2), q(n+1) = sqrt(p(n+1)*q(n)) (geometric mean, i.e., log(q(n+1)) = (log(p(n+1)) + log(q(n)))/2), for n >= 0. The error of p(n) and q(n) decreases by a factor of approximately 4 each iteration, i.e., approximately 2 bits are gained by each iteration. Set r(n) = (2*q(n) + p(n))/3, the error decreases by a factor of approximately 16 for each iteration, i.e., approximately 4 bits are gained by each iteration. For a similar scheme see also A244644. - A.H.M. Smeets, Jul 12 2018
The area of a circle circumscribing a unit-area regular octagon. - Amiram Eldar, Nov 05 2020

Examples

			1.11072073453959156175397...
From _Peter Bala_, Mar 03 2015: (Start)
Asymptotic expansion at n = 5000.
The truncated series Sum_{k = 0..5000 - 1} (-1)^floor(k/2)/(2*k + 1) = 1.110(6)207345(42)591561(18)3970(5238)1.... The bracketed digits show where this decimal expansion differs from that of Pi/(2*sqrt(2)). The numbers 1, -3, 57, -2763 must be added to the bracketed numbers to give the correct decimal expansion to 30 digits: Pi/(2*sqrt(2)) = 1.110(7)207345(39)591561(75)3970 (2475)1.... (End)
From _Peter Bala_, Nov 24 2016: (Start)
Case m = 1, n = 1:
Pi/(2*sqrt(2)) = 4*Sum_{k >= 0} (-1)^(1 + floor(k/2))/((2*k - 1)*(2*k + 1)*(2*k + 3)).
We appear to have the following asymptotic expansion for the tails of this series: for N divisible by 4, Sum_{k >= N/2} (-1)^floor(k/2)/((2*k - 1)*(2*k + 1)*(2*k + 3)) ~ 1/N^3 - 14/N^5 + 691/N^7 - 62684/N^9 - ..., where the coefficient sequence [1, 0, -14, 0, 691, 0, -62684, ...] appears to come from the e.g.f. (1/2!)*cosh(x)/cosh(2*x)*sinh(x)^2 = x^2/2! - 14*x^4/4! + 691*x^6/6! - 62684*x^8/8! + .... Cf. A019670.
For example, take N = 10^5. The truncated series Sum_{k = 0..N/2 -1} (-1)^(1+floor(k/2))/((2*k - 1)*(2*k + 1)*(2*k + 3)) = 0.27768018363489(8)89043849(11)61878(80026)6163(351171)58.... The bracketed digits show where this decimal expansion differs from that of (1/4)*Pi/(2*sqrt(2)). The numbers -1, 14, -691, 62684 must be added to the bracketed numbers to give the correct decimal expansion: (1/4)*Pi/(2*sqrt(2)) = 0.27768018363489(7) 89043849(25)61878(79335)6163(413855)58... (End)
		

References

  • J. M. Arnaudiès, P. Delezoide et H. Fraysse, Exercices résolus d'Analyse du cours de mathématiques - 2, Dunod, 1993, Exercice 5, p. 240.
  • George Boros and Victor H. Moll, Irresistible integrals, Cambridge University Press, 2006, p. 149.
  • Steven R. Finch, Mathematical Constants, Encyclopedia of Mathematics and its Applications, vol. 94, Cambridge University Press, 2003, Section 1.4.1, p. 20.
  • L. B. W. Jolley, Summation of Series, Dover, 1961, eq. 76, page 16.
  • Joel L. Schiff, The Laplace Transform: Theory and Applications, Springer-Verlag New York, Inc. (1999). See p. 149.
  • David Wells, The Penguin Dictionary of Curious and Interesting Numbers. Penguin Books, NY, 1986, Revised edition 1987. See p. 53.

Crossrefs

Programs

  • Maple
    simplify( sum((cos((1/2)*k*Pi)+sin((1/2)*k*Pi))/(2*k+1), k = 0 .. infinity) );  # Peter Bala, Mar 09 2015
  • Mathematica
    RealDigits[Pi/Sqrt@8, 10, 111][[1]] (* Michael De Vlieger, Sep 23 2016 and slightly modified by Robert G. Wilson v, Jul 23 2018 *)
  • PARI
    default(realprecision, 20080); x=Pi*sqrt(2)/4; for (n=1, 20000, d=floor(x); x=(x-d)*10; write("b093954.txt", n, " ", d)); \\ Harry J. Smith, Jun 17 2009

Formula

Equals 1/A112628.
Equals Integral_{x=0..oo} 1/(x^4+1) dx. - Jean-François Alcover, Apr 29 2013
From Peter Bala, Feb 05 2015: (Start)
Pi/(2*sqrt(2)) = Sum_{k >= 0} binomial(2*k,k)*1/(2*k + 1)*(1/8)^k.
The integer sequences A(n) := 2^n*(2*n + 1)! and B(n) := A(n)*( Sum {k = 0..n} binomial(2*k,k)*1/(2*k + 1)*(1/8)^k ) both satisfy the second order recurrence equation u(n) = (12*n^2 + 1)*u(n-1) - 4*(n - 1)*(2*n - 1)^3*u(n-2). From this observation we can obtain the continued fraction expansion Pi/(2*sqrt(2)) = 1 + 1/(12 - 4*3^3/(49 - 4*2*5^3/(109 - 4*3*7^3/(193 - ... - 4*(n - 1)*(2*n - 1)^3/((12*n^2 + 1) - ... ))))). Cf. A002388 and A019670. (End)
From Peter Bala, Mar 03 2015: (Start)
Pi/(2*sqrt(2)) = Sum_{k >= 0} (-1)^floor(k/2)/(2*k + 1) = limit (n -> infinity) Sum_{k = -n .. n - 1} (-1)^k/(4*k + 1). See Wells.
We conjecture the asymptotic expansion Pi/(2*sqrt(2)) - Sum {k = 0..n - 1} (-1)^floor(k/2)/(2*k + 1) ~ 1/(2*n) - 3/(2*n)^3 + 57/(2*n)^5 - 2763/(2*n)^7 + ..., where n is a multiple of 4 and the sequence of unsigned coefficients [1, 3, 57, 2763, ...] is A000281. An example with n = 5000 is given below. (End)
From Peter Bala, Sep 21 2016: (Start)
c = 2 * Sum_{k >= 0} (-1)^k * (4*k + 2)/((4*k + 1)*(4*k + 3)) = A181048 + A181049. The asymptotic expansion conjectured above follows from the asymptotic expansions given in A181048 and A181049.
c = 1/2 * Integral_{x = 0..Pi/2} sqrt(tan(x)) dx. (End)
From Peter Bala, Nov 24 2016: (Start)
Let m be an odd integer and n a nonnegative integer. Then Pi/(2*sqrt(2)) = 2^n*m^(2*n)*(2*n)!*Sum_{k >= 0} (-1)^(n+floor(k/2)) * 1/Product_{j = -n..n} (2*k + 1 + 2*m*j). Cf. A003881.
In the particular case m = 1 the result has the equivalent form: for n a nonnegative integer, Pi/(2*sqrt(2)) = 2^n*(2*n)!*Sum_{k >= 0} (-1)^(n+k)*(8*k + 4)* 1/Product_{j = -n..n+1} (4*k + 2*j + 1). The case m = 1, n = 1 is considered in the Example section below.
Let m be an odd integer and n a nonnegative integer. Then Pi/(2*sqrt(2)) = 4^n*m^(2*n)*(2*n)!*Sum_{k >= 0} (-1)^(n+floor(k/2)) * 1/Product_{j = -n..n} (2*k + 1 + 4*m*j). (End)
Equals Integral_{x = 0..oo} cosh(x)/cosh(2*x) dx. - Peter Bala, Nov 01 2019
Equals Sum_{k>=1} A188510(k)/k = Sum_{k>=1} Kronecker(-8,k)/k = 1 + 1/3 - 1/5 - 1/7 + 1/9 + 1/11 - 1/13 - 1/15 + ... - Jianing Song, Nov 16 2019
From Amiram Eldar, Jul 16 2020: (Start)
Equals Product_{k>=1} (1 - (-1)^k/(2*k+1)).
Equals Integral_{x=0..oo} dx/(x^2 + 2).
Equals Integral_{x=0..Pi/2} dx/(sin(x)^2 + 1). (End)
Equals Integral_{x=0..oo} x^2/(x^4 + 1) dx (Arnaudiès). - Bernard Schott, May 19 2022
Equals Integral_{x = 0..1} 1/(2*x^2 + (1 - x)^2) dx. - Peter Bala, Jul 22 2022
Equals Integral_{x = 0..1} 1/(1 - x^4)^(1/4) dx. - Terry D. Grant, Mar 17 2023
Equals 1/Product_{p prime} (1 - Kronecker(-8,p)/p), where Kronecker(-8,p) = 0 if p = 2, 1 if p == 1 or 3 (mod 8) or -1 if p == 5 or 7 (mod 8). - Amiram Eldar, Dec 17 2023
Equals A068465*A068467. - R. J. Mathar, Jun 27 2024
From Stefano Spezia, Jun 05 2025: (Start)
Equals Sum_{k>=1} (-1)^(k+1)(1/(4*k - 3) + 1/(4*k - 1)).
Equals Product_{k=0..oo} (1 + (-1)^k/(2*k + 3)).
Equals Integral_{x=0..oo} 1/(2*x^2 + 1).
Equals Integral_{x=0..1} 1/((1 + x^2)*sqrt(1 - x^2)). (End)

A092425 Decimal expansion of Pi^4.

Original entry on oeis.org

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

Views

Author

Mohammad K. Azarian, Mar 22 2004

Keywords

Examples

			97.40909103400243723644033268870511124972758567268542169146785938997085...
		

Crossrefs

Cf. A000796 (Pi), A002388 (Pi^2), A091925 (Pi^3), A092731 (Pi^5), A092732 (Pi^6), A092735 (Pi^7), A092736 (Pi^8), A058286 (continued fraction), A013662.

Programs

  • Magma
    R:= RealField(150); (Pi(R))^4; // G. C. Greubel, Mar 09 2018
    
  • Magma
    R:=RealField(110); SetDefaultRealField(R); n:=Pi(R)^4; Reverse(Intseq(Floor(10^98*n))); // Bruno Berselli, Mar 12 2018
  • Mathematica
    RealDigits[Pi^4, 10, 100][[1]] (* G. C. Greubel, Mar 09 2018 *)
  • PARI
    default(realprecision, 20080); x=Pi^4/10; for (n=2, 20000, d=floor(x); x=(x-d)*10; write("b092425.txt", n, " ", d)); \\ Harry J. Smith, Jun 22 2009
    

Formula

Equals 120 * Sum_{j>=1} Sum_{i=1..j-1} 1/(i*j)^2. - Enrique Pérez Herrero, Jun 29 2012
Equals Sum_{k>=1} k*(k+1)*(k+2)*zeta(k+3)/2^(k-1). - Amiram Eldar, May 21 2021
From Peter Bala, Oct 21 2023: (Start)
Pi^4 = 90*Sum_{n >= 1} 1/n^4 (Euler).
The following faster converging series representations for the constant Pi^4 may be easily verified using partial fraction expansions of the summands of the series. Presumably, these are the first three cases of an infinite family of similar results.
Let P(n) = n*(n + 1)*(n + 2)/2!. Then Pi^4 = 1575/16 - 15*Sum_{n >= 1} d/dn(P(n))/P(n)^4.
Let Q(n) = n*(n + 1)*(n + 2)*(n + 3)*(n + 4)/4!. Then Pi^4 = 673165/6912 + Sum_{n >= 1} d/dn(Q(n))/Q(n)^4.
Let R(n) = n*(n + 1)*(n + 2)*(n + 3)*(n + 4)*(n + 5)*(n + 6)/6!. Then Pi^4 = 5610787/57600 - (3/56)*Sum_{n >= 1} d/dn(R(n))/R(n)^4.
Taking 10 terms of the last series gives the approximation Pi^4 = 97.4090910340
024372(50...), correct to 16 decimal places. (End)

A229099 Decimal expansion of 1 - 6/Pi^2.

Original entry on oeis.org

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

Views

Author

Keywords

Comments

Probability that a random number is not squarefree; probability that two random numbers have a common divisor greater than 1.

Examples

			0.39207289814597337133672322074163416657384735196652070692634580863496...
		

Crossrefs

Programs

  • Mathematica
    RealDigits[1-1/Zeta[2],10,90][[1]] (* Stefano Spezia, Feb 21 2025 *)
  • PARI
    1-6/Pi^2

Formula

Equals 1 - 1/zeta(2). - Stefano Spezia, Feb 21 2025

A004011 Theta series of D_4 lattice; Fourier coefficients of Eisenstein series E_{gamma,2}.

Original entry on oeis.org

1, 24, 24, 96, 24, 144, 96, 192, 24, 312, 144, 288, 96, 336, 192, 576, 24, 432, 312, 480, 144, 768, 288, 576, 96, 744, 336, 960, 192, 720, 576, 768, 24, 1152, 432, 1152, 312, 912, 480, 1344, 144, 1008, 768, 1056, 288, 1872, 576, 1152, 96, 1368, 744, 1728, 336
Offset: 0

Views

Author

Keywords

Comments

D_4 is also the Barnes-Wall lattice in 4 dimensions.
E_{gamma,2} is the unique normalized modular form for Gamma_0(2) of weight 2.
Cubic AGM theta functions: a(q) (see A004016), b(q) (A005928), c(q) (A005882).
Ramanujan's Eisenstein series: P(q) (see A006352), Q(q) (A004009), R(q) (A013973).
Convolution square is A008658. - Michael Somos, Aug 21 2014
Expansion of 2*P(x^2) - P(x) in powers of x where P() is a Ramanujan Eisenstein series. - Michael Somos, Feb 16 2015
a(n) is the number of Hurwitz quaternions of norm n. - Michael Somos, Feb 16 2015

Examples

			G.f. = 1 + 24*x + 24*x^2 + 96*x^3 + 24*x^4 + 144*x^5 + 96*x^6 + 192*x^7 + 24*x^8 + ...
G.f. = 1 + 24*q^2 + 24*q^4 + 96*q^6 + 24*q^8 + 144*q^10 + 96*q^12 + 192*q^14 + 24*q^16 + ...
		

References

  • Bruce C. Berndt, Ramanujan's Notebooks Part V, Springer-Verlag, 1998, see p. 148 Eq. (9.11).
  • Harvey Cohn, Advanced Number Theory, Dover Publications, Inc., 1980, p. 89. Eq. (1).
  • J. H. Conway and N. J. A. Sloane, "Sphere Packings, Lattices and Groups", Springer-Verlag, p. 119.
  • S. Ramanujan, Notebooks, Tata Institute of Fundamental Research, Bombay 1957 Vol. 1, see page 214.
  • N. J. A. Sloane, Seven Staggering Sequences, in Homage to a Pied Puzzler, E. Pegg Jr., A. H. Schoen and T. Rodgers (editors), A. K. Peters, Wellesley, MA, 2009, pp. 93-110.
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Partial sums give A046949.
Cf. A108092 (convolution fourth root).

Programs

  • Magma
    Basis( ModularForms( Gamma0(2), 2), 54) [1]; /* Michael Somos, May 27 2014 */
    
  • Maple
    readlib(ifactors): with(numtheory): for n from 1 to 100 do if n mod 2 = 0 then m := n/ifactors(n)[2][1][1]^ifactors(n)[2][1][2] else m := n fi: printf(`%d,`,24*sigma(m)) od: # James Sellers, Dec 07 2000
  • Mathematica
    a[ n_] := If[ n < 0, 0, With[ {m = Floor @ Sqrt[4 n]}, SeriesCoefficient[ Sum[ q^( x^2 + y^2 + z^2 + t^2 + (x + y + z) t ), {x, -m, m}, {y, -m, m}, {z, -m, m}, {t, -m, m}] + O[q]^(n + 1), n]]]; (* Michael Somos, Jan 11 2011 *)
    a[n_] := 24*Total[ Select[ Divisors[n], OddQ]]; a[0]=1; Table[a[n], {n, 0, 52}] (* Jean-François Alcover, Sep 12 2012 *)
    a[ n_] := With[{m = InverseEllipticNomeQ @q}, SeriesCoefficient[ (1 + m) (EllipticK[ m] / (Pi/2))^2, {q, 0, n}]]; (* Michael Somos, Jun 04 2013 *)
    a[ n_] := With[{m = InverseEllipticNomeQ @q}, SeriesCoefficient[ (1 - m/2) (EllipticK[ m] / (Pi/2))^2, {q, 0, 2 n}]]; (* Michael Somos, Jun 04 2013 *)
    a[ n_] := SeriesCoefficient[ EllipticTheta[ 3, 0, q]^4 + EllipticTheta[ 2, 0, q]^4, {q, 0, n}]; (* Michael Somos, Jun 04 2013 *)
    a[ n_] := SeriesCoefficient[ (EllipticTheta[ 3, 0, q]^4 + EllipticTheta[ 4, 0, q]^4) / 2, {q, 0, 2 n}]; (* Michael Somos, Jun 04 2013 *)
  • PARI
    {a(n) = if( n<1, n==0, 24 * sumdiv( n, d, d%2 * d))}; /* Michael Somos, Apr 17 2000 */
    
  • PARI
    {a(n) = my(G); if( n<0, 0, G = [2, 1, 1, 1; 1, 2, 0, 0; 1, 0, 2, 0; 1, 0, 0, 2]; polcoeff( 1 + 2 * x * Ser(qfrep( G, n, 1)), n))}; /* Michael Somos, Sep 11 2007 */
    
  • PARI
    {a(n) = my(A); if( n<0, 0, A = x * O(x^n); polcoeff( eta(x^2 + A)^20 / (eta(x + A) * eta(x^4 + A))^8 + 16 * x * eta(x^4 + A)^8 / eta(x^2 + A)^4, n))}; /* Michael Somos, Oct 21 2017 */
    
  • Python
    from sympy import divisors
    def a(n): return 1 if n==0 else 24*sum(d for d in divisors(n) if d%2)
    print([a(n) for n in range(101)]) # Indranil Ghosh, Jun 24 2017
    
  • Python
    from math import prod
    from sympy import factorint
    def A004011(n): return 24*prod((p**(e+1)-1)//(p-1) for p, e in factorint(n).items() if p > 2) if n else 1 # Chai Wah Wu, Nov 13 2022
  • Sage
    ModularForms( Gamma0(2), 2, prec=54).0; # Michael Somos, Jun 04 2013
    

Formula

a(0) = 1; if n>0 then a(n) = 24 (Sum_{d|n, d odd, d>0} d) = 24 * A000593(n).
G.f.: 1 + 24 Sum_{n>0} n x^n /(1 + x^n). a(n) = A000118(2*n) = A096727(2*n).
G.f.: (1/2) * (theta_3(z)^4 + theta_4(z)^4) = theta_3(2z)^4 + theta_2(2z)^4 = Sum_{k>=0} a(k) * x^(2*k).
G.f.: Sum_{a, b, c, d in Z} x^(a^2 + b^2 + c^2 + d^2 + a*d + b*d + c*d). - Michael Somos, Jan 11 2011
Expansion of (1 + k^2) * K(k^2)^2 / (Pi/2)^2 in powers of nome q. - Michael Somos, Jun 10 2006
Expansion of (1 - k^2/2) * K(k^2)^2 / (Pi/2)^2 in powers of nome q^2. - Michael Somos, Mar 14 2012
Expansion of b(x) * b(x^2) + 3 * c(x) * c(x^2) in powers of x where b(), c() are cubic AGM theta functions. - Michael Somos, Jan 11 2011
Expansion of b(x) * b(x^2) + c(x) * c(x^2) / 3 in powers of x^3 where b(), c() are cubic AGM theta functions. - Michael Somos, Mar 14 2012
G.f. is a period 1 Fourier series which satisfies f(-1 / (2 t)) = 2 (t/i)^2 f(t) where q = exp(2 Pi i t). - Michael Somos, Sep 11 2007
G.f. A(x) satisfies 0 = f(A(x), A(x^2), A(x^4)) where f(u, v, w) = u^2 - 2*u*v - 7*v^2 - 8*v*w + 16*w^2. - Michael Somos, May 29 2005
G.f. A(x) satisfies 0 = f(A(x), A(x^2), A(x^3), A(x^6)) where f(u1, u2, u3, u6) = u1^2 + 4*u2^2 + 9*u3^2 + 36*u6^2 - 2*u1*u2 - 10*u1*u3 + 10*u1*u6 + 10*u2*u3 - 40*u2*u6 - 18*u3*u6. - Michael Somos, Sep 11 2007
Sum_{k=1..n} a(k) ~ c * n^2, where c = Pi^2 = 9.869604... (A002388). - Amiram Eldar, Dec 29 2023

Extensions

Additional comments from Barry Brent (barryb(AT)primenet.com)
Showing 1-10 of 85 results. Next