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

A349264 Generalized Euler numbers, a(n) = n!*[x^n](sec(4*x)*(sin(4*x) + 1)).

Original entry on oeis.org

1, 4, 16, 128, 1280, 16384, 249856, 4456448, 90767360, 2080374784, 52975108096, 1483911200768, 45344872202240, 1501108249821184, 53515555843342336, 2044143848640217088, 83285910482761809920, 3605459138582973251584, 165262072909347030040576, 7995891855149741436305408
Offset: 0

Views

Author

Peter Luschny, Nov 20 2021

Keywords

Examples

			Exponential generating functions of generalized Euler numbers in context:
egf1 = sec(1*x)*(sin(x) + 1).
   [A000111, A000364, A000182]
egf2 = sec(2*x)*(sin(x) + cos(x)).
   [A001586, A000281, A000464]
egf3 = sec(3*x)*(sin(2*x) + cos(x)).
   [A007289, A000436, A000191]
egf4 = sec(4*x)*(sin(4*x) + 1).
   [A349264, A000490, A000318]
egf5 = sec(5*x)*(sin(x) + sin(3*x) + cos(2*x) + cos(4*x)).
   [A349265, A000187, A000320]
egf6 = sec(6*x)*(sin(x) + sin(5*x) + cos(x) + cos(5*x)).
   [A001587, A000192, A000411]
egf7 = sec(7*x)*(-sin(2*x) + sin(4*x) + sin(6*x) + cos(x) + cos(3*x) - cos(5*x)).
   [A349266, A064068, A064072]
egf8 = sec(8*x)*2*(sin(4*x) + cos(4*x)).
   [A349267, A064069, A064073]
egf9 = sec(9*x)*(4*sin(3*x) + 2)*cos(3*x)^2.
   [A349268, A064070, A064074]
		

Crossrefs

Programs

  • Maple
    sec(4*x)*(sin(4*x) + 1): series(%, x, 20): seq(n!*coeff(%, x, n), n = 0..19);
  • Mathematica
    m = 19; CoefficientList[Series[Sec[4*x] * (Sin[4*x] + 1), {x, 0, m}], x] * Range[0, m]! (* Amiram Eldar, Nov 20 2021 *)
  • PARI
    seq(n)={my(x='x + O('x^(n+1))); Vec(serlaplace((sin(4*x) + 1)/cos(4*x)))} \\ Andrew Howroyd, Nov 20 2021

A002437 a(n) = A000364(n) * (3^(2*n+1) + 1)/4.

Original entry on oeis.org

1, 7, 305, 33367, 6815585, 2237423527, 1077270776465, 715153093789687, 626055764653322945, 698774745485355051847, 968553361387420436695025, 1632180870878422847476890007, 3286322019402928956112227932705, 7791592461957309952817483706344167, 21485762937086358457367440231243675985
Offset: 0

Views

Author

Keywords

Comments

The terms are multiples of the Euler numbers (A000364).

Examples

			a(4) = A000364(4) * (3^(2*4+1)+1)/4 = 1385 * (3^9+1)/4 = 1385 * 4921 = 6815585.
		

References

  • A. Fletcher, J. C. P. Miller, L. Rosenhead and L. J. Comrie, An Index of Mathematical Tables. Vols. 1 and 2, 2nd ed., Blackwell, Oxford and Addison-Wesley, Reading, MA, 1962, Vol. 1, p. 75.
  • J. W. L. Glaisher, Messenger of Math., 28 (1898), 36-79, see esp. p. 51.
  • L. B. W. Jolley, Summation of Series, Dover, 2nd ed. (1961)
  • 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

Bisections: A156168, A156169.
Cf. other sequences with a g.f. of the form cos(x)/(1 - k*sin^2(x)): A012494 (k=-1), A001209 (k=1/2), A000364 (k=1), A000281 (k=2), A156134 (k=3).

Programs

  • Maple
    Q:=proc(n) option remember; if n=0 then RETURN(1); else RETURN(expand((u^2+1)*diff(Q(n-1),u)+u*Q(n-1))); fi; end;
    [seq(subs(u=sqrt(3),Q(2*n)),n=0..25)];
  • Mathematica
    Table[Abs[EulerE[2 n]] (3^(2 n + 1) + 1) / 4, {n, 0, 30}] (* Vincenzo Librandi, Feb 07 2017 *)

Formula

A000364(n) * (3^(2*n+1) + 1)/4.
Q_2n(sqrt(3)), where the polynomials Q_n() are defined in A104035. - N. J. A. Sloane, Nov 06 2009
a(n) = (-1)^n*Sum_{k = 0..2*n-1} w^(2*n+k)*Sum_{j = 1..2*n-1} (-1)^(k-j)*binomial(2*n-1,k-j)*(2*j - 1)^(2*n-2), where w = exp(2*Pi*i/6) (i = sqrt(-1)). Cf. A002439. - Peter Bala, Jan 21 2011
Sum_{n>=1} (-1)^floor((n-1)/2) 1/A007310(n)^s = r_s with r_{2s+1} = 2 *(Pi/6)^(2s+1) *a(s) /(2s)!. [Jolley eq (315)]. - R. J. Mathar, Mar 24 2011
From Peter Bala, Feb 06 2017: (Start)
E.g.f.: cos(x)^2/cos(3*x) = cos(x)/(1 - 4*sin(x)^2) = 1 + 7*x^2/2! + 305*x^4/4! + 33367*x^6/6! + .... This is the even part of (1/2)*sec(x + Pi/3). Cf. A000191. (End)
a(n) = (1/2)*Integral_{x = 0..inf} x^(2*n)*cosh(Pi*x/3)/cosh(Pi*x/2) dx. - Cf. A000281. - Peter Bala, Nov 08 2019

Extensions

More terms from Herman P. Robinson
Further terms from N. J. A. Sloane, Nov 06 2009

A352977 Expansion of e.g.f. cos(2x) cos(3x) / cos(6x) (even powers only).

Original entry on oeis.org

1, 23, 3985, 1743623, 1424614945, 1870693029623, 3602792061891505, 9566946196183630823, 33500193836861731481665, 149565522713623779723211223, 829235405016410370201483113425, 5589623533324449496004527793434823, 45017811997394066193946619670380594785
Offset: 0

Views

Author

F. Chapoton, Apr 13 2022

Keywords

Comments

Only terms of even index are given. Terms of odd index are zero.

Crossrefs

Intermediate case between A002437 and A349429.
Cf. A000192.

Programs

  • Maple
    egf := (cos(x) + cos(5*x))*sec(6*x) / 2: ser := series(egf, x, 32):
    seq(n!*coeff(ser, x ,n), n = 0..24, 2); # Peter Luschny, Apr 13 2022
  • PARI
    my(x='x+O('x^30)); select(x->(x>0), Vec(serlaplace(cos(2*x)*cos(3*x)/cos(6*x)))) \\ Michel Marcus, Apr 13 2022
  • Sage
    x = PowerSeriesRing(QQ, 'x', default_prec=30).gen()
    f = cos(2*x) * cos(3*x) / cos(6*x)
    [cf for cf in f.egf_to_ogf() if cf]
    

Formula

E.g.f.: cos(2*x) * cos(3*x) / cos(6*x).
From Peter Luschny, Apr 13 2022: (Start)
E.g.f.: (cos(x) + cos(5*x))*sec(6*x) / 2, even powers only.
a(n) = A000192(n)/2. (End)
a(n) ~ 2^(6*n + 3/2) * 3^(2*n + 1/2) * n^(2*n + 1/2) / (Pi^(2*n + 1/2) * exp(2*n)). - Vaclav Kotesovec, Apr 15 2022
Showing 1-3 of 3 results.