login
Triangle read by rows: coefficients in the expansion cos(x)^n = (1/2)^n * Sum_{k=0..n} T(n,k) * cos(k*x).
16

%I #35 Aug 27 2018 10:15:53

%S 1,0,2,2,0,2,0,6,0,2,6,0,8,0,2,0,20,0,10,0,2,20,0,30,0,12,0,2,0,70,0,

%T 42,0,14,0,2,70,0,112,0,56,0,16,0,2,0,252,0,168,0,72,0,18,0,2,252,0,

%U 420,0,240,0,90,0,20,0,2

%N Triangle read by rows: coefficients in the expansion cos(x)^n = (1/2)^n * Sum_{k=0..n} T(n,k) * cos(k*x).

%C These coefficients are especially useful when integrating powers of cosine x (see examples).

%C Nonzero, even elements of the first column are given by A000984; T(2n,0) = binomial(2n,n).

%C For the rational triangles for even and odd powers of cos(x) see A273167/A273168 and A244420/A244421, respectively. - _Wolfdieter Lang_, Jun 13 2016

%C Mathematica needs no TrigReduce to integrate Cos[x]^k. See link. - _Zak Seidov_, Jun 13 2016

%H Zak Seidov, <a href="/A273496/a273496.pdf">No Need For TrigReduce</a>

%F From _Robert Israel_, May 24 2016: (Start)

%F T(n,k) = 0 if n-k is odd.

%F T(n,0) = binomial(n,n/2) if n is even.

%F T(n,k) = 2*binomial(n,(n-k)/2) otherwise. (End)

%e n/k| 0 1 2 3 4 5 6

%e -------------------------------

%e 0 | 1

%e 1 | 0 2

%e 2 | 2 0 2

%e 3 | 0 6 0 2

%e 4 | 6 0 8 0 2

%e 5 | 0 20 0 10 0 2

%e 6 | 20 0 30 0 12 0 2

%e -------------------------------

%e cos(x)^4 = (1/2)^4 (6 + 8 cos(2x) + 2 cos(4x)).

%e I4 = Int dx cos(x)^4 = (1/2)^4 Int dx ( 6 + 8 cos(2x) + 2 cos(4x) ) = C + 3/8 x + 1/4 sin(2x) + 1/32 sin(4x).

%e Over range [0,2Pi], I4 = (3/4) Pi.

%t T[MaxN_] := Function[{n}, With[

%t {exp = Expand[Times[ 2^n, TrigReduce[Cos[x]^n]]]},

%t Prepend[Coefficient[exp, Cos[# x]] & /@ Range[1, n],

%t exp /. {Cos[_] -> 0}]]][#] & /@ Range[0, MaxN];Flatten@T[10]

%t (* alternate program *)

%t T2[MaxN_] := Function[{n}, With[{exp = Expand[(Exp[I x] + Exp[-I x])^n]}, Prepend[2 Coefficient[exp, Exp[I # x]] & /@ Range[1, n], exp /. {Exp[_] -> 0}]]][#] & /@ Range[0, MaxN]; T2[10] // ColumnForm (* _Bradley Klee_, Jun 13 2016 *)

%Y Cf. A000984, A001790, A046161, A038533, A038534, A273506, A273507, A273167, A273168, A244420, A244421.

%K nonn,tabl

%O 0,3

%A _Bradley Klee_, May 23 2016