[go: up one dir, main page]

login
Revisions by Johannes W. Meijer (See also Johannes W. Meijer's wiki page)

(Bold, blue-underlined text is an addition; faded, red-underlined text is a deletion.)

Showing entries 1-10 | older changes
Eulerian numbers (Euler's triangle: column k=2 of A008292, column k=1 of A173018).
(history; published version)
#385 by Johannes W. Meijer at Mon Jun 03 16:48:21 EDT 2019
STATUS

editing

proposed

#384 by Johannes W. Meijer at Mon Jun 03 16:47:51 EDT 2019
COMMENTS

The utility of large networks, particularly social networks, with n participants is given by the terms a(n) of this sequence. This assertion is known as Reed’s Law, see the Wikipedia link. - Johannes W. Meijer, Jun 03 2019

LINKS

Wikipedia, <a href="https://en.wikipedia.org/wiki/Reed%27s_law">Reed’s Law</a>

STATUS

approved

editing

Gijswijt's sequence: a(1) = 1; for n>1, a(n) = largest integer k such that the word a(1)a(2)...a(n-1) is of the form xy^k for words x and y (where y has positive length), i.e., the maximal number of repeating blocks at the end of the sequence so far.
(history; published version)
#91 by Johannes W. Meijer at Sun Feb 03 10:08:46 EST 2019
STATUS

editing

proposed

#90 by Johannes W. Meijer at Sun Feb 03 10:08:10 EST 2019
LINKS

D. C. Gijswijt, <a href="httphttps://www.pyth.eu/jaargangenuploads/user/ArchiefPDF/Pyth55-3.pdf">Krulgetallen</a>, Pythagoras, 55ste Jaargang, Nummer 3, Jan 2016, (Article in Dutch about this sequence, see pages 10-13, cover and back cover).

STATUS

approved

editing

Discussion
Sun Feb 03
10:08
Johannes W. Meijer: Link repaired.
Decimal expansion of the higher-order exponential integral E(x, m=2, n=1) at x=1.
(history; published version)
#30 by Johannes W. Meijer at Fri Sep 28 10:08:55 EDT 2018
STATUS

editing

proposed

#29 by Johannes W. Meijer at Fri Sep 28 10:08:16 EDT 2018
COMMENTS

We define the higher-order exponential integrals by E(x,m,n) = x^(n-1)*int(E(t,m-1,n)/t^n, t=x..infinity) for m => 0 =1 and n => =1 with E(x,m=0,n) = exp(-x), see Meijer and Baken.

FORMULA

E(x=1,m=2,n=1) = gamma^2/2 + Pi^2/12 + Sum_{k=1..infinity} ((-1)^k/(k^2*k!)).

E(x=0,n,m) = (1/(n-1))^m for n=>=2.

E(x,m,n) = (-1)^m * ((-x)^(n-1)/(n-1)!) *sum Sum_{kz=0..floor(m/2)}(alpha (kz, n)*(G(m-2*kz, n) ) + (-1) ^m * ((-x)^(n-1)/(n-1)!) * Sum_{kz=0..floor(m/2)} (Sum_{i=1..m-2*kz}(alpha (kz, n) *G(m-2*kz-i, n)*logln(x)^i/i!,i=)) + (-1..)^m-2 *kz))) + Sum_{ kx=0..(n-2)} ((-x)^kx/((kx-n+1)^m*kx!)) + (-1)^m * Sum_{ky=n..infinity}((-x)^ky /(( ky-n+1)^m*ky!))).

AUTHOR

Johannes W. Meijer & and _Nico Baken (n.h.g.baken(AT)tudelft.nl), _, Aug 13 2009, Aug 17 2009

STATUS

approved

editing

Discussion
Fri Sep 28
10:08
Johannes W. Meijer: Formula repaired, see style sheet, plus some minor edits.
#14 by Johannes W. Meijer at Fri Sep 28 10:05:07 EDT 2018
STATUS

editing

proposed

#13 by Johannes W. Meijer at Fri Sep 28 10:04:44 EDT 2018
COMMENTS

The higher order exponential integrals, see A163931, are defined by E(x,m,n) = x^(n-1)*int(E(t,m-1,n)/t^n, t=x..infinity) for m =>0 =1 and n>=1, with E(x,m=0,n) = exp(-x).

FORMULA

G(k,n) = (1/k)*((gamma-sum(p^*G(k-1,n)) - (1/k),*Sum_{p=1..n-1}(p^(-1))* G(k-1,n) +sum (1/k)* Sum_{i=0..k-2}(Zeta(k-i) -sum(p^* G(i,n)) - (1/k-)*Sum_{i)),=0..k-2}(Sum_{p=1..n-1}(p^(i-k)) * G(i,n),i=0..k-2)) with G(0,n) = 1 for k => =0 and n => =1.

GF(z,n) = GAMMA(n-z)/GAMMA(n).

CROSSREFS

(gamma - G(1,n)) equals A001008(n-1)/A002805(n-1) for n =>=2.

The structure of the G(k,n=1) formula leads formulas lead (replace gamma by G and Zeta by Z) to A036039.

AUTHOR

Johannes W. Meijer & and _Nico Baken (n.h.g.baken(AT)tudelft.nl), _, Aug 13 2009, Aug 17 2009

STATUS

approved

editing

Discussion
Fri Sep 28
10:05
Johannes W. Meijer: Formula rewritten, see the style sheet, plus some minor edits.
Numerators of the higher order exponential integral constants alpha(k,4).
(history; published version)
#10 by Johannes W. Meijer at Fri Sep 28 10:00:45 EDT 2018
STATUS

editing

proposed

#9 by Johannes W. Meijer at Fri Sep 28 09:59:05 EDT 2018
COMMENTS

The higher order exponential integrals, see A163931, are defined by E(x,m,n) = x^(n-1)*int(E(t,m-1,n)/t^n, t=x..infinity) for m >=0 1 and n>=1, with E(x,m=0,n) = exp(-x).

The first Maple program uses the alpha(k,n) formula and the second the GF(z,n) to generate the alpha(k,n) coefficients in each column.

FORMULA

alpha(k,n) = (1/k)*sum(sum(p^(-2 *( Sum_{i=0..k-i)),1} (Sum_{p = 0..n-1}(p^(2*i-2*k))*alpha(i, n), i = 0..k-1) with alpha(0,n) = 1, k >= 0 and n >= 1.

GF(z,n) = product((1-(z/k)^2)^(-1), k = 1..n-1) = (Pi*z/sin(Pi*z))/(Beta(n+z,n-z)/Beta(n,n))

MAPLE

coln:=4; nmax1 := nmax+1: 16; for n from 0 to nmax1 do A008955(n, 0):=1 end do: for n from 0 to nmax1 do A008955(n, n) := (n!)^2 end do: for n from 1 to nmax1 do for m from 1 to n-1 do A008955(n, m) := A008955(n-1, m-1)*n^2 + A008955(n-1, m) end do: end do: m:=coln-1: f(m):=0: for n from 0 to m do f(m) := f(m) + (-1)^(n + m)*A008955(m, n)*z^(2*m-2*n) od: GF(z, coln) := m!^2/f(m): GF(z, coln):=series(GF(z, coln), z, nmaxnmax1);

CROSSREFS

Cf. A163931 (E(x, m, n)), A163930 (gamma(k, n)).

The GF(z,n) lead to A008955.

The denominators of a(1,n), n >= 2, lead to A007407.

AUTHOR

Johannes W. Meijer & and _Nico Baken (n.h.g.baken(AT)tudelft.nl), _, Aug 13 2009, Aug 17 2009

EXTENSIONS

Edited by Johannes W. Meijer, Sep 22 2012

STATUS

approved

editing

Discussion
Fri Sep 28
10:00
Johannes W. Meijer: Formula rewritten, see the style sheet, error corrected, some minor edits and link added.