OFFSET
0,3
COMMENTS
a(46) is negative. - Vaclav Kotesovec, Oct 03 2020
LINKS
Seiichi Manyama, Table of n, a(n) for n = 0..400 (terms n = 0..100 from Vaclav Kotesovec)
FORMULA
E.g.f.: exp( Sum_{k>0} sigma(k)*sin(x)^k/k ).
MATHEMATICA
nmax = 25; CoefficientList[Series[Product[1/(1 - Sin[x]^k), {k, 1, nmax}], {x, 0, nmax}], x] * Range[0, nmax]! (* Vaclav Kotesovec, Oct 03 2020 *)
PROG
(PARI) N=40; x='x+O('x^N); Vec(serlaplace(1/eta(sin(x))))
(PARI) N=40; x='x+O('x^N); Vec(serlaplace(1/prod(k=1, N, 1-sin(x)^k)))
(PARI) N=40; x='x+O('x^N); Vec(serlaplace(exp(sum(k=1, N, sigma(k)*sin(x)^k/k))))
CROSSREFS
KEYWORD
sign
AUTHOR
Seiichi Manyama, Oct 02 2020
STATUS
approved