A012279 Expansion of e.g.f. sec(exp(x)*log(x+1)).

1, 0, 1, 3, 16, 70, 580, 3850, 39940, 350544, 4408460, 48605788, 715662176, 9518173080, 160469679136, 2503106438040, 47524013744272, 851125166412928, 17967763294589776, 363412179220061872

Offset

0,4

Links

G. C. Greubel, Table of n, a(n) for n = 0..441

Example

E.g.f. = 1 + x^2/2! + 3*x^3/3! + 16*x^4/4! + 70*x^5/5! + ...

Maple

seq(coeff(series(factorial(n)*sec(exp(x)*log(x+1)), x, n+1), x, n), n = 0 .. 20); # Muniru A Asiru, Oct 28 2018

Mathematica

With[{nn = 30}, CoefficientList[Series[Sec[Exp[x]*Log[x+1]], {x, 0, nn}], x] Range[0, nn]!] (* G. C. Greubel, Oct 28 2018 *)

Prog

(PARI) x='x+O('x^30); Vec(serlaplace(1/cos(exp(x)*log(x+1)))) \\ G. C. Greubel, Oct 28 2018
(Magma) m:=30; R<x>:=PowerSeriesRing(Rationals(), m); b:=Coefficients(R!( 1/Cos(Exp(x)*Log(x+1)) )); [Factorial(n-1)*b[n]: n in [1..m-1]]; // G. C. Greubel, Oct 28 2018

Crossrefs

Sequence in context: A370248 A370274 A015524 * A037098 A316170 A038602

Adjacent sequences: A012276 A012277 A012278 * A012280 A012281 A012282

Keyword

nonn

Author

Patrick Demichel (patrick.demichel(AT)hp.com)

Status

approved