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

 Hints (Greetings from The On-Line Encyclopedia of Integer Sequences!)
 A012091 cos(arcsin(arctan(x)))=1-1/2!*x^2+5/4!*x^4-109/6!*x^6+4521/8!*x^8... 1

%I #11 Feb 07 2015 16:58:07

%S 1,-1,5,-109,4521,-330681,36468717,-5721101541,1205289090513,

%T -328897528901361,112769846051251797,-47472161654530776285,

%U 24068680135210975714425,-14466929567049590650593705

%N cos(arcsin(arctan(x)))=1-1/2!*x^2+5/4!*x^4-109/6!*x^6+4521/8!*x^8...

%F a(n) = (2*n)!*sum(k=1..2*n, ((2*k)!*binomial(2*k-2,k-1)*(-1)^(n+k+1)*sum(i=2*k..2*n, (2^(i+1)*stirling1(i,2*k)*binomial(2*n-1,i-1))/i!))/(k*2^(4*k))) with n>0, a(0)=1. [_Vladimir Kruchinin_, Oct 08 2012]

%t With[{nn=30},Take[CoefficientList[Series[Cos[ArcSin[ArcTan[x]]],{x,0,nn}],x] Range[0,nn]!,{1,-1,2}]] (* _Harvey P. Dale_, Feb 07 2015 *)

%o (Maxima) a[n]:=if n=0 then 1 else (2*n)!*sum(((2*k)!*binomial(2*k-2,k-1)*(-1)^(n+k+1)*sum((2^(i+1)*stirling1(i,2*k)*binomial(2*n-1,i-1))/i!, i, 2*k, 2*n))/(k*2^(4*k)), k, 1, 2*n); makelist(a[n], n, 0, 13); [_Vladimir Kruchinin_, Oct 08 2012]

%Y Cf. A012254.

%K sign

%O 0,3

%A Patrick Demichel (patrick.demichel(AT)hp.com)

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.

Last modified June 10 14:31 EDT 2023. Contains 363205 sequences. (Running on oeis4.)