login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A225109 E.g.f. (sin(3x) + cos x) / cos(4x). 4
1, 3, 15, 117, 1185, 15123, 230895, 4116837, 83860545, 1921996323, 48942778575, 1370953667157, 41893214676705, 1386843017916723, 49441928730798255, 1888542637550347077, 76946148390480577665, 3331009898404800736323, 152682246738275154625935, 7387240827905368219116597 (list; graph; refs; listen; history; text; internal format)
OFFSET
0,2
LINKS
FORMULA
a(n) = Im(2*((1-I)/(1+I))^n*(1+sum_{j=0..n}(binomial(n,j)*Li_{-j}(I)* 4^j))). - Peter Luschny, Apr 29 2013
a(n) ~ n! * sqrt(2+sqrt(2)) * 2^(3*n+1)/Pi^(n+1). - Vaclav Kotesovec, Jun 02 2013
MAPLE
per4 := proc(n) local j; 2*((1-I)/(1+I))^n*(1+add(binomial(n, j)* polylog(-j, I)*4^j, j=0..n)) end: A225109 := n -> Im(per4(n));
seq(A225109(i), i=0..11); # Peter Luschny, Apr 29 2013
MATHEMATICA
mx = 17; Range[0, mx]! CoefficientList[ Series[ (Sin[3x] + Cos[x])/Cos[4x], {x, 0, mx}], x] (* Robert G. Wilson v, Apr 28 2013 *)
PROG
(PARI) v=Vec((sin(3*x) + cos(x)) / cos(4*x)); vector(#v, i, v[i]*(i-1)!)
(PARI) x='x+O('x^66); Vec(serlaplace((sin(3*x)+cos(x))/cos(4*x))) \\ Joerg Arndt, Apr 28 2013
CROSSREFS
Sequence in context: A074596 A087806 A245043 * A080290 A365777 A372161
KEYWORD
nonn
AUTHOR
M. F. Hasler, Apr 28 2013
STATUS
approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified April 22 01:34 EDT 2024. Contains 371887 sequences. (Running on oeis4.)