Ich habe ein Spezialverfahren entwickelt, wie du cos ^ n als Fourierreihe ausdrücken kannst. Hier kennst du den Spruch
   " Der kürzeste Umweg zur reellen Analysis führt immer noch über die komplexe Ebene. "
   Eulersatz; setze
    a  :=  exp  (  +  i  x  )  ;  b  :=  exp  (  -  i  x  )        (  1a  )
    
     Beachte insbesondere
        a  b  =  1     (  1b  )
    Dann konstituiert der Eulersatz das  LGS
      a  =  cos  (  x  )  +  i  sin  (  x  )           (  2a  )
      b  =  cos  (  x  )  -  i  sin  (  x  )          (  2b  )
   Rein formal juristisch lassen sich ( 2ab ) nach den beiden Unbekannten cos und sin auflösen:
          cos  (  x  )  =  1/2  (  a  +  b  )       (  3  )
   Alles was du jetzt noch brauchst, ist der binomische Lehrsatz; beachte jedoch insbesondere ( 1b ) .   Da jedoch hinterher etwas Reelles raus kommen soll, musst du die Summe so umsortieren, dass jeweils komplex konjugierte Pärchen nebeneinander stehen:
   cos  ^ n  (  x  )  =  ( 1 / 2 ^ n )    SUMME  ( n k )   [ a ^ ( n - 2 k ) + b ^ ( n - 2 k ) ]   ( 4a )
   Für gerade n erstreckt sich die binomische Summe wegen der Symmetrie nur bis  k = n/2 ,  für ungerade n bis k = ( n - 1 ) / 2 , so dass du bekommst
   cos  ^ n  (  x  )  =  [ 1 / 2 ^ ( n - 1 ) ]  SUMME  ( n  k  )  cos  (  n  -  2  k  )  x    (  4b  )
   als Fourierreihe der Potenzen der Winkelfunktionen.  Diese Oberwellen lassen sich aber elementar integrieren.