OFFSET
0,4
LINKS
Seiichi Manyama, Table of n, a(n) for n = 0..800
FORMULA
a(n+4) = 2*A249059(n) for n >= 0.
E.g.f.: exp(x*(2+x)/2) * (exp(1/2) * sqrt(2*Pi) * (erf((1+x)/sqrt(2)) - erf(1/sqrt(2))) - 1). - Vaclav Kotesovec, Sep 05 2017
a(n) ~ (sqrt(Pi) * exp(1/2) * (1 - erf(1/sqrt(2))) - sqrt(2)/2) * n^(n/2) * exp(sqrt(n) - n/2 - 1/4). - Vaclav Kotesovec, Sep 05 2017
MATHEMATICA
RecurrenceTable[{a[0] == -1, a[1] == 1, a[n] == a[n-1] + (n-1)*a[n-2]}, a[n], {n, 0, 20}] (* Vaclav Kotesovec, Sep 04 2017 *)
CoefficientList[Series[E^(x*(2 + x)/2) * (E^(1/2)*Sqrt[2*Pi]*(Erf[(1 + x)/Sqrt[2]] - Erf[1/Sqrt[2]]) - 1), {x, 0, 20}], x]*Range[0, 20]! (* Vaclav Kotesovec, Sep 05 2017 *)
PROG
(GAP)
a:=[-1, 1];; for n in [3..10^2] do a[n]:=a[n-1]+(n-2)*a[n-2]; od; a; # Muniru A Asiru, Sep 07 2017
CROSSREFS
KEYWORD
sign
AUTHOR
Seiichi Manyama, Sep 04 2017
STATUS
approved