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-2 of 2 results.

A256500 Decimal expansion of the positive solution to x = 2*(1-exp(-x)).

Original entry on oeis.org

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

Views

Author

Stanislav Sykora, Mar 31 2015

Keywords

Comments

Each of the positive solutions to x = q*(1-exp(-x)) obtained for q = 2, 3, 4, and 5, appears in several formulas pertinent to Planck's black-body radiation law. For a given q, the solution can be also written as q+W(-q/exp(q)), where W is the Lambert function. Here q = 2.
The constant appears in asymptotic formula for A007820. - Vladimir Reshetnikov, Oct 10 2016

Examples

			1.5936242600400400923230418758751602417890024248188593649995...
		

Crossrefs

Cf. A194567 (q=3), A256501 (q=4), A256502 (q=5).

Programs

  • Mathematica
    RealDigits[2 + LambertW[-2 Exp[-2]], 10, 100][[1]] (* Vladimir Reshetnikov, Oct 10 2016 *)
  • PARI
    a2=solve(x=0.1,10,x-2*(1-exp(-x))) \\ Use real precision in excess

Formula

Equals 2*(1-A106533). - Miko Labalan, Dec 18 2024
Equals log(A229553). - Hugo Pfoertner, Dec 19 2024

A226762 Greatest k such that 1/k >= mean of {1, 1/2, 1/3, ..., 1/n}.

Original entry on oeis.org

1, 1, 1, 1, 2, 2, 2, 2, 3, 3, 3, 3, 4, 4, 4, 4, 4, 5, 5, 5, 5, 5, 6, 6, 6, 6, 6, 7, 7, 7, 7, 7, 8, 8, 8, 8, 8, 8, 9, 9, 9, 9, 9, 10, 10, 10, 10, 10, 10, 11, 11, 11, 11, 11, 11, 12, 12, 12, 12, 12, 12, 13, 13, 13, 13, 13, 13, 14, 14, 14, 14, 14, 14, 15, 15
Offset: 1

Views

Author

Clark Kimberling, Jun 19 2013

Keywords

Comments

Largest integer not exceeding the harmonic mean of the first n numbers.

Examples

			1/4 < mean{1,1/2,1/3,...,1/9} < 1/3, so that a(9) = 3.
		

Crossrefs

Programs

  • Mathematica
    f[n_] := Mean[Table[1/k, {k, 1, n}]]
    Table[Floor[1/f[n]], {n, 1, 120}]   (* this sequence *)
    Table[Ceiling[1/f[n]], {n, 1, 120}] (* A226763 *)
  • PARI
    \\ This uses only precision-independent integer operations:
    a(n)=(n*n!)\sum(k=1,n,n!\k)  \\ Stanislav Sykora, Apr 08 2015

Formula

a(n) = floor(n/(Sum_{k=1..n} 1/k)).

Extensions

Name corrected by Jason Yuen, Nov 02 2024
Showing 1-2 of 2 results.