積分 文字式 フーリエ級数展開 python

Question

pythonを使ってフーリエ級数展開の答えを求めようとしています。
フーリエ級数展開では以下の式を解かなければいけません。
an=1/π*∫(π −π)f(x)cosnx dx (n=0,1,2,...)
その時に文字式のcosnxをそのまま積分する方法が分かりません。

from scipy.integrate import quad
import numpy as np

def fourier(fun):
    a = []
    b = []
    res,err = quad(lambda x:fun(x)*cos(n*x),-pi,pi)
    a.append(res/pi)
    res,err = quad(lambda x:fun(x)*sin(n*x),-pi,pi)
    b.append(res/pi)
    return a,b

と大体で書いてみたのですが、うまくいきませんでした。
解決策やコードが分かる方助けてください。

Answer 1

from pylab import *
を追加する。cos, pi にnp.を付ける必要がなくなる。

Answer 2

だいたいあってると思います。
cosやpiはnp.をつけるとよいのではないでしょうか。

余談ですが、MathematicaやSymPy使わないのですか