OFFSET
0,2
LINKS
G. C. Greubel, Table of n, a(n) for n = 0..385
FORMULA
E.g.f.: 1/(cos(2*x) - sin(2*x)).
a(n) = 2^n * A001586(n).
a(n) = | 2*8^n*lerchphi(-1,-n,1/4) |. - Peter Luschny, Apr 27 2013
G.f.: 1/Q(0), where Q(k) = 1 - 2*x*(2*k+1) - 8*x^2*(k+1)^2/Q(k+1) ; (continued fraction). - Sergei N. Gladkovskii, Sep 27 2013
a(n) ~ 4 * n^(n+1/2) * (8/Pi)^n / (sqrt(Pi)*exp(n)). - Vaclav Kotesovec, Oct 07 2013
E.g.f.: 1/(1-2*x)*(1 + 2*x^2/((1-2*x)*W(0) - x )), where W(k) = x + (k+1)/( 1 - 2*x/( 2*k+3 - x*(2*k+3)/W(k+1) )); (continued fraction ). - Sergei N. Gladkovskii, Dec 27 2013
MATHEMATICA
CoefficientList[Series[1/(Cos[2*x]-Sin[2*x]), {x, 0, 20}], x]* Range[0, 20]! (* Vaclav Kotesovec, Oct 07 2013 *)
PROG
(PARI) a(n)=if(n<0, 0, n!*polcoeff(1/(cos(2*x+x*O(x^n))-sin(2*x+x*O(x^n))), n))
(Sage)
from mpmath import mp, lerchphi
mp.dps = 32; mp.pretty = True
def A079858(n): return abs(2*8^n*lerchphi(-1, -n, 1/4))
[int(A079858(n)) for n in (0..17)] # Peter Luschny, Apr 27 2013
(PARI) x='x+O('x^66); v=Vec(serlaplace( 1/(cos(2*x)-sin(2*x)) ) ) \\ Joerg Arndt, Apr 27 2013
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Michael Somos, Jan 20 2003
STATUS
approved