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.

Previous Showing 21-30 of 48 results. Next

A011379 a(n) = n^2*(n+1).

Original entry on oeis.org

0, 2, 12, 36, 80, 150, 252, 392, 576, 810, 1100, 1452, 1872, 2366, 2940, 3600, 4352, 5202, 6156, 7220, 8400, 9702, 11132, 12696, 14400, 16250, 18252, 20412, 22736, 25230, 27900, 30752, 33792, 37026, 40460, 44100, 47952, 52022, 56316, 60840
Offset: 0

Views

Author

Glen Burch (gburch(AT)erols.com), Felice Russo

Keywords

Comments

(1) a(n) = sum of second string of n triangular numbers - sum of first n triangular numbers, or the 2n-th partial sum of triangular numbers (A000217) - the n-th partial sum of triangular numbers (A000217). The same for natural numbers gives squares. (2) a(n) = (n-th triangular number)*(the n-th even number) = n(n+1)/2 * (2n). - Amarnath Murthy, Nov 05 2002
Let M(n) be the n X n matrix m(i,j)=1/(i+j+x), let P(n,x) = (Product_{i=0..n-1} i!^2)/det(M(n)). Then P(n,x) is a polynomial with integer coefficients of degree n^2 and a(n) is the coefficient of x^(n^2-1). - Benoit Cloitre, Jan 15 2003
Y values of solutions of the equation: (X-Y)^3-X*Y=0. X values are a(n)=n*(n+1)^2 (see A045991) - Mohamed Bouhamida, May 09 2006
a(2d-1) is the number of self-avoiding walk of length 3 in the d-dimensional hypercubic lattice. - Michael Somos, Sep 06 2006
a(n) mod 10 is periodic 5: repeat [0, 2, 2, 6, 0]. - Mohamed Bouhamida, Sep 05 2009
This sequence is related to A005449 by a(n) = n*A005449(n)-sum(A005449(i), i=0..n-1), and this is the case d=3 in the identity n^2*(d*n+d-2)/2 - Sum_{k=0..n-1} k*(d*k+d-2)/2 = n*(n+d)*(2*d*n+d-3)/6. - Bruno Berselli, Nov 18 2010
Using (n, n+1) to generate a primitive Pythagorean triangle, the sides will be 2*n+1, 2*(n^2+n), and 2*n^2+2*n+1. Inscribing the largest rectangle with integral sides will have sides of length n and n^2+n. Side n is collinear to side 2*n+1 of the triangle and side n^2+n is collinear to side 2*(n^2+n) of the triangle. The areas of theses rectangles are a(n). - J. M. Bergot, Sep 22 2011
a(n+1) is the sum of n-th row of the triangle in A195437. - Reinhard Zumkeller, Nov 23 2011
Partial sums of A049450. - Omar E. Pol, Jan 12 2013
From Jon Perry, May 11 2013: (Start)
Define a 'stable brick triangle' as:
-----
| c |
---------
| a | | b |
----------
with a, b, c > 0 and c <= a + b. This can be visualized as two bricks with a third brick on top. The third brick can only be as strong as a+b, otherwise the wall collapses - for example, (1,2,4) is unstable.
a(n) gives the number of stable brick triangles that can be formed if the two supporting bricks are 1 <= a <= n and 1 <= b <= n: a(n) = Sum_{a=1..n} Sum_{b=1..n} Sum_c 1 = n^3 + n^2 as given in the Adamchuk formula.
So for i=j=n=2 we have 4:
1 2 3 4
2 2 2 2 2 2 2 2
For example, n=2 gives 2 from [a=1,b=1], 3 from both [a=1,b=2] and [a=2,b=1] and 4 from [a=2,b=2] so a(2) = 2 + 3 + 3 + 4 = 12. (End)
Define the infinite square array m(n,k) by m(n,k) = (n-k)^2 if n >= k >= 0 and by m(n,k) = (k+n)*(k-n) if 0 <= n <= k. This contains A120070 below the diagonal. Then a(n) = Sum_{k=0..n} m(n,k) + Sum_{r=0..n} m(r,n), the "hook sum" of the terms to the left of m(n,n) and above m(n,n) with irrelevant (vanishing) terms on the diagonal. - J. M. Bergot, Aug 16 2013
a(n) is the sum of all pairs with repetition drawn from the set of odd numbers 2*n-3. This is similar to A027480 but using the odd integers instead. Example using n=3 gives the odd numbers 1,3,5: 1+1, 1+3, 1+5, 3+3, 3+5,5+5 having a total of 36=a(3). - J. M. Bergot, Apr 05 2016
a(n) is the first Zagreb index of the complete graph K[n+1]. The first Zagreb index of a simple connected graph is the sum of the squared degrees of its vertices. Alternately, it is the sum of the degree sums d(i)+d(j) over all edges ij of the graph. - Emeric Deutsch, Nov 07 2016
a(n-2) is the maximum sigma irregularity over all trees with n vertices. The extremal graphs are stars. (The sigma irregularity of a graph is the sum of squares of the differences between the degrees over all edges of the graph.) - Allan Bickle, Jun 14 2023

Examples

			a(3) = 3^2+3^3 = 36.
		

References

  • L. B. W. Jolley, "Summation of Series", Dover Publications, 1961, pp. 50, 64.

Crossrefs

Cf. A011379, A181617, A270205 (sigma irregularities of maximal k-degenerate graphs).

Programs

Formula

a(n) = 2*A002411(n).
a(n) = Sum_{j=1..n} (Sum_{i=1..n} (i+j)), row sums of A126890 skipping numbers in the first column. - Alexander Adamchuk, Oct 12 2004
Sum_{n>0} 1/a(n) = (Pi^2 - 6)/6 = 0.6449340... [Jolley eq 272] - Gary W. Adamson, Dec 22 2006
a(n) = 2*n*binomial(n+1,2) = 2*n*A000217(n). - Arkadiusz Wesolowski, Feb 10 2012
G.f.: 2*x*(1 + 2*x)/(1 - x)^4. - Arkadiusz Wesolowski, Feb 11 2012
a(n) = A000330(n) + A002412(n) = A000292(n) + A002413(n). - Omar E. Pol, Jan 11 2013
a(n) = A245334(n+1,2), n > 0. - Reinhard Zumkeller, Aug 31 2014
Sum_{n>=1} 1/a(n) = A013661-1. - R. J. Mathar, Oct 18 2019 [corrected by Jason Yuen, Aug 04 2024]
Sum_{n>=1} (-1)^(n+1)/a(n) = 1 + Pi^2/12 - 2*log(2). - Amiram Eldar, Jul 04 2020
E.g.f.: exp(x)*x*(2 + 4*x + x^2). - Stefano Spezia, May 20 2021
a(n) = n*A002378(n) = A000578(n) + A000290(n). - J.S. Seneschal, Jun 18 2024

A028560 a(n) = n*(n + 6).

Original entry on oeis.org

0, 7, 16, 27, 40, 55, 72, 91, 112, 135, 160, 187, 216, 247, 280, 315, 352, 391, 432, 475, 520, 567, 616, 667, 720, 775, 832, 891, 952, 1015, 1080, 1147, 1216, 1287, 1360, 1435, 1512, 1591, 1672, 1755, 1840, 1927, 2016, 2107, 2200, 2295, 2392
Offset: 0

Views

Author

Keywords

Comments

Nonnegative X values of solutions to the equation X + (X + 3)^2 + (X + 6)^3 = Y^2. To prove that X = n^2 + 6n: Y^2 = X + (X + 3)^2 + (X + 6)^3 = X^3 + 19*X^2 + 115X + 225 = (X + 9)*(X^2 + 10X + 25) = (X + 9)*(X + 5)^2 it means: (X + 9) must be a perfect square, so X = k^2 - 9 with k>=3. we can put: k = n + 3, which gives: X = n^2 + 6n and Y = (n + 3)*(n^2 + 6n + 5). - Mohamed Bouhamida, Nov 12 2007
a(m) where m is a positive integer are the only positive integer values of t for which the Binet-de Moivre Formula of the recurrence b(n)=6*b(n-1)+t*b(n-2) with b(0)=0 and b(1)=1 has a root which is a square. In particular, sqrt(6^2+4*t) is an integer since 6^2+4*t=6^2+4*a(m)=(2*m+6)^2. Thus, the charcteristic roots are k1=6+m and k2=-m. - Felix P. Muga II, Mar 27 2014
Also, numbers k such that k + 9 is a perfect square.

Crossrefs

a(n-3), n>=4, third column (used for the Paschen series of the hydrogen atom) of triangle A120070.
Cf. A005563.

Programs

Formula

a(n) = (n+3)^2 - 3^2 = n*(n+6).
G.f.: x*(7-5*x)/(1-x)^3.
a(n) = 2*n + a(n-1) + 5. - Vincenzo Librandi, Aug 05 2010
Sum_{n>=1} 1/a(n) = 49/120 = 0.4083333... - R. J. Mathar, Mar 22 2011
a(n) = A028884(n) - 1. - Reinhard Zumkeller, Apr 07 2013
E.g.f.: x*(x+7)*exp(x). - G. C. Greubel, Aug 19 2017
Sum_{n>=1} (-1)^(n+1)/a(n) = 37/360. - Amiram Eldar, Nov 04 2020
a(n) = A056220(n+1) - A000290(n-1). - Leo Tavares, Sep 29 2022
From Amiram Eldar, Feb 05 2024: (Start)
Product_{n>=1} (1 - 1/a(n)) = -4*sqrt(10)*sin(sqrt(10)*Pi)/(3*Pi).
Product_{n>=1} (1 + 1/a(n)) = 45*sqrt(2)*sin(2*sqrt(2)*Pi)/(7*Pi). (End)

Extensions

Edited by Robert G. Wilson v, Feb 06 2002

A016061 a(n) = n*(n+1)*(4*n+5)/6.

Original entry on oeis.org

0, 3, 13, 34, 70, 125, 203, 308, 444, 615, 825, 1078, 1378, 1729, 2135, 2600, 3128, 3723, 4389, 5130, 5950, 6853, 7843, 8924, 10100, 11375, 12753, 14238, 15834, 17545, 19375, 21328, 23408, 25619, 27965, 30450, 33078, 35853, 38779, 41860
Offset: 0

Views

Author

Keywords

Comments

Number of ZnS molecules in cluster of n layers in zinc blende crystal.
(Zinc sulfide crystallizes in two different forms: wurtzite and zinc blende, the latter is also spelled zincblende.) - Jonathan Vos Post, Jan 22 2013
The Kn4 triangle sums of the Connell-Pol triangle A159797 lead to the sequence given above. For the definitions of the Kn4 and other triangle sums see A180662. - Johannes W. Meijer, May 20 2011
If one generated primitive Pythagorean triangles (2n+1, 2n+3) the collective sum of their perimeters for each n is four times the numbers listed in this sequence. - J. M. Bergot, Jul 18 2011
a(n) is the number of 3-tuples (w,x,y) having all terms in {0,...,n} and nA000292(n)+A000292(n+1)=n^3. - Clark Kimberling, Jun 04 2012
Degrees of the Hilbert polynomials for B_3 and C_3, per p. 13 of Gashi et al. - Jonathan Vos Post, Dec 14 2013
Number of solutions to a + b = c + d when 0 < a <= k, 0 <= b, c, d <= k, k = 0, 1, 2, 3.... Taken from Step 1 2007 problem #1(i) using 4 digit balanced numbers. - Bobby Milazzo, Mar 09 2013
From J. M. Bergot, Jun 18 2013: (Start)
Consider the lower half, including the main diagonal, of the array in A144216 as a triangle. The rows begin:
0;
1, 2;
3, 4, 6;
6, 7, 9, 12, ...
The sum of the terms in row(n) is a(n). (End)
This sequence is related to A008865 by a(n) = n*A008865(n+1) - Sum_{i=1..n} A008865(i) for n>0. - Bruno Berselli, Aug 06 2015

References

  • P. Jena and S. N. Behera, Clusters and Nanostructured Materials, Nova Science Publishers, 1996.

Crossrefs

Bisection of A002623.
Row sums of triangle A120070.

Programs

  • Magma
    I:=[0,3,13,34]; [n le 4 select I[n] else 4*Self(n-1)-6*Self(n-2)+4*Self(n-3)-Self(n-4): n in [1..40]]; // Vincenzo Librandi, Jul 25 2013
  • Maple
    A016061 := proc(n)
        n*(n+1)*(4*n+5)/6 ;
    end proc: # R. J. Mathar, Sep 26 2013
  • Mathematica
    CoefficientList[Series[x (3 + x) / (1 - x)^4, {x, 0, 40}], x] (* Vincenzo Librandi, Jul 25 2013 *)
    Table[n(n+1)(4*n+5)/6, {n,0,100}] (* Wesley Ivan Hurt, Sep 25 2013 *)
  • PARI
    v=vector(40,i,t(i)); s=0; forstep(i=2,40,2,s+=v[i]; print1(s","))
    

Formula

G.f.: x*(3+x)/(1-x)^4. - Paul Barry, Feb 27 2003
Partial sums of A014105. - Jon Perry, Jul 23 2003
a(n) = Sum_{i=0..n-1} 2*i^2 + i. - Jani Nurminen (slinky(AT)iki.fi), May 14 2006
a(n) = 2*n^3/3 +3*n^2/2 + 5*n/6. - Jonathan Vos Post, Dec 14 2013
a(n) = (4*n+5)/(2*n+1)*A000330(n). - Alexander R. Povolotsky, Mar 09 2013
a(n) = 4*a(n-1) -6*a(n-2) +4*a(n-3) -a(n-4). - Bobby Milazzo, Mar 10 2013
Sum_{n>=1} 1/a(n) = 12*Pi/5 + 72*log(2)/5 - 426/25. - Amiram Eldar, Jan 04 2022
E.g.f.: exp(x)*x*(18 + 21*x + 4*x^2)/6. - Stefano Spezia, Jul 31 2022

A024352 Numbers which are the difference of two positive squares, c^2 - b^2 with 1 <= b < c.

Original entry on oeis.org

3, 5, 7, 8, 9, 11, 12, 13, 15, 16, 17, 19, 20, 21, 23, 24, 25, 27, 28, 29, 31, 32, 33, 35, 36, 37, 39, 40, 41, 43, 44, 45, 47, 48, 49, 51, 52, 53, 55, 56, 57, 59, 60, 61, 63, 64, 65, 67, 68, 69, 71, 72, 73, 75, 76, 77, 79, 80, 81, 83, 84, 85, 87, 88, 89, 91, 92, 93, 95, 96
Offset: 1

Views

Author

Keywords

Comments

These are the solutions to the equation x^2 + xy = n where y mod 2 = 0, y is positive and x is any positive integer. - Andrew S. Plewe, Oct 19 2007
Ordered different terms of A120070 = 3, 8, 5, 15, 12, 7, ... (which contains two 15's, two 40's, and two 48's). Complement: A139544. (See A139491.) - Paul Curtz, Sep 01 2009
A024359(a(n)) > 0. - Reinhard Zumkeller, Nov 09 2012
If a(n) mod 6 = 3, n > 1, then a(n) = c^2 - f(a(n))^2 where f(n) = (floor(4*n/3) - 3 - n)/2. For example, 171 = 30^2 - 27^2 and f(171) = 27. - Gary Detlefs, Jul 15 2014

Crossrefs

Same as A042965 except for initial terms. - Michael Somos, Jun 08 2000
Different from A020884.

Programs

  • Haskell
    a024352 n = a024352_list !! (n-1)
    a024352_list = 3 : drop 4 a042965_list
    -- Reinhard Zumkeller, Nov 09 2012
    
  • Magma
    [3] cat [4 +Floor((4*n-3)/3): n in [2..100]]; // G. C. Greubel, Apr 22 2023
    
  • Mathematica
    Union[Flatten[Table[Select[Table[b^2 - c^2, {c, b-1}], # < 100 &], {b, 100}]]] (* Robert G. Wilson v, Jun 05 2004 *)
    LinearRecurrence[{1,0,1,-1},{3,5,7,8,9},70] (* Harvey P. Dale, Dec 20 2021 *)
  • PARI
    is(n)=(n%4!=2 && n>4) || n==3 \\ Charles R Greathouse IV, May 31 2013
    
  • Python
    def A024352(n): return 3 if n==1 else 3+(n<<2)//3 # Chai Wah Wu, Feb 10 2025
  • SageMath
    def A024352(n): return 4 + ((4*n-3)//3) - int(n==1)
    [A024352(n) for n in range(1,101)] # G. C. Greubel, Apr 22 2023
    

Formula

Consists of all positive integers except 1, 4 and numbers == 2 (mod 4).
a(n) = a(n-3) + 4, n > 4.
G.f.: (3 + 2*x + 2*x^2 - 2*x^3 - x^4)/(1 - x - x^3 + x^4). - Ralf Stephan, before May 13 2008
a(n) = a(n-1) + a(n-3) - a(n-4), for n > 5. - Ant King, Oct 03 2011
a(n) = 4 + floor((4*n-3)/3), n > 1. - Gary Detlefs, Jul 15 2014

Extensions

Edited by N. J. A. Sloane, Sep 19 2008

A028566 a(n) = n*(n+8).

Original entry on oeis.org

0, 9, 20, 33, 48, 65, 84, 105, 128, 153, 180, 209, 240, 273, 308, 345, 384, 425, 468, 513, 560, 609, 660, 713, 768, 825, 884, 945, 1008, 1073, 1140, 1209, 1280, 1353, 1428, 1505, 1584, 1665, 1748, 1833, 1920, 2009, 2100, 2193, 2288, 2385
Offset: 0

Views

Author

Keywords

Comments

a(m) where m is a positive integer are the only positive integer values of t for which the Binet-de Moivre Formula of the recurrence b(n) = 8*b(n-1) + t*b(n-2) with b(0) = 0 and b(1) = 1 has a root which is a square. In particular, sqrt(8^2 + 4*t) is a positive integer since 8^2 + 4*t = 8^2 + 4*a(m) = (2*m + 8)^2. Thus, the characteristics roots are r1 = 8 + m and r2 = -m. - Felix P. Muga II, Mar 28 2014

Crossrefs

a(n - 4), n >= 5, fourth column (used for the Brackett series of the hydrogen atom) of triangle A120070.

Programs

  • Magma
    [n*(n+8): n in [0..50]]; // G. C. Greubel, Jul 31 2022
  • Mathematica
    Table[n (n+8), {n, 0, 50}] (* Bruno Berselli, Apr 06 2014 *)
  • PARI
    a(n)=n*(n+8)
    
  • Sage
    [n*(n+8) for n in [0..50]] # Bruno Berselli, Apr 06 2014
    
  • Scala
    (0 to 49).map { n: Int => n * n + 8 * n } // Alonso del Arte, Nov 10 2019
    

Formula

a(n) = (n+4)^2 - 4^2 = n*(n+8), n >= 0.
G.f.: x*(9 - 7*x)/(1 - x)^3.
a(n) = 2*n + a(n-1) + 7. - Vincenzo Librandi, Aug 05 2010
Sum_{n >= 1} 1/a(n) = 761/2240 = 0.3397321... - R. J. Mathar, Mar 22 2011
Sum_{n>=1} (-1)^(n+1)/a(n) = 533/6720. - Amiram Eldar, Jan 15 2021
E.g.f.: x*(9 + x)*exp(x). - G. C. Greubel, Jul 31 2022
From Amiram Eldar, Feb 12 2024: (Start)
Product_{n>=1} (1 - 1/a(n)) = 315*sin(sqrt(17)*Pi)/(13*sqrt(17)*Pi).
Product_{n>=1} (1 + 1/a(n)) = -32*sqrt(15)*sin(sqrt(15)*Pi)/(11*Pi). (End)

A098603 a(n) = n*(n+10).

Original entry on oeis.org

0, 11, 24, 39, 56, 75, 96, 119, 144, 171, 200, 231, 264, 299, 336, 375, 416, 459, 504, 551, 600, 651, 704, 759, 816, 875, 936, 999, 1064, 1131, 1200, 1271, 1344, 1419, 1496, 1575, 1656, 1739, 1824, 1911, 2000, 2091, 2184, 2279, 2376, 2475, 2576, 2679, 2784
Offset: 0

Views

Author

Eugene McDonnell (eemcd(AT)mac.com), Nov 04 2004

Keywords

Comments

These are the only positive integer values of t for which the Binet-de Moivre formula for the recurrence b(n) = 10*b(n-1)+t*b(n-2) with b(0)=0 and b(1)=1 has a root which is a square. In particular, sqrt(10^2+4*t) is a positive integer since 10^2+4*t = 10^2+4*a(m) = (2*m+10)^2. Thus the characteristic roots are r1=10+m and r2 = -m. - Felix P. Muga II, Mar 28 2014

Crossrefs

Cf. A098832.
a(n-5), n>=6, fifth column (used for the Pfund series of the hydrogen atom) of triangle A120070.

Programs

Formula

a(n) = (n+5)^2 - 5^2 = n*(n+10), n>=0.
G.f.: x*(11-9*x)/(1-x)^3.
a(n) = a(n-1) + 2*n + 9, (with a(0)=0). - Vincenzo Librandi, Nov 17 2010
Sum_{n>=1} 1/a(n) = 7381/25200 via sum_{n>=0} 1/((n+x)*(n+y)) = (psi(x)-psi(y))/(x-y). - R. J. Mathar, Jul 14 2012
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3), with a(0)=0, a(1)=11, a(2)=24. - Harvey P. Dale, Jul 26 2014
Sum_{n>=1} (-1)^(n+1)/a(n) = 1627/25200. - Amiram Eldar, Jan 15 2021
E.g.f.: x*(11 + x)*exp(x). - G. C. Greubel, Jul 31 2022
From Amiram Eldar, Feb 12 2024: (Start)
Product_{n>=1} (1 - 1/a(n)) = -18144*sqrt(2/13)*sin(sqrt(26)*Pi)/(935*Pi).
Product_{n>=1} (1 + 1/a(n)) = 126*sqrt(6)*sin(2*sqrt(6)*Pi)/(23*Pi). (End)

Extensions

More terms from Emeric Deutsch, Mar 11 2005

A098849 a(n) = n*(n + 16).

Original entry on oeis.org

0, 17, 36, 57, 80, 105, 132, 161, 192, 225, 260, 297, 336, 377, 420, 465, 512, 561, 612, 665, 720, 777, 836, 897, 960, 1025, 1092, 1161, 1232, 1305, 1380, 1457, 1536, 1617, 1700, 1785, 1872, 1961, 2052, 2145, 2240, 2337, 2436, 2537, 2640, 2745, 2852, 2961
Offset: 0

Views

Author

Eugene McDonnell (eemcd(AT)mac.com), Nov 04 2004

Keywords

Crossrefs

a(n-8), n>=9, eighth column (used for the n=8 series of the hydrogen atom) of triangle A120070.

Programs

Formula

a(n) = (n+8)^2 - 8^2 = n*(n + 16), n>=0.
G.f.: x*(17 - 15*x)/(1-x)^3.
a(n) = a(n-1) + 2*n + 15 (with a(0)=0). - Vincenzo Librandi, Nov 17 2010
From G. C. Greubel, Jul 29 2016: (Start)
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3).
E.g.f.: x*(17 + x)*exp(x). (End)
From Amiram Eldar, Jan 15 2021: (Start)
Sum_{n>=1} 1/a(n) = H(16)/16 = A001008(16)/A102928(16) = 2436559/11531520, where H(k) is the k-th harmonic number.
Sum_{n>=1} (-1)^(n+1)/a(n) = 95549/2306304. (End)

Extensions

More terms from Emeric Deutsch, Mar 26 2005

A098847 a(n) = n*(n + 12).

Original entry on oeis.org

0, 13, 28, 45, 64, 85, 108, 133, 160, 189, 220, 253, 288, 325, 364, 405, 448, 493, 540, 589, 640, 693, 748, 805, 864, 925, 988, 1053, 1120, 1189, 1260, 1333, 1408, 1485, 1564, 1645, 1728, 1813, 1900, 1989, 2080, 2173, 2268, 2365, 2464, 2565, 2668, 2773
Offset: 0

Views

Author

Eugene McDonnell (eemcd(AT)mac.com), Nov 04 2004

Keywords

Crossrefs

a(n-6), n>=7, sixth column (used for the n=6 series of the hydrogen atom) of triangle A120070.

Programs

Formula

a(n) = (n+6)^2 - 6^2 = n*(n + 12), n>=0.
G.f.: x*(13 - 11*x)/(1-x)^3.
a(n) = 2*n + a(n-1) + 11 (with a(0)=0). - Vincenzo Librandi, Nov 17 2010
a(0)=0, a(1)=13, a(2)=28, a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3). - Harvey P. Dale, May 24 2012
Sum_{n>=1} 1/a(n) = 86021/332640 = 0.258600... via Sum_{n>=0} 1/((n+x)(n+y)) = (psi(x)-psi(y))/(x-y). - R. J. Mathar, Jul 14 2012
E.g.f.: x*(13 + x)*exp(x). - G. C. Greubel, Jul 29 2016
Sum_{n>=1} (-1)^(n+1)/a(n) = 18107/332640. - Amiram Eldar, Jan 15 2021

Extensions

More terms from Robert G. Wilson v, Jul 14 2005

A098848 a(n) = n*(n + 14).

Original entry on oeis.org

0, 15, 32, 51, 72, 95, 120, 147, 176, 207, 240, 275, 312, 351, 392, 435, 480, 527, 576, 627, 680, 735, 792, 851, 912, 975, 1040, 1107, 1176, 1247, 1320, 1395, 1472, 1551, 1632, 1715, 1800, 1887, 1976, 2067, 2160, 2255, 2352, 2451, 2552, 2655, 2760, 2867
Offset: 0

Views

Author

Eugene McDonnell (eemcd(AT)mac.com), Nov 04 2004

Keywords

Crossrefs

Cf. A098832.
a(n-7), n>=8, seventh column (used for the n=7 series of the hydrogen atom) of triangle A120070.

Programs

Formula

a(n) = (n+7)^2 - 7^2 = n*(n + 14), n>=0.
G.f.: x*(15 - 13*x)/(1-x)^3.
a(n) = 2*n + a(n-1) + 13 (with a(0)=0). - Vincenzo Librandi, Nov 16 2010
Sum_{n>=1} 1/a(n) = 1171733/5045040 = 0.2322544518... via Sum_{n>=0} 1/((n+x)(n+y)) = (psi(x)-psi(y))/(x-y). - R. J. Mathar, Jul 14 2012
From G. C. Greubel, Jul 29 2016: (Start)
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3).
E.g.f.: x*(15 + x)*exp(x). (End)
Sum_{n>=1} (-1)^(n+1)/a(n) = 237371/5045040. - Amiram Eldar, Jan 15 2021

Extensions

More terms from Robert G. Wilson v, Jul 14 2005

A120072 Numerator triangle for hydrogen spectrum rationals.

Original entry on oeis.org

3, 8, 5, 15, 3, 7, 24, 21, 16, 9, 35, 2, 1, 5, 11, 48, 45, 40, 33, 24, 13, 63, 15, 55, 3, 39, 7, 15, 80, 77, 8, 65, 56, 5, 32, 17, 99, 6, 91, 21, 3, 4, 51, 9, 19, 120, 117, 112, 105, 96, 85, 72, 57, 40, 21, 143, 35, 5, 1, 119, 1, 95, 5, 7, 11, 23
Offset: 2

Views

Author

Wolfdieter Lang, Jul 20 2006

Keywords

Comments

Frequencies or energies of the spectral lines of the hydrogen (H) atom are given, according to quantum theory, by r(m,n)*3.287*PHz (1 Peta Hertz= 10^15 s^{-1}) or r(m,n)*13.599 eV (electron Volts), respectively. The wave lengths are lambda(m,n) = (1/r(m,n))* 91.196 nm (all decimals rounded). See the W. Lang link for more details.
The spectral series for n=1,2,...,7, m>=n+1, are named after Lyman, Balmer, Paschen, Brackett, Pfund, Humphreys, Hansen-Strong, respectively.
The corresponding denominator triangle is A120073.
The rationals are r(m,n):= a(m,n)/A120073(m,n) = A120070(m,n)/(m^2*n^2) = 1/ n^2 - 1/m^2 and they are given in lowest terms.

Examples

			For the rational triangle see W. Lang link.
Numerator triangle begins as:
   3;
   8,  5;
  15,  3,  7;
  24, 21, 16,  9;
  35,  2,  1,  5, 11;
  48, 45, 40, 33, 24, 13;
  63, 15, 55,  3, 39,  7, 15;
  80, 77,  8, 65, 56,  5, 32, 17;
  99,  6, 91, 21,  3,  4, 51,  9, 19;
		

Crossrefs

Row sums give A120074.
Row sums of r(m, n) triangle give A120076(m)/A120077(m), m>=2.

Programs

  • Magma
    [Numerator(1/k^2 - 1/n^2): k in [1..n-1], n in [2..18]]; // G. C. Greubel, Apr 24 2023
    
  • Mathematica
    Table[1/n^2 - 1/m^2, {m,2,12}, {n,m-1}]//Flatten//Numerator (* Jean-François Alcover, Sep 16 2013 *)
  • SageMath
    def A120072(n,k): return numerator(1/k^2 - 1/n^2)
    flatten([[A120072(n,k) for k in range(1,n)] for n in range(2,19)]) # G. C. Greubel, Apr 24 2023

Formula

a(m,n) = numerator(r(m,n)) with r(m,n) = 1/n^2 - 1/m^2, m>=2, n=1..m-1.
The g.f.s for the columns n=1,..,10 of triangle r(m,n) = a(m, n) / A120073(m, n), m >= 2, 1 <= n <= m-1, are given in the W. Lang link.
Previous Showing 21-30 of 48 results. Next