login
The OEIS Foundation is supported by donations from users of the OEIS and by a grant from the Simons Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A318007 E.g.f. A(x) satisfies: A(x) = sin(x) + cos(x)*A(x)^2 with A(0)=0. 4
1, 2, 11, 100, 1321, 22622, 474851, 11786920, 337650001, 10962556442, 397804232891, 15954966065740, 700861670953081, 33464274335643062, 1725656338796874131, 95578727117480612560, 5658893822397686566561, 356659432612090891074482, 23841281202421071709150571, 1684762749472172141605523380 (list; graph; refs; listen; history; text; internal format)
OFFSET

1,2

LINKS

Robert Israel, Table of n, a(n) for n = 1..368

FORMULA

E.g.f. A(x) satisfies:

(1) A(x) = sin(x) + cos(x)*A(x)^2.

(2) A(x) = sin(x) * Sum_{n>=0} binomial(2*n,n)/(n+1) * sin(2*x)^n/2^n.

(3) A(x) = (1 - sqrt(1 - 2*sin(2*x))) / (2*cos(x)).

(4) A(x) = 2*sin(x) / (1 + sqrt(1 - 2*sin(2*x))).

a(n) ~ (sqrt(3) - 1) * 2^(2*n - 3/2) * 3^(n - 1/4) * n^(n-1) / (Pi^(n - 1/2) * exp(n)). - Vaclav Kotesovec, Oct 13 2020

EXAMPLE

E.g.f.: A(x) = x + 2*x^2/2! + 11*x^3/3! + 100*x^4/4! + 1321*x^5/5! + 22622*x^6/6! + 474851*x^7/7! + 11786920*x^8/8! + 337650001*x^9/9! + 10962556442*x^10/10! + ...

such that A(x) = sin(x) + cos(x)*A(x)^2.

MAPLE

E:= (1 - sqrt(1 - 2*sin(2*x))) / (2*cos(x)):

S:= series(E, x, 31):

seq(coeff(S, x, j)*j!, j=1..30); # Robert Israel, Aug 29 2018

MATHEMATICA

m = 21; A[x_] = (1 - Sqrt[1 - 2 Sin[2 x]] )/(2 Cos[x]); Rest[Range[0, m - 1]! * CoefficientList[A[x] + O[x]^m, x]] (* Jean-Fran├žois Alcover, Apr 29 2019 *)

PROG

(PARI) {a(n) = my(A = 2*sin(x +x^2*O(x^n)) / (1 + sqrt(1 - 2*sin(2*x +x^2*O(x^n)))) ); n!*polcoeff(A, n)}

for(n=1, 25, print1(a(n), ", "))

CROSSREFS

Cf. A318003, A318004, A318599.

Sequence in context: A143135 A205806 A220433 * A243950 A056732 A157715

Adjacent sequences:  A318004 A318005 A318006 * A318008 A318009 A318010

KEYWORD

nonn

AUTHOR

Paul D. Hanna, Aug 28 2018

STATUS

approved

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

License Agreements, Terms of Use, Privacy Policy. .

Last modified January 25 07:24 EST 2021. Contains 340416 sequences. (Running on oeis4.)